标签:
一开始想DP做法。。发现不会QAQ
暴力费用流挺好想的。。就是拆点限制经过次数,然后每个点的出点往能到达的点的入点连边。。跑个最大费用最大流。
但显然边数能达到n^2。。显而易见的优化就是,如果存在路径点a->b->c,那么a就没必要连往c了。(虽然优化完极限情况下边数也差不多是n^2了。。。)
但是有可能出现借道的情况。就从入点往出点连一条流量1,费用0的边。(不会被借道两次)
有一个神坑。。。这题用spfa跑最长路的时候,队列数组要开到n^2。。。。。因为这个原因卡了一个晚上
1 #include<cstdio> 2 #include<iostream> 3 #include<cstring> 4 #include<algorithm> 5 #include<cmath> 6 #define ll long long 7 using namespace std; 8 const int maxn=4023,inf=1002333333; 9 struct poi{int x,y;}a[maxn]; 10 struct zs{int pre;short too,flow,dis;}e[4233333];int tot,last[maxn]; 11 short dl[4233333]; 12 int dis[maxn]; 13 bool u[maxn],ins[maxn]; 14 int i,j,k,n,m,s,t,ss; 15 int ans; 16 17 int ra;char rx; 18 inline int read(){ 19 rx=getchar(),ra=0; 20 while(rx<‘0‘||rx>‘9‘)rx=getchar(); 21 while(rx>=‘0‘&&rx<=‘9‘)ra*=10,ra+=rx-48,rx=getchar();return ra; 22 } 23 24 bool spfa(){ 25 memset(dis,255,(t+1)<<2),dis[ss]=0,dl[1]=ss,u[ss]=1; 26 int l=0,r=1,i,now; 27 while(l<r&&!dis[ss]) 28 for(i=last[now=dl[++l]],u[now]=0;i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]<dis[now]+e[i].dis){ 29 dis[e[i].too]=dis[now]+e[i].dis; 30 if(!u[e[i].too])u[e[i].too]=1,dl[++r]=e[i].too; 31 } 32 return dis[t]>-1; 33 } 34 inline int min(int a,int b){return a<b?a:b;} 35 int dfs(int x,int mx){ 36 if(x==t)return mx; 37 int used=0,i,w;ins[x]=1; 38 for(i=last[x];i;i=e[i].pre)if(e[i].flow&&dis[e[i].too]==dis[x]+e[i].dis&&!ins[e[i].too]) 39 if((w=dfs(e[i].too,min(mx-used,(int)e[i].flow)))){ 40 e[i].flow-=w,e[i^1].flow+=w,ans+=w*e[i].dis,used+=w; 41 if(used==mx){ins[x]=0;return mx;} 42 } 43 dis[x]=-1,ins[x]=0;return used; 44 } 45 46 inline void insert(int a,int b,int c,int d){ 47 e[++tot].too=b,e[tot].flow=c,e[tot].dis= d,e[tot].pre=last[a],last[a]=tot, 48 e[++tot].too=a,e[tot].flow=0,e[tot].dis=-d,e[tot].pre=last[b],last[b]=tot; 49 } 50 bool cmp(poi a,poi b){return a.x<b.x||(a.x==b.x&&a.y<b.y);} 51 int main(){ 52 n=read(); 53 for(i=1;i<=n;i++) 54 a[i].x=read(),a[i].y=read(); 55 sort(a+1,a+1+n,cmp); 56 ss=0,s=n<<1|1,t=s+1,tot=1; 57 insert(ss,s,2,0); 58 for(i=1;i<=n;i++){ 59 insert(i,i+n,1,0),insert(i,i+n,1,1),insert(s,i,1,0),insert(i+n,t,1,0); 60 int mn=inf; 61 for(j=i+1;j<=n;j++) 62 if(a[j].y<mn&&a[j].y>=a[i].y)mn=a[j].y,insert(i+n,j,2,0); 63 } 64 while(spfa())dfs(ss,inf); 65 printf("%d\n",ans); 66 }
[bzoj1930] [Shoi2003]pacman 吃豆豆
标签:
原文地址:http://www.cnblogs.com/czllgzmzl/p/5641741.html