BZOJ 2597 [Wc2007]剪刀石头布

题解：

http://wenku.baidu.com/view/0ad00abec77da26925c5b01c.html

 

吐槽这题数据，本地测全wa，就一直对拍，最后毛了，交了一发，ac了。。。

难道这个题是spj？

 

View Code

  1 #include <iostream>

  2 #include <cstring>

  3 #include <cstdio>

  4 #include <cstdlib>

  5 #include <algorithm>

  6 #include <ctime>

  7  

  8 #define N 20000

  9 #define M 2001000

 10  

 11 using namespace std;

 12  

 13 int head[N],to[M],next[M],len[M],pr[M];

 14 int pre[N],n,m,p,S,T,dis[N],q[M<<1],cnt;

 15 int kr[200][200],cha[N],fx[N],rd[N],num;

 16 int ts[200][200];

 17 bool vis[N];

 18 double ttt;

 19 

 20 inline void add(int u,int v,int r,int w)

 21 {

 22     to[cnt]=v; len[cnt]=r; pr[cnt]=w; next[cnt]=head[u]; head[u]=cnt++;

 23     to[cnt]=u; len[cnt]=0; pr[cnt]=-w; next[cnt]=head[v]; head[v]=cnt++;

 24 }

 25 

 26 inline void read()

 27 {

 28     scanf("%d",&n);

 29     memset(rd,0,sizeof rd);

 30     memset(head,-1,sizeof head); cnt=0;

 31     m=n*(n-1)/2;

 32     S=0; T=m+n+1; num=0;

 33     for(int i=1;i<=n;i++)

 34         for(int j=1;j<=n;j++)

 35         {

 36             scanf("%d",&kr[i][j]);

 37             if(j<=i) continue;

 38             num++;

 39             if(kr[i][j]==2) add(num,i+m,1,0),add(num,j+m,1,0),rd[i]++,rd[j]++;

 40             else if(kr[i][j]==1) add(num,i+m,1,0),rd[i]++;

 41             else add(num,j+m,1,0),rd[j]++;

 42         }

 43     for(int i=1;i<=m;i++) add(S,i,1,0);

 44     for(int i=1;i<=n;i++)

 45         for(int j=1;j<=rd[i];j++)

 46             add(i+m,T,1,cha[j]);

 47 }

 48  

 49 inline bool spfa()

 50 {

 51     memset(dis,0x3f,sizeof dis);

 52     memset(pre,-1,sizeof pre);

 53     int h=1,t=2,sta;

 54     q[1]=S; vis[S]=true; dis[S]=0;

 55     while(h<t)

 56     {

 57         sta=q[h++]; vis[sta]=false;

 58         for(int i=head[sta];~i;i=next[i])

 59             if(len[i]&&dis[to[i]]>dis[sta]+pr[i])

 60             {

 61                 dis[to[i]]=dis[sta]+pr[i];

 62                 pre[to[i]]=i;

 63                 if(!vis[to[i]]) vis[to[i]]=true,q[t++]=to[i];

 64             }

 65     }

 66     return pre[T]!=-1;

 67 }

 68  

 69 inline void updata()

 70 {

 71     for(int i=pre[T];~i;i=pre[to[i^1]])

 72     {

 73         len[i]-=1; len[i^1]+=1;

 74     }

 75 }

 76  

 77 inline int findwho(int u)

 78 {

 79     for(int i=head[u];~i;i=next[i])

 80         if(!len[i]) return to[i];

 81 }

 82  

 83 inline void go()

 84 {

 85     int ans=0;

 86     while(spfa()) ans+=dis[T],updata();

 87     printf("%d\n",(n*(n-1)*(n-2)/3+m-ans)/2);

 88     

 89     int now=0;

 90     for(int i=1,t;i<=n;i++)

 91         for(int j=i+1;j<=n;j++)

 92         {

 93             now++; t=findwho(now);

 94             ts[t-m][i+j-t+m]=1;

 95         }

 96 

 97     for(int i=1;i<=n;i++)

 98     {

 99         for(int j=1;j<=n;j++)

100             printf("%d ",ts[i][j]);

101         puts("");

102     }

103 } 

104  

105 inline void prep()

106 {

107     for(int i=1;i<=110;i++) fx[i]=i*i,cha[i]=fx[i]-fx[i-1];

108 }

109  

110 int main()

111 {

112     prep();

113     read();

114     go();

115     return 0;

116 }