首页 写随笔

cdcq(本博客废弃!现用博客:https://www.cnblogs.com/cdcq/)

本博客废弃!现用博客:https://www.cnblogs.com/cdcq/

导航

【BZOJ3144】【HNOI2013】切糕

总算做了一道2011以后的省选题了……
原题:

图片题面好评!

P,Q,R≤40,0≤D≤R,给出的所有的不和谐值不超过1000。

文本样例好评!

 

恩这个是听妹主席讲过后会写的,首先把每个点拆成链,那么割掉这个链上的某条边就表示这个点选了某个权值,边的流量就是这个点设成这个值的花费,表示要花费掉这些代价来把这条边割掉

相邻两个点的值之差<=d怎么搞呐,假设两个相邻的点x,y,那么x的第i个点往y的i-d个点连oo的边,即可,大概就像酱紫:

 

 如果有两个被割掉的线段距离超过d,就依然存在一条路径从s到t,就要继续割

差不多是酱紫:

恩就酱

代码:

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #include<cstring>
 5 #include<cmath>
 6 using namespace std;
 7 const int oo=168430090;
 8 int rd(){int z=0,mk=1;  char ch=getchar();
 9     while(ch<'0'||ch>'9'){if(ch=='-')mk=-1;  ch=getchar();}
10     while(ch>='0'&&ch<='9'){z=(z<<3)+(z<<1)+ch-'0';  ch=getchar();}
11     return z*mk;
12 }
13 struct ddd{int nxt,y,v,rvs;}e[2100000];  int lk[66100],ltp=0;
14 inline void ist(int x,int y,int z){
15     e[++ltp].nxt=lk[x],lk[x]=ltp,e[ltp].y=y,e[ltp].v=z,e[ltp].rvs=ltp+1;
16     e[++ltp].nxt=lk[y],lk[y]=ltp,e[ltp].y=x,e[ltp].v=0,e[ltp].rvs=ltp-1;
17 }
18 const int fx[4]={1,-1,0,0},fy[4]={0,0,1,-1};
19 int n,m,h,d;  int N;  int s,t;
20 int lvl[66100];
21 int q[66100],hd=0;
22 bool gtlvl(){
23     memset(lvl,0,sizeof(lvl));
24     q[hd=1]=s,lvl[s]=1;
25     for(int k=1;k<=hd;++k)
26         for(int i=lk[q[k]];i;i=e[i].nxt)if(e[i].v && !lvl[e[i].y])
27             lvl[e[i].y]=lvl[q[k]]+1,q[++hd]=e[i].y;
28     return lvl[t];
29 }
30 int mxflw(int x,int y){
31     if(x==t)  return y;
32     int bwl=0,flw=0;
33     for(int i=lk[x];i && bwl<y;i=e[i].nxt)if(e[i].v && lvl[e[i].y]==lvl[x]+1)
34         if((flw=mxflw(e[i].y,min(y-bwl,e[i].v)))){
35             bwl+=flw;
36             e[i].v-=flw,e[e[i].rvs].v+=flw;
37         }
38     if(!bwl)  lvl[x]=0;
39     return bwl;
40 }
41 int dnc(){
42     int bwl=0,flw=0;
43     while(gtlvl())while((flw=mxflw(s,oo)))  bwl+=flw;
44     return bwl;
45 }
46 bool chck(int x,int y,int z){
47     return x+fx[z]>=1 && x+fx[z]<=n && y+fy[z]>=1 && y+fy[z]<=m;
48 }
49 inline int gtid(int x,int y,int z){  return z*N+(x-1)*m+y;}
50 int main(){//freopen("ddd.in","r",stdin);
51     cin>>n>>m>>h>>d;  N=n*m;  s=0,t=N*(h+1)+1;
52     for(int k=1;k<=h;++k)for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)
53         ist(gtid(i,j,k-1),gtid(i,j,k),rd());
54     for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)for(int k=0;k<4;++k)if(chck(i,j,k))
55         for(int p=1;p+d<=h;++p){
56             ist(gtid(i,j,p+d),gtid(i+fx[k],j+fy[k],p),oo);
57             ist(gtid(i+fx[k],j+fy[k],p+d),gtid(i,j,p),oo);
58         }
59     for(int i=1;i<=n;++i)for(int j=1;j<=m;++j)  ist(s,gtid(i,j,0),oo),ist(gtid(i,j,h),t,oo);
60     cout<<dnc()<<endl;
61     return 0;
62 }
View Code

 

posted on 2017-03-13 15:27  cdcq_old  阅读(177)  评论(0编辑  收藏  举报