BZOJ 2597: [Wc2007]剪刀石头布|费用流

神题吓死人- -！费用流神题！！！
n个点考虑最多有C（n，3）个题意中描述的3元环
考虑每个人赢的次数
赢一次对这种3圆环的破坏个数为0，两次为1，三次再加2，四次再加3……
ans=C（n,3）-破坏个数
显然要让破坏个数最少，那就是费用流的模型了
原图中的每条边抽象为点
S连每条边流量为1费用0
每条边若胜负不定（2） 连两个人流量为1费用为0
若胜负确定（0,1）连胜利的人流量为1费用为0
然后每个人连到汇点n-1条边（其实这样显然多了，可以在前面计算一下连边的数目）流量都为1，费用依次为0,1,2,3,……

#include<cstdio>#include<cstdlib>#include<cstring>#include<queue>#include<cmath>#include<vector>#include<set>#include<map>#include<iostream>#include<algorithm>#define MX 1e9#define T 111#define E 11111using namespace std;int sc(){    int i=0; char c=getchar();    while(c>'9'||c<'0')c=getchar();    while(c>='0'&&c<='9')i=i*10+c-'0',c=getchar();    return i;}int dis[E],vis[E],q[E],pre[E],from[E];int head[E],lst[E*100],nxt[E*100],v[E*100],c[E*100];int s[T][T],d[T];int n,m,tot=1,num,ans;int S,W;void insert(int x,int y,int _c,int _v){    lst[++tot]=y; nxt[tot]=head[x]; c[tot]=_c; v[tot]=_v; head[x]=tot;    lst[++tot]=x; nxt[tot]=head[y]; c[tot]=0; v[tot]=-_v; head[y]=tot;}bool spfa(){    for(int i=1;i<=W;i++) dis[i]=MX; dis[S]=0;    int l=1,r=2; q[1]=S;    while(l!=r)    {        int x=q[l++]; l%=E; vis[x]=0;        for(int i=head[x];i;i=nxt[i])            if(c[i]&&dis[lst[i]]>dis[x]+v[i])            {                dis[lst[i]]=dis[x]+v[i];                from[lst[i]]=x;                pre[lst[i]]=i;                if(!vis[lst[i]])                {                    vis[lst[i]]=1;                    q[r++]=lst[i];                    r%=E;                }            }    }    return dis[W]!=MX;}void mcf(){    int cc=MX;    for(int i=W;i!=S;i=from[i])cc=min(cc,c[pre[i]]);    for(int i=W;i!=S;i=from[i])    {        int k=pre[i];        ans+=v[k]*cc;        c[k]-=cc;        c[k^1]+=cc;    }}int find(int x){    for(int i=head[x];i;i=nxt[i])        if(!c[i])return lst[i];}int main(){    n=sc();m=n*(n-1)/2;    S=n+m+1,W=S+1;num=n;    for(int i=1;i<=n;i++)        for(int j=1;j<=n;j++)        {            s[i][j]=sc();            if(i>=j)continue;            num++; insert(S,num,1,0);            if(s[i][j]==1)insert(num,i,1,0),d[i]++;            else if(s[i][j]==0)insert(num,j,1,0),d[j]++;            else insert(num,i,1,0),insert(num,j,1,0),d[i]++,d[j]++;        }    for(int i=1;i<=n;i++)        for(int j=0;j<d[i];j++)            insert(i,W,1,j);    while(spfa())mcf();    cout<<n*(n-1)*(n-2)/6-ans<<endl;    num=n;    for(int i=1;i<=n;i++)    {        for(int j=1;j<=n;j++)        {            if(i==j){printf("0 ");continue;}            if(i>j){printf("%d ",!s[j][i]);continue;}            num++; int x=find(num);            if(x==i)s[i][j]=1;else s[i][j]=0;            printf("%d ",s[i][j]);        }        puts("");    }           return 0;}