网友投稿 528 2022-10-30
迭代法解线性方程组
#include #include #include #include#include#include#include #include typedef long long ll;typedef unsigned long long ull;using namespace std;typedef pair pii;#define#define#define#define#defineconst int maxn=1e7+1;#define#define#define#define#defineconst int mod=998244353;const int MOD=10007;ll n;double a[100][100];double b[maxn],x_k[maxn],x_k1[maxn];#defineint main(){ cin>>n; for(int i=1;i<=n;i++){ for(int j=1;j<=n;j++){ cin>>a[i][j]; } } for(int i=1;i<=n;i++) { cin>>b[i]; x_k[i]=x_k1[i]=0; } int i,j,k=0; double sum,epsilon; ll flag=0; while(1) { epsilon=0.0; k=k+1; for(int i=1;i<=n;i++){ sum=0; for(int j=1;j<=n;j++){ if(j!=i) sum+=a[i][j]*x_k[j]; } x_k1[i]=(b[i]-sum)*1.0/a[i][i]; epsilon+=fabs(x_k1[i]-x_k[i]); } for(int i=1;i<=n;i++) x_k[i]=x_k1[i]; if(epsilon < EPS) { flag=1; break; } if(k>=MAX){ flag=0; break; } } if(flag){ for(int i=1;i<=n;i++) cout< 版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。
版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。
发表评论
暂时没有评论，来抢沙发吧~
暂时没有评论，来抢沙发吧~