http://acm.timus.ru/problem.aspx?space=1&num=1580
简单题 对于每个联通块进行搜索 必须所以连通块都可以求出解 才可以
对于某个联通块 可以设起点为 a[]+b[]*x (0+1*x ) 然后不断搜索 搜到没有搜到的
计算这个点的a[]+b[]*x 如果搜到已经搜过的 如果b[i]+b[j]等于0 则看a[i]+a[j]是否符合所给的边 不符合则无解
否则求出 x 的值 如果一个连通块中多次求得 x 的值 则每次必须一样
代码:
#include<iostream> #include<cstdio> #include<cstring> #include<string> #include<map> #include<vector> #include<stack> #include<set> #include<map> #include<queue> #include<algorithm> #include<cmath> #define LL long long #define sint short int //#pragma comment(linker, "/STACK:1024000000,1024000000") using namespace std; const int N=1005; const double ex=1e-6; const int INF=0x3f3f3f3f; int head[N],I; struct node { int j,next; int c; }edge[N*N]; bool visited[N]; int a[N],b[N]; double c[N]; void add(int i,int j,int c) { edge[I].j=j; edge[I].c=c; edge[I].next=head[i]; head[i]=I++; } bool equals(double x1,double x2) { if(abs(x1-x2)<ex) return true; return false; } bool bfs(int l) { queue<int>qt; stack<int>st; a[l]=0; b[l]=1; qt.push(l); st.push(l); visited[l]=true; double X; bool flag=false; while(!qt.empty()) { int i=qt.front(); qt.pop(); for(int t=head[i];t!=-1;t=edge[t].next) { int j=edge[t].j; if(!visited[j]) { visited[j]=true; a[j]=edge[t].c-a[i]; b[j]=-b[i]; qt.push(j); st.push(j); }else { if(b[i]+b[j]==0) { if(a[i]+a[j]!=edge[t].c) {return false;} }else { double tmpx=1.0*(edge[t].c-(a[i]+a[j]))/(b[i]+b[j]); if(flag&&!equals(tmpx,X)) {return false;} flag=true;X=tmpx; } } } } if(flag==false) return false; while(!st.empty()) { int k=st.top(); st.pop(); c[k]=a[k]+X*b[k]; } return true; } int main() { //freopen("data.in","r",stdin); int n,m; while(cin>>n>>m) { memset(head,-1,sizeof(head)); I=0; while(m--) { int l,r; int k; cin>>l>>r>>k; add(l,r,k); add(r,l,k); } memset(visited,false,sizeof(visited)); int l; for(l=1;l<=n;++l) { if(!visited[l]) { if(!bfs(l)) break; } } if(l<=n) cout<<"IMPOSSIBLE"<<endl; else { for(int i=1;i<=n;++i) printf("%.2lf\n",c[i]); } } return 0; }