标签：style blog http color io os ar for sp
某天，外星人展开了对地球的侵略，OIer们开始与之周旋。。。
外星人派出了172849个外星人乘着UFO来到地球，curimit神new带领着OIer们奋力抵抗。curimit神new觉得，仅仅抵抗外星人是不行的，因为外星人还有172849架UFO，必须从根源阻止外星人！
于是curimit神new交给小z一个任务：潜入外星人的基地，摧毁外星人的燃料库。
小z拿着curimit神new给他的地图，来到了外星人的基地。
外星人的燃料库由一排N个燃料筒组成，每个燃料筒中装有一些燃料。小z需要用炸弹炸毁燃料筒。有的燃料筒里面有许多燃料，炸毁它价值就比较高。但有的燃料筒里面只有很少的燃料，甚至没有燃料，所以炸毁它价值很低，甚至不值得去炸毁！于是curimit神new对每个燃料筒有一个估价，有正有负。但是外星人也不是没有防范。在燃料库中有172849个外星人在不停的巡逻，如果小z把任意连续M个燃料筒炸毁K个以上，外星人就会察觉，小z就遭殃了！
于是小z希望你能够帮他定出一个方案，在不被外星人察觉的情况下，能够使自己炸毁的燃料筒的估价和最大（如果所有筒的估价都是负的，那么一个不炸的估价和最大，是0）。
第一行有三个数，N，M，K。N<=1000，k<=m<=100。描述见题。
第二行有N个整数，绝对值小于1000，第i个数表示curimit神new对第i个燃料筒的估价。
一个数，为满足条件的最大估价和。
1 #include<cstdio> 2 #include<cstdlib> 3 #include<cmath> 4 #include<cstring> 5 #include<algorithm> 6 #include<iostream> 7 #include<vector> 8 #include<map> 9 #include<set> 10 #include<queue> 11 #include<string> 12 #define inf 1000000000 13 #define maxn 100000 14 #define maxm 100000 15 #define eps 1e-10 16 #define ll long long 17 #define pa pair<int,int> 18 #define for0(i,n) for(int i=0;i<=(n);i++) 19 #define for1(i,n) for(int i=1;i<=(n);i++) 20 #define for2(i,x,y) for(int i=(x);i<=(y);i++) 21 #define for3(i,x,y) for(int i=(x);i>=(y);i--) 22 #define mod 1000000007 23 using namespace std; 24 inline int read() 25 { 26 int x=0,f=1;char ch=getchar(); 27 while(ch<‘0‘||ch>‘9‘){if(ch==‘-‘)f=-1;ch=getchar();} 28 while(ch>=‘0‘&&ch<=‘9‘){x=10*x+ch-‘0‘;ch=getchar();} 29 return x*f; 30 } 31 int n,m,k,mincost,tot=1,s,t,a[maxn],head[maxn],q[maxm],d[maxn],from[2*maxm]; 32 bool v[maxn]; 33 struct edge{int from,next,go,v,c;}e[2*maxm]; 34 void ins(int x,int y,int z,int w) 35 { 36 e[++tot].go=y;e[tot].from=x;e[tot].v=z;e[tot].c=w;e[tot].next=head[x];head[x]=tot; 37 } 38 void insert(int x,int y,int z,int w) 39 { 40 ins(x,y,z,w);ins(y,x,0,-w); 41 } 42 bool spfa() 43 { 44 for (int i=s;i<=t;i++){v[i]=0;d[i]=inf;} 45 int l=0,r=1,y;q[1]=s;d[s]=0;v[0]=1; 46 while(l!=r) 47 { 48 int x=q[++l];if(l==maxn)l=0;v[x]=0; 49 for (int i=head[x];i;i=e[i].next) 50 if(e[i].v&&d[x]+e[i].c<d[y=e[i].go]) 51 { 52 d[y]=d[x]+e[i].c;from[y]=i; 53 if(!v[y]){v[y]=1;q[++r]=y;if(r==maxn)r=0;} 54 } 55 } 56 return d[t]!=inf; 57 } 58 void mcf() 59 { 60 while(spfa()) 61 { 62 int tmp=inf; 63 for(int i=from[t];i;i=from[e[i].from]) tmp=min(tmp,e[i].v); 64 mincost+=d[t]*tmp; 65 for(int i=from[t];i;i=from[e[i].from]){e[i].v-=tmp;e[i^1].v+=tmp;} 66 } 67 } 68 int main() 69 { 70 freopen("input.txt","r",stdin); 71 freopen("output.txt","w",stdout); 72 n=read();m=read();k=read(); 73 s=0;t=n+1; 74 for1(i,n)a[i]=read(); 75 insert(s,1,k,0);insert(n-m+2,t,k,0); 76 for1(i,m)insert(1,i+1,1,-a[i]); 77 for2(i,m+1,n-m)insert(i-m+1,i+1,1,-a[i]); 78 for2(i,n-m+1,n)insert(i-m+1,n-m+2,1,-a[i]); 79 for1(i,n-m+1)insert(i,i+1,inf,0); 80 mcf(); 81 printf("%d\n",-mincost); 82 return 0; 83 }
标签：style blog http color io os ar for sp
原文地址：http://www.cnblogs.com/zyfzyf/p/4035093.html
