1580. Dean's Debts

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;
}