Poj 1926 Pollution

//这道题目很简单,矩阵相乘即可,只是在精度上要控制好

//虽然简单却暴露了自己是多么马虎,以此为戒吧 :(

//1.首先是sizeof问题(代码中有标明),2.然后是看题,比如说这组数据:

1
3 2
1
2
0.6
2 3
1 2
它是不能平衡的,,但是题目说一定会保证input数据都能平衡

3.最后就是致命的,,忘了孤立点的事情!,,蛋蛋疼!

#include <iostream> #include <queue> #include <stack> #include <string> #include <map> #include <vector> #include <cstring> #include <cstdio> #include <cmath> #include <algorithm> using namespace std; const int MAX=105,INF=1<<30; int n,m,di; bool g[MAX][MAX]; double tank[MAX],mat[MAX][MAX],odd[MAX][MAX], I[MAX][MAX]; int d[MAX]; inline void mul (double m1[MAX][MAX],double m2[MAX][MAX]){ double ans[MAX][MAX]; for(int i=0;i<n;i++){ for(int j=0;j<n;j++){ ans[i][j]=0; for(int k=0;k<n;k++){ ans[i][j]+=m1[i][k]*m2[k][j]; } } } memcpy(m1,ans,sizeof(double)*MAX*MAX);//不能sizeof(m1), //因为它在函数里面只是一个指针 return ; } void Pow(double a[MAX][MAX],int k){ memcpy(odd,I,sizeof(odd)); while(k){ if(k&1) mul(odd,a); mul(a,a); k>>=1; } return ; } int main() { #ifndef ONLINE_JUDGE freopen("i.txt", "r", stdin); //freopen("e://t1.txt", "w", stdout); #endif int T,u,v; double ans[MAX]; cin>>T; memset(I,0,sizeof(I)); for(int i=0;i<MAX;i++) I[i][i]=1.0; while(T--){ memset(g,0,sizeof(g)); memset(d,0,sizeof(d)); memset(mat,0,sizeof(mat)); cin>>n>>m; for(int i=0;i<n;i++) cin>>tank[i]; for(int i=0;i<m;i++){ cin>>u>>v; --u; --v; g[u][v]=1; g[v][u]=1; } for(int i=0;i<n;i++){ di=0; for(int j=0;j<n;j++){ if(g[i][j]) ++di; } d[i]=di; if(!di) continue; for(int j=0;j<n;j++){ if(g[i][j]){ mat[i][j]=1/(double)di; } } } Pow(mat,2*m*m); memset(ans,0,sizeof(ans)); for(int j=0;j<n;j++){ for(int k=0;k<n;k++){ ans[j]+=tank[k]*odd[k][j]; } } for(int i=0;i<n;i++){ if(d[i]) printf("%.3lf/n",ans[i]); else printf("%.3lf/n",tank[i]); } printf("/n"); } return 0; }