BZOJ 3135: [Baltic2013]pipes

来源:互联网 发布:tcp 51 udp 端口 编辑:程序博客网 时间:2024/06/07 00:53

首先如果m>n 未知数就比方程多了,所以肯定是无解或者多解 自然就输出0了

这样子就把问题变成了 n-1<=m<=n
treeDP就好 一个是树 一个是环套树

嗯 我觉得我的代码还是挺好看的
代码如果有误请指出(因为我跟网上的一个AC代码对拍的时候有数据答案不一样 虽说那组我对了),谢谢!

#include<bits/stdc++.h>using namespace std;const int N=100002,inf=1e9+7;inline int read(){    char ch=getchar(); int x=0,f=1;    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0'; ch=getchar();}    return x*f;}struct node{int y,next,id;}e[N<<1]; int len,first[N];void ins(int x,int y,int d){    e[++len]=(node){y,first[x],d},first[x]=len;}int n,m,s[N],w[N],c[N],cl,b[N],t;bool in[N],v[N];int dfs(int x,int fa){    for(int k=first[x];k;k=e[k].next){        int y=e[k].y; if(y==fa || in[y])continue;        s[x]-= (w[e[k].id]=dfs(y,x));    }    return s[x];}bool fc(int x,int p){    v[x]=1;    for(int k=first[x];k;k=e[k].next){        if(k==(p^1))continue; int y=e[k].y;        if(!v[y]){            if(fc(y,k)){                c[++cl]=x,b[cl]=e[k].id,in[x]=1;                if(x!=t)return 1;                return 0;            }        }        else{            t=y,c[++cl]=x,b[cl]=e[k].id,in[x]=1;            return 1;        }    }    return 0;}int main(){    n=read(),m=read(); len=1;    if(n<m) return printf("0\n"),0;    for(int i=1;i<=n;++i)s[i]=read();    for(int i=1;i<=m;++i){        int x=read(),y=read();        ins(x,y,i),ins(y,x,i);    }    if(n==m){        fc(1,0);        if(!(cl&1)) return printf("0\n"),0;        for(int i=1;i<=cl;++i)dfs(c[i],0);        int x=0;        for(int i=1;i<=cl;++i)            x+=((i&1)?1:-1)*s[c[i]];        if(x&1) return printf("0\n"),0;        x=x>>1; w[b[1]]=x;        for(int i=2;i<=cl;++i)w[b[i]]=s[c[i-1]]-w[b[i-1]];        for(int i=1;i<=m;++i)printf("%d\n",w[i]*2);    }    else{        if(dfs(1,0))printf("0\n");        else            for(int i=1;i<=m;++i)printf("%d\n",w[i]*2);    }    return 0;}
原创粉丝点击