bzoj4823: [Cqoi2017]老C的方块(最小割)

Posted CHerish_OI

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了bzoj4823: [Cqoi2017]老C的方块(最小割)相关的知识,希望对你有一定的参考价值。

4823: [Cqoi2017]老C的方块

题目:传送门 

 


 

 

题解:

   毒瘤题ORZ....

   太菜了看出来是最小割啥边都不会建...狂%大佬强强强

   

    黑白染色?不!是四个色一起染,四层图跑最小割。。。

   很惊奇的发现染完色之后只要是不喜欢的图形都一定可以由黄-->黑-->红-->绿 组成

   那就很nice啦...兴高采烈的去敲代码...结果10^5*10^5???搞毛线...太弱了ORZ,又看了一波大佬的操作,用map存!

   woc...不谈了不谈了...撸撸撸(分情况分到想屎...虽然不多)

   注意一下...黑到红连边的时候是不可以连inf的...然而我傻逼逼的打了...改成min(w[黑],w[红])...

 


 

 

代码:

  1 #include<cstdio>
  2 #include<cstring>
  3 #include<cstdlib>
  4 #include<cmath>
  5 #include<algorithm>
  6 #include<map>
  7 #define inf 999999999
  8 #define qread(x) x=read()
  9 using namespace std;
 10 typedef long long LL;
 11 inline LL read()
 12 {
 13     LL x=0,f=1;char ch;
 14     while(ch<\'0\' || ch>\'9\'){if(ch==\'-\')f=-1;ch=getchar();}
 15     while(ch>=\'0\' && ch<=\'9\'){x=x*10+ch-\'0\';ch=getchar();}
 16     return f*x;
 17 }
 18 struct node
 19 {
 20     int x,y,c,next,other;
 21 }a[1110000];int len,last[410000];
 22 int n,m,T,st,ed;
 23 void ins(int x,int y,int c)
 24 {
 25     int k1,k2;
 26     k1=++len;
 27     a[len].x=x;a[len].y=y;a[len].c=c;
 28     a[len].next=last[x];last[x]=len;
 29      
 30     k2=++len;
 31     a[len].x=y;a[len].y=x;a[len].c=0;
 32     a[len].next=last[y];last[y]=len;
 33      
 34     a[k1].other=k2;
 35     a[k2].other=k1;
 36 }
 37 int head,tail;
 38 int list[410000],h[410000];
 39 bool bt_h()
 40 {
 41     memset(h,0,sizeof(h));h[st]=1;
 42     list[1]=st;head=1;tail=2;
 43     while(head!=tail)
 44     {
 45         int x=list[head];
 46         for(int k=last[x];k;k=a[k].next)
 47         {
 48             int y=a[k].y;
 49             if(!h[y] && a[k].c)
 50             {
 51                 h[y]=h[x]+1;
 52                 list[tail++]=y;
 53             }
 54         }
 55         head++;
 56     }
 57     if(h[ed])return true;
 58     return false;
 59 }
 60 int find_flow(int x,int flow)
 61 {
 62     int s=0,t;
 63     if(x==ed)return flow;
 64     for(int k=last[x];k;k=a[k].next)
 65     {
 66         int y=a[k].y;
 67         if(h[y]==h[x]+1 && a[k].c && s<flow)
 68         {
 69             s+=t=find_flow(y,min(a[k].c,flow-s));
 70             a[k].c-=t;a[a[k].other].c+=t;
 71         }
 72     }
 73     if(!s)h[x]=0;
 74     return s;
 75 }
 76 int x[110000],y[110000],w[110000];//第x列,第y行,金币数 
 77 map<pair<int,int>,int>id;
 78 int main()
 79 {
 80     len=0;memset(last,0,sizeof(last));
 81     memset(w,0,sizeof(w));
 82     qread(n);qread(m);qread(T);st=T+1;ed=st+1;
 83     for(int i=1;i<=T;i++){qread(x[i]);qread(y[i]);qread(w[i]);id[make_pair(x[i],y[i])]=i;}
 84     //黄-->黑-->红-->绿 
 85     for(int i=1;i<=T;i++)
 86     {
 87         if(x[i]%4==0)//黑或黄 
 88         {
 89             if(y[i]%2==1)ins(st,i,w[i]);//
 90             else//
 91             {
 92                 if(id[make_pair(x[i]+1,y[i])])ins(id[make_pair(x[i]+1,y[i])],i,inf);
 93                 if(id[make_pair(x[i],y[i]+1)])ins(id[make_pair(x[i],y[i]+1)],i,inf);
 94                 if(id[make_pair(x[i],y[i]-1)])ins(id[make_pair(x[i],y[i]-1)],i,inf); 
 95             }
 96         }
 97         if(x[i]%4==1)//黑或黄 
 98         {
 99             if(y[i]%2==0)ins(st,i,w[i]);//
100             else//
101             {
102                 //黑-->黄 
103                 if(id[make_pair(x[i],y[i]-1)])ins(id[make_pair(x[i],y[i]-1)],i,inf);
104                 if(id[make_pair(x[i],y[i]+1)])ins(id[make_pair(x[i],y[i]+1)],i,inf);
105                 if(id[make_pair(x[i]-1,y[i])])ins(id[make_pair(x[i]-1,y[i])],i,inf);
106                 
107                 //黑-->红(计算w) 
108                 if(id[make_pair(x[i]+1,y[i])])
109                     ins(i,id[make_pair(x[i]+1,y[i])],min(w[i],w[id[make_pair(x[i]+1,y[i])]]));
110             }
111         }
112         if(x[i]%4==2)//红或绿
113         {
114             if(y[i]%2==1)//
115             {
116                 if(id[make_pair(x[i],y[i]-1)])ins(i,id[make_pair(x[i],y[i]-1)],inf);
117                 if(id[make_pair(x[i],y[i]+1)])ins(i,id[make_pair(x[i],y[i]+1)],inf);
118                 if(id[make_pair(x[i]+1,y[i])])ins(i,id[make_pair(x[i]+1,y[i])],inf);
119             }
120             else ins(i,ed,w[i]);//绿 
121         }
122         if(x[i]%4==3)//红或绿
123         {
124             if(y[i]%2==0)//
125             {
126                 //红-->绿 
127                 if(id[make_pair(x[i],y[i]-1)])ins(i,id[make_pair(x[i],y[i]-1)],inf);
128                 if(id[make_pair(x[i],y[i]+1)])ins(i,id[make_pair(x[i],y[i]+1)],inf);
129                 if(id[make_pair(x[i]-1,y[i])])ins(i,id[make_pair(x[i]-1,y[i])],inf);
130                 
131                 //黑-->红
132                 if(id[make_pair(x[i]+1,y[i])])
133                     ins(id[make_pair(x[i]+1,y[i])],i,min(w[i],w[id[make_pair(x[i]+1,y[i])]])); 
134             }
135             else ins(i,ed,w[i]);//绿 
136         }
137     }
138     int ans=0;
139     while(bt_h())ans+=find_flow(st,inf);
140     printf("%d\\n",ans);
141     return 0;
142 }

 

以上是关于bzoj4823: [Cqoi2017]老C的方块(最小割)的主要内容,如果未能解决你的问题,请参考以下文章

Bzoj4823 [Cqoi2017]老C的方块

BZOJ 4823: [Cqoi2017]老C的方块

[bzoj4823][洛谷P3756][Cqoi2017]老C的方块

bzoj 4823: [Cqoi2017]老C的方块最大权闭合子图

省选之前的未完成的计划(截至到省选)

@loj - 3022@ 「CQOI2017」老 C 的方块