BZOJ2597: [Wc2007]剪刀石头布

BZOJ2597: [Wc2007]剪刀石头布

费用流 之 特殊性质建图

题解：

感觉数据范围像网络流之类的，但是根本不会做。。。

当直观建模失效的时候，就弃疗吧，一般要么就是线性规划建图，要么就是有特殊的性质。

我们注意到图中的边是完全的，这时候应该想到总共有C(n,3)个三角形，正难则反，减去不符合条件的就是答案。

定义a败给b等价于a->b，当一个三角形中的一个点x有两条出边时，这个三角形就是不符合条件的。

换言之，对于一个点x，它任选两条边都构成一个三角形，其中不符合条件的有C(out[x],2)个。

$$ans=C(n,3)-\sum{C(out[x],2)}$$

拆开组合数得：

$$ans=n(n-1)(n-2)/6-out[x](out[x]-1)/2$$

out[x]每增加1，ans减少out[x] (原)

因此转化成建图就是S->x，cap=1，cost依次为0、1、2……n-2，代表出度为1的时候损耗0，出度为2的时候损耗0+1，出度为3的时候损耗0+1+2，以此类推。

补全剩下的边：

x->可以或必须成为x的出边的“边点”，cap=1，cost=0

“边点”->T，cap=1，cost=0

跑费用流即可。

输出方案：x->连向“边点”的边满载则为0，否则为1.

Code：

#include <iostream>
#include <cstring>
#include <cstdio>
#define D(x) cout<<#x<<" = "<<x<<"  "
#define E cout<<endl
using namespace std;
const int N = 20005;
const int INF = 0x3f3f3f3f;

int n,f[105][105],S,T,mark[105][105],q[N],front,tail;

void read(int &x){
    x=0; char c;
    while(!isdigit(c=getchar()));
    x=c-'0';
    while(isdigit(c=getchar()))x=x*10+c-'0';
}

struct Edge{
    int from,to,next,cap,flow,cost;
} e[200005];
int head[N], ec=1;
void add(int a,int b,int cap,int cost){
    ec++; e[ec].cap=cap; e[ec].flow=0; e[ec].cost=cost; 
    e[ec].from=a; e[ec].to=b; e[ec].next=head[a]; head[a]=ec;
}
void add2(int a,int b,int cap,int cost){
    add(a,b,cap,cost); add(b,a,0,-cost);
}

int d[N],pre[N]; bool vis[N]; 
bool spfa(){
    memset(d,0x3f,sizeof(d));
    front=tail=0; q[tail++]=S; d[S]=0;
    while(front!=tail){
        int u=q[front++]; vis[u]=false;
        for(int i=head[u];i;i=e[i].next){
            if(e[i].cap==e[i].flow)continue;
            int v=e[i].to;
            if(d[v]>d[u]+e[i].cost){
                d[v]=d[u]+e[i].cost; pre[v]=i;
                if(!vis[v]){ vis[v]=true; if(tail==N)tail=0; q[tail++]=v; }
            }
        }
    }
    return d[T]!=INF;
}

int mxf(){
    int ans=0, flow=0, a;
    while(spfa()){
        a=INF;
        for(int i=pre[T];i;i=pre[e[i].from]) a=min(a,e[i].cap-e[i].flow);
        for(int i=pre[T];i;i=pre[e[i].from]){
            e[i].flow+=a; e[i^1].flow-=a; ans+=e[i].cost*a;
        }
        flow+=a;
//      D(flow); D(ans); E;
    }
    return ans;
}

void build(){
    memset(head,0,sizeof(head)); ec=1;
    S=N-5; T=S+1;
    for(int i=1;i<=n;i++)
        for(int j=0;j<=n-2;j++)
            add2(S,i,1,j);
    int cnt=n;
    for(int i=1;i<=n;i++)
        for(int j=1;j<i;j++){
            cnt++;
            add2(cnt,T,1,0);
            if(f[i][j]==0 || f[i][j]==2) add2(i,cnt,1,0), mark[i][j]=ec-1;
            if(f[i][j]==1 || f[i][j]==2) add2(j,cnt,1,0), mark[j][i]=ec-1;
        }
}

int main(){
    freopen("a.in","r",stdin);
    freopen("a.out","w",stdout);
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            read(f[i][j]);
    build(); int ans=n*(n-1)*(n-2)/6 - mxf();
    printf("%d\n",ans);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=n;j++){
            if(i==j)printf("0%c",j==n?'\n':' ');
            else printf("%d%c",e[mark[i][j]].flow?0:1,j==n?'\n':' ');
        }
    }
}