bzoj3135: [Baltic2013]pipesd

Description

有n个水库，m条管道。Jester会在某些管道中间凿开一个洞，让水流出来或者用水泵把水打进去。保证这个流速是偶数。对于一条管道(u, v)，如果在中间凿开了一个洞让水流出来，流速是2d m^3/s，那么水库u和v失水速度为d m^3/s。同理，如果往一条管道(u, v)注水，流速为2p m^3/s，那么u和v得到水的速度是p m^3/s。

给定图的构造以及每个水库的水流的变化，问每条边的方案是否唯一。

Input

第一行：水库数量n，管道数量m （1 <= n <= 100 000, 1 <= m <= 500 000）下面n行：第i个水库的变化速度ci （-10^9 <= ci <= 10^9）接下来m行：(u, v)，保证没有重边

Output

如果方案唯一，输出方案，每行一个数xi（-10^9 <= xi <= 10^9）表示第i条管道的流量变化。放水为负数，灌水为正数。否则输出0。

边看做未知数则每个点构成一个方程，保证有解，因此n<m时方案不唯一

对于一个度为1的点，与其相邻的边的答案可以被确定，因此可以以此为突破口消元，若n=m-1则最后可以全部解出，若n=m则最后剩下一个环

对于偶环，显然方案不唯一，对于奇环则可以简单地推出公式O(n)解出唯一解

#include<cstdio>
typedef long long i64;
char buf[],*ptr=buf,*pmx=buf+;
inline int g(){
    if(ptr==pmx)fread(ptr=buf,,,stdin);
    return *(ptr++);
}
int _(){
    int x=,c=g(),f=;
    while(c<)c=='-'&&(f=-),c=g();
    while(c>)x=x*+c-,c=g();
    return x*f;
}
int n,m;
int et[],enx[],e0[],ep=,deg[],q[],ql=,qr=,v[],e[],p=;
i64 f[],ans[],s[];
void dfs(int w,int pa){
    deg[w]=-;
    v[++p]=w;
    for(int i=e0[w];i;i=enx[i]){
        int u=et[i];
        if(deg[u]==)e[p]=i>>,dfs(u,w);
        else if(u!=pa&°[u]==-)e[p]=i>>;
    }
}
int main(){
    n=_();m=_();
    if(m>n)return puts(""),;
    for(int i=;i<=n;++i)f[i]=_();
    for(int i=,a,b;i<m;++i){
        a=_();b=_();
        ++deg[a];++deg[b];
        et[ep]=b;enx[ep]=e0[a];e0[a]=ep++;
        et[ep]=a;enx[ep]=e0[b];e0[b]=ep++;
    }
    for(int i=;i<=n;++i)if(deg[i]==)q[++qr]=i;
    while(ql!=qr){
        int w=q[++ql];
        deg[w]=;
        for(int i=e0[w];i;i=enx[i]){
            int u=et[i];
            if(!deg[u])continue;
            ans[i>>]=f[w]*;
            f[u]-=f[w];
            if(==--deg[u])q[++qr]=u;
        }
    }
    for(int w=;w<=n;++w)if(deg[w]==){
        dfs(w,);
        if(p%==)return puts(""),;
        v[p+]=v[];
        for(int i=;i<=p;++i)s[i]=s[i-]+(i&?f[v[i+]]:-f[v[i+]]);
        for(int i=;i<=p;++i)ans[e[i]]=(s[p]-s[i-]*)*(i&?:-);
        break;
    }
    for(int i=;i<=m;++i)printf("%lld\n",ans[i]);
    return ;
}