申奥成功后，布布经过不懈努力，终于成为奥组委下属公司人力资源部门的主管。布布刚上任就遇到了一个难题：为即将启动的奥运新项目招募一批短期志愿者。经过估算，这个项目需要N 天才能完成，其中第i 天至少需要Ai 个人。 布布通过了解得知，一共有M 类志愿者可以招募。其中第i 类可以从第Si 天工作到第Ti 天，招募费用是每人Ci 元。新官上任三把火，为了出色地完成自己的工作，布布希望用尽量少的费用招募足够的志愿者，但这并不是他的特长！于是布布找到了你，希望你帮他设计一种最优的招募方案。
标签：
申奥成功后，布布经过不懈努力，终于成为奥组委下属公司人力资源部门的主管。布布刚上任就遇到了一个难题：为即将启动的奥运新项目招募一批短期志愿者。经过估算，这个项目需要N 天才能完成，其中第i 天至少需要Ai 个人。 布布通过了解得知，一共有M 类志愿者可以招募。其中第i 类可以从第Si 天工作到第Ti 天，招募费用是每人Ci 元。新官上任三把火，为了出色地完成自己的工作，布布希望用尽量少的费用招募足够的志愿者，但这并不是他的特长！于是布布找到了你，希望你帮他设计一种最优的招募方案。
第一行包含两个整数N, M，表示完成项目的天数和可以招募的志愿者的种类。 接下来的一行中包含N 个非负整数，表示每天至少需要的志愿者人数。 接下来的M 行中每行包含三个整数Si, Ti, Ci，含义如上文所述。为了方便起见，我们可以认为每类志愿者的数量都是无限多的。
仅包含一个整数，表示你所设计的最优方案的总费用。
招募第一类志愿者3名，第三类志愿者4名 30%的数据中，1 ≤ N, M ≤ 10，1 ≤ Ai ≤ 10； 100%的数据中，1 ≤ N ≤ 1000，1 ≤ M ≤ 10000，题目中其他所涉及的数据均 不超过2^31-1。
题解：网络流毕业题的难度啊。。。让我跟所有题解一样指向这里啦，膜BYVOID神犇。。。
具体代码实现的时候动一下脑筋，观察一下就好咯。
1 #include<iostream> 2 #include<cstdio> 3 #include<cmath> 4 #include<algorithm> 5 #include<stack> 6 #include<queue> 7 #include<cstring> 8 #define PAU putchar(‘ ‘) 9 #define ENT putchar(‘\n‘) 10 #define MSE(a,b) memset(a,b,sizeof(a)) 11 #define REN(x) for(ted*e=fch[x];e;e=e->nxt) 12 #define TIL(x) for(int i=1;i<=x;i++) 13 using namespace std; 14 const int maxn=1000+10,maxm=400000+10,inf=1e9; 15 struct zkw{ 16 struct ted{int x,y,w,c;ted*nxt,*re;}adj[maxm],*fch[maxn],*ms; 17 int n,S,T,d[maxn],cost,ans;bool inq[maxn],vis[maxn]; 18 void init(int n){this->n=n;ms=adj;MSE(vis,false);MSE(inq,false);return;} 19 void add(int x,int y,int w,int c){ 20 *ms=(ted){x,y,w,c,fch[x],ms+1};fch[x]=ms++; 21 *ms=(ted){y,x,0,-c,fch[y],ms-1};fch[y]=ms++; 22 return; 23 } 24 bool bfs(){ 25 for(int i=1;i<=n;i++)d[i]=inf;queue<int>Q;Q.push(T);d[T]=0; 26 while(!Q.empty()){ 27 int x=Q.front();Q.pop();inq[x]=false;REN(x){ 28 int v=e->y;if(e->re->w&&d[v]>d[x]+e->re->c){ 29 d[v]=d[x]+e->re->c;if(!inq[v])inq[v]=true,Q.push(v); 30 } 31 } 32 }for(ted*e=adj;e!=ms;e++)e->c+=d[e->y]-d[e->x];cost+=d[S];return d[S]!=inf; 33 } 34 int dfs(int x,int delta){ 35 if(x==T||!delta)return(ans+=delta*cost,delta);int flow=0,k;vis[x]=true;REN(x){ 36 int v=e->y;if(e->w&&!e->c&&!vis[v]&&(k=dfs(v,min(delta,e->w)))){ 37 e->w-=k;e->re->w+=k;flow+=k;delta-=k;if(!delta)break; 38 } 39 }return flow; 40 } 41 int mcmf(int S,int T){ 42 this->S=S;this->T=T;while(bfs())do MSE(vis,false);while(dfs(S,inf));return ans; 43 } 44 }sol; 45 inline int read(){ 46 int x=0;bool sig=true;char ch=getchar(); 47 for(;!isdigit(ch);ch=getchar())if(ch==‘-‘)sig=false; 48 for(;isdigit(ch);ch=getchar())x=10*x+ch-‘0‘;return sig?x:-x; 49 } 50 inline void write(int x){ 51 if(x==0){putchar(‘0‘);return;}if(x<0)putchar(‘-‘),x=-x; 52 int len=0;static int buf[20];while(x)buf[len++]=x%10,x/=10; 53 for(int i=len-1;i>=0;i--)putchar(buf[i]+‘0‘);return; 54 } 55 int n,m,need[maxn],S,T; 56 int main(){ 57 n=read();m=read();S=n+2;T=n+3;sol.init(n+3); 58 TIL(n)need[i]=read();int x,y,w; 59 TIL(m)x=read(),y=read(),w=read(),sol.add(x,y+1,inf,w),sol.add(y+1,x,0,-w); 60 n++;TIL(n){ 61 int t=need[i]-need[i-1]; 62 if(t>=0)sol.add(S,i,t,0),sol.add(i,S,0,0); 63 else sol.add(i,T,-t,0),sol.add(T,i,0,0); 64 if(i>1)sol.add(i,i-1,inf,0),sol.add(i-1,i,0,0); 65 }write(sol.mcmf(S,T)); 66 return 0; 67 }
标签：
原文地址：http://www.cnblogs.com/chxer/p/4740883.html
