【BZOJ2597】【Wc2007】剪刀石头布 费用流，没写zkw卡时过
题解，比较常规式是费用流，没写那个神贪心。
首先是三元环需要取补集，先C(n,3)算出总环数，然后减去失败的三元环。
我们发现在一个三元环中，如果有某个选手入度（or 出度）为2，那么就会破坏这个环。
所以当选手有x个度时，就会破坏∑(i∈[0,x])个三元环（0+....+x-1），P.S.这个是C(x,2)推出来的。
建图：
第一层源点
到下一层每个点若干条边
条数： 那个选手可能赢的次数,
费用： 第i条边费用i-1，
流量： 1
第二层选手
到可能赢的每个对应竞赛建边
条数： 1
费用： 0
流量： 1
第三层比赛
到汇点连边
条数： 1
费用： 0
流量： 1
第四层汇点
完事。
除了1、2层中间的多建边有点思想成分，其它都很渣。
不多说。
至于那个方案数，就是check2、3层间边的满流和对应。
呃，来个图告诉一下大家zkw费用流和裸费用流的差距。
zkw费用流是wyfcyx写的。
然后popoqqq和eolv写的也是费用流，……连小常数君eolv都16s……醉了啊！！！
他当前博客http://wyfcyx.is-programmer.com/
代码：
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define N 5500
#define M 50000
#define P 105
#define inf 0x3f3f3f3f
using namespace std;
struct KSD
{
int u,v,len,fee,next;
}e[M];
int head[N],cnt;
void add(int u,int v,int w,int len)
{
cnt++;
e[cnt].u=u;
e[cnt].v=v;
e[cnt].len=len;
e[cnt].fee=w;
e[cnt].next=head[u];
head[u]=cnt;
}
int s,t,dist[N];
int pre[N],lim[N];
bool in[N];
queue<int>q;
int spfa()
{
memset(dist,0x3f,sizeof(dist));
while(!q.empty())q.pop();
int i,u,v;
dist[s]=0,in[s]=1;
lim[s]=inf;
q.push(s);
while(!q.empty())
{
u=q.front(),q.pop(),in[u]=0;
for(i=head[u];i;i=e[i].next)
{
v=e[i].v;
if(!e[i].len)continue;
if(dist[v]>dist[u]+e[i].fee)
{
dist[v]=dist[u]+e[i].fee;
lim[v]=min(e[i].len,lim[u]);
pre[v]=i;
if(!in[v])
{
in[v]=1;
q.push(v);
}
}
}
}
return dist[t];
}
void handle(int flow)
{
for(int i=pre[t];i;i=pre[e[i].u])
{
e[i].len-=flow;
e[i^1].len+=flow;
}
}
int map[P][P],ans[P][P],minfee;
int n,m,sum[P],id[P][P],idx[N],idy[N];
int win[P];
int main()
{
// freopen("hand.in","r",stdin);
// freopen("hand.out","w",stdout);
cnt=1;
int i,j,k,cnt,v;
scanf("%d",&n);
m=n*(n-1)/2,minfee=n*(n-1)*(n-2)/6;
s=n+m+1,t=n+m+2,cnt=n;
for(i=1;i<=n;i++)for(j=i+1;j<=n;j++)id[i][j]=++cnt,idx[cnt]=i,idy[cnt]=j,add(cnt,t,0,1),add(t,cnt,0,0);
for(i=1;i<=n;i++)for(j=1;j<=n;j++)scanf("%d",&map[i][j]);
for(i=1;i<n;i++)for(j=i+1;j<=n;j++)
{
if(map[i][j])add(i,id[i][j],0,1),add(id[i][j],i,0,0),win[i]++;
if(map[i][j]-1)add(j,id[i][j],0,1),add(id[i][j],j,0,0),win[j]++;
}
for(i=1;i<=n;i++)for(j=1;j<=win[i];j++)add(s,i,j-1,1),add(i,s,1-j,0);
while(spfa(),dist[t]<inf)
{
minfee-=dist[t]*lim[t];
handle(lim[t]);
}
printf("%d\n",minfee);
for(k=n+1;k<s;k++)
{
for(i=head[k];i;i=e[i].next)
{
v=e[i].v;
if(v!=t&&e[i].len)
{
if(v==idx[k])ans[idx[k]][idy[k]]=1,ans[idy[k]][idx[k]]=0;
else ans[idx[k]][idy[k]]=0,ans[idy[k]][idx[k]]=1;
}
}
}
for(i=1;i<=n;i++)
{
for(j=1;j<n;j++)printf("%d ",ans[i][j]);
printf("%d\n",ans[i][n]);
}
return 0;
}