高斯消元

线性代数的方法,利用矩阵来求解多元一次线性方程组的根。

https://blog.youkuaiyun.com/just_sort/article/details/60327014

众所周知,柱爷的数学非常好,尤其擅长概率论!

某日柱爷在喵蛤蛤村散步,无意间踏入了远古法阵!

法阵很奇怪,是一个长度为N

的走廊,初始时柱爷在最左边,现在柱爷要到最右边去!

柱爷的行动方式如下:

每个回合柱爷会投一次骰子,根据骰子上的点数
1
X,柱爷会相应的往右边移动X步.

骰子的数值是
1到6,取到每面的概率相同

在某些位置可能有传送门,一旦柱爷在该回合结束后在这个位置上,会被强制传送到传送门的另外一边

传送门是单向的,同时每个位置不会有超过1个传送门,同时不会存在a→b,b→c这种情况

在任意时刻柱爷都必须保证在法阵内,也就说如果在这一回合结束后柱爷的位置在法阵外,那么这回合柱爷将什么都不做

那么请问柱爷到达最右边的期望回合数是多少呢?或者是永远都无法到达?
Input

第一行两个整数N
,M

,分别表示法阵的长度和传送门的数量

接下来M
行,每行两个整数u,v,表示从u到v

有一扇传送门

数据保证:

1≤N≤300
1
0≤M≤[N−22]

1

 

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define maxn 305
#define esp (1e-14)
using namespace std;
int n,m,f[maxn];
long double a[maxn][maxn];//构造的高斯消元的矩阵,代表第i个方程式的第j个系数是多少 ,精度要求很高 
 
void read(int& x){
    x=0;
    char c=getchar();
    for(;c<'0'||c>'9';c=getchar());
    for(;c>='0' && c <= '9'; c = getchar()  ){
        x = x*10 + c - '0';
    }
}
 
 
int main(){
    read( n );read(m);//读入优化 
    int x,y;
    for(int i=1;i <= m ; i++){
        read(x);read(y);
        f [ x ] = y;//如果有传送的话,到哪里 
    }
    
    for(int i = 1 ; i < n ; i++ ){
        a[ i ][ i ]=6;//第一个方程  
        if( f[ i ] )a[ i ][ f[ i ] ]=-6.0;//如果有传送门 系数直接抵消 x-y=0 相当于 x=y 
        else{
            a[i][n+1]=6.0;//方程右边的常数 
            for(int j = 1 ;j <= 6; j++){
                if(i+j>n)a[i][i]-=1.0;//另外一个方程 
                else a[i][i+j]-=1.0;
            }
        }
    }
    a[n][n]=1.0;//最后的方程 
    a[n][n+1]=0;
    //高斯消元 
    for(int i=1;i<=n;i++){
        int p=i;
        for(int j=i+1;j<=n;j++){//向下查找第j个系数不为0的方程 
            if(fabs(a[j][i])> esp)p=j;
        }
        if(fabs(a[p][i])>esp){
            for(int j=i;j<=n+1;j++)swap(a[i][j],a[p][j]);//把方程移上来 
            for(int j=i+1;j<=n;j++){//向下消元 同时除去其他的系数 
                if(fabs(a[j][i])>esp){
                    long double k=a[j][i]/a[i][i];//消元 
                    for(int t=i;t<=n+1;t++)a[j][t]-=a[i][t]*k;//系数相减 
                }
            }
        }
    }  
    //回带                 a[i][n+1]就是第i个未知数的解 
    for(int i=n;i>=1;i--){
        for(int j = i+1;j <= n; j++)if(fabs(a[i][j])>esp){
            a[i][n+1]-=a[j][n+1]*a[i][j];//用已知的解求未知解 
        }
        if(abs(a[i][i])<=esp&&abs(a[i][n+1])>esp){//如果出现矛盾 
            printf("-1");
            return 0;
        }
        a[i][1+n]/=a[i][i];//求出当前的解 
    }
    
    printf("%.12f",(double)a[1][1+n]);//输出第一个未知数的解 就是答案 
    return 0;
}

 

 

https://blog.youkuaiyun.com/pengwill97/article/details/77282328

 

转载于:https://www.cnblogs.com/downrainsun/p/10980105.html

内容概要:本文介绍了DeepSeek在职场中的应用,从提示语技巧到多场景应用,涵盖了DeepSeek的基础模型(V3)、深度思考模型(R1)及其联网搜索功能。文中详细描述了DeepSeek的模型对比,包括操作规范、结果导向、路径灵活性、响应模式和风险特征等方面。此外,还探讨了DeepSeek在制作可视化图表、PPT、海报、视频以及批量生成新媒体文案等具体应用场景中的使用方法和技巧。最后,文章展示了DeepSeek在市场调查、AI应用开发等方面的应用实例,强调了其在人机协同和共生领域的潜力。 适用人群:适用于希望提升工作效率和创新能力的职场人士,特别是从事数据分析、内容创作、市场营销、AI开发等领域的专业人士。 使用场景及目标:①通过DeepSeek的基础模型(V3)和深度思考模型(R1)进行高效的任务处理和复杂推理;②利用DeepSeek制作可视化图表、PPT、海报和视频,提高内容创作的质量和效率;③通过DeepSeek进行市场调查和AI应用开发,优化业务流程并推动创新。 其他说明:DeepSeek不仅提供了强大的AI工具,还强调了人机协同的重要性。用户在使用过程中应注意操作规范,结合实际需求选择合适的模型,并充分利用DeepSeek的各项功能来实现高效的职场应用。文中还提到了多个国际竞赛中的获奖情况,展示了DeepSeek团队在AI领域的卓越实力。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值