3280: 小R的烦恼
Time Limit: 10 Sec Memory Limit: 512 MBSubmit: 332 Solved: 169
[Submit][Status][Discuss]
Description
小R最近遇上了大麻烦，他的程序设计挂科了。于是他只好找程设老师求情。善良的程设老师答应不挂他，但是要求小R帮助他一起解决一个难题。
问题是这样的，程设老师最近要进行一项邪恶的实验来证明P=NP，这个实验一共持续n天，第i天需要a[i]个研究生来给他搬砖。研究生毕竟也是人，所以雇佣研究生是需要钱的，机智的程设老师已经联系好了m所大学，第j所大学共有l[j]个研究生，同时雇佣这所大学的一个研究生需要p[j]元钱。
本来程设老师满心欢喜的以为，这样捡最便宜的max{a[i]}个研究生雇来，就可以完成实验；结果没想到，由于他要求硕士生们每天工作25个小时不许吃饭睡觉上厕所喝水说话咳嗽打喷嚏呼吸空气，因此一天下来给他搬砖的所有研究生都会进入濒死状态。濒死状态的研究生，毫无疑问，就不能再进行工作了。但是机智的老师早早联系好了k家医院，第i家医院医治一个濒死的研究生需要d[i]天,并且需要q[i]元钱。
现在，程设老师想要知道，最少花多少钱，能够在这n天中满足每天的需要呢？若无法满足，则请输出”impossible”。注意，由于程设老师良心大大的坏，所以他是可以不把濒死的研究生送去医院的！
Input
本题包含多组数据；第一行是一个数T(T<=11)，表示数据组数，以下T组数据。
对于每一组数据，第一行三个数，n,m,k;
以下一行n个数，表示a[1]…a[n]
接着一行2m个数，表示l[1],p[1]…l[n],p[n]
接着一行2k个数，表示d[1],q[1]…d[n],q[n]
Output
对于每组数据以样例的格式输出一行，两个数分别表示第几组数据和最少钱数。
Sample Input
2
3 2 1
10 20 30
40 90 15 100
1 5
3 2 1
10 20 30
40 90 15 100
2 5
3 2 1
10 20 30
40 90 15 100
1 5
3 2 1
10 20 30
40 90 15 100
2 5
Sample Output
Case 1: 4650
Case 2: impossible
Case 2: impossible
HINT
样例解释：买下90块钱的那40个研究生，另外再买10个100块钱的。这样，第一天用完的10个人全部送到医院，那么他们在第三天可以继续使用；同时，第二天和第三天都用新的研究生来弥补，这样一共需要花费40*90 + 10*100 + 5*10 = 4650元。
数据规模：
对于30%的数据中的每组数据,
满足n<=5,m,k<=2,其余数均小于等于100或者
n<=10,m,k<=10,其余数均小于等于20.
对于100%的数据
n,m,k<=50,其余数均小于等于100.
似乎这种题有大概2~3道一样的23333
1 #include<bits/stdc++.h> 2 #define rep(i,l,r) for(int i=l;i<=r;++i) 3 using namespace std; 4 const int N=102333,inf=2147483647; 5 struct Edge{ 6 int to,next,from,c,w; 7 }e[1000000]; 8 int head[N],tot=1,ans,dis[N],from[N],T,x,y,n,m,K; 9 bool used[N]; 10 inline void ins(int u,int v,int w,int cost) { 11 e[++tot].to=v; e[tot].next=head[u]; head[u]=tot; e[tot].w=w; e[tot].c=cost; e[tot].from=u; 12 } 13 14 inline bool spfa() { 15 queue<int> q; for(int i=0;i<=T;i++) dis[i]=inf; dis[0]=0; q.push(0); used[0]=1; 16 while(!q.empty()) { 17 int x=q.front(); q.pop(); 18 for(int k=head[x];k;k=e[k].next) 19 if(e[k].w>0&&dis[x]+e[k].c<dis[e[k].to]){ 20 dis[e[k].to]=dis[x]+e[k].c; from[e[k].to]=k; 21 if(!used[e[k].to]) { 22 used[e[k].to]=1; q.push(e[k].to); 23 } 24 } 25 used[x]=0; 26 } 27 if(dis[T]==inf) return 0;else return 1; 28 } 29 30 inline void run() { 31 int x=inf; 32 for(int k=from[T];k;k=from[e[k].from]) x=min(x,e[k].w); 33 for(int k=from[T];k;k=from[e[k].from]) { 34 e[k].w-=x; e[k^1].w+=x; ans+=e[k].c*x; 35 } 36 } 37 38 inline void insert(int u,int v,int w,int c){ 39 ins(u,v,w,c); ins(v,u,0,-c); 40 } 41 int main(){ 42 int t; 43 scanf("%d",&t); 44 rep(l,1,t){ 45 memset(head,0,sizeof head); tot=1; ans=0; 46 scanf("%d%d%d",&n,&m,&K); 47 T=1+n+n; 48 printf("Case %d: ",l); 49 rep(i,1,n) { 50 scanf("%d",&x); 51 insert(0,i+n,x,0); 52 insert(i,T,x,0); 53 if(i<n) insert(i,i+1,inf,0); 54 } 55 rep(i,1,m) scanf("%d%d",&x,&y),insert(0,1,x,y); 56 rep(i,1,K){ 57 scanf("%d%d",&x,&y); 58 rep(j,1,n) if(j+x+1<=n) insert(j+n,j+x+1,inf,y); 59 } 60 while(spfa()) run(); 61 bool flag=1; 62 for(int i=2;i<=tot;i+=2) if(e[i].to==T && e[i].w) {puts("impossible"),flag=0; break;} 63 if(flag) printf("%d ",ans); 64 } 65 }