1. 程式人生 > >1061: [Noi2008]誌願者招募

1061: [Noi2008]誌願者招募

span col ins 每天 pan time memset stat noi

Time Limit: 20 Sec Memory Limit: 162 MB
Submit: 5742 Solved: 3449
[Submit][Status][Discuss]

Description

  申奧成功後,布布經過不懈努力,終於成為奧組委下屬公司人力資源部門的主管。布布剛上任就遇到了一個難 題:為即將啟動的奧運新項目招募一批短期誌願者。經過估算,這個項目需要N 天才能完成,其中第i 天至少需要 Ai 個人。 布布通過了解得知,一共有M 類誌願者可以招募。其中第i 類可以從第Si 天工作到第Ti 天,招募費用 是每人Ci 元。新官上任三把火,為了出色地完成自己的工作,布布希望用盡量少的費用招募足夠的誌願者,但這
並不是他的特長!於是布布找到了你,希望你幫他設計一種最優的招募方案。

Input

  第一行包含兩個整數N, M,表示完成項目的天數和可以招募的誌願者的種類。 接下來的一行中包含N 個非負 整數,表示每天至少需要的誌願者人數。 接下來的M 行中每行包含三個整數Si, Ti, Ci,含義如上文所述。為了 方便起見,我們可以認為每類誌願者的數量都是無限多的。

Output

  僅包含一個整數,表示你所設計的最優方案的總費用。

Sample Input

3 3
2 3 4
1 2 2
2 3 5
3 3 2

Sample Output

14

HINT

1 ≤ N ≤ 1000,1 ≤ M ≤ 10000,題目中其他所涉及的數據均 不超過2^31-1。

網絡流解法(因為有板子):

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<queue>
 5 using namespace std;
 6 
 7 const int INF=0x7f7f7f7f;
 8 const int T=1004;
 9 const int N=1005
; 10 const int M=50005; 11 12 struct Edge 13 { 14 int from,to,v,c,next;//from點到to點 v代表容量 c代表花費 next下一條邊 15 }E[M]; 16 int node=1; 17 int head[N],from[N],dis[N],vis[N]; 18 19 int n,m,ans; 20 21 void ins(int from,int to,int v,int c) 22 { 23 node++; 24 E[node]=(Edge){from,to,v,c,head[from]}; 25 head[from]=node; 26 } 27 28 void insert(int from,int to,int v,int c) 29 { 30 ins(from,to,v,c);ins(to,from,0,-c); 31 } 32 33 bool spfa() 34 { 35 queue<int> Q; 36 memset(dis,0x7f,sizeof(dis)); 37 Q.push(0);dis[0]=0;vis[0]=1; 38 while(!Q.empty()) 39 { 40 int q=Q.front();Q.pop(); 41 for(int i=head[q];i;i=E[i].next) 42 if(E[i].v>0&&dis[q]+E[i].c<dis[E[i].to]) 43 { 44 dis[E[i].to]=dis[q]+E[i].c; 45 from[E[i].to]=i; 46 if(!vis[E[i].to]) 47 { 48 Q.push(E[i].to); 49 vis[E[i].to]=1; 50 } 51 } 52 vis[q]=0; 53 } 54 return dis[T]!=INF; 55 } 56 57 void mcf() 58 { 59 int x=INF; 60 for(int i=from[T];i;i=from[E[i].from]) 61 x=min(E[i].v,x); 62 for(int i=from[T];i;i=from[E[i].from]) 63 { 64 ans+=x*E[i].c; 65 E[i].v-=x;E[i^1].v+=x; 66 } 67 } 68 69 int main() 70 { 71 scanf("%d %d",&n,&m); 72 int l=0,r; 73 for(int i=1;i<=n;i++) 74 { 75 scanf("%d",&r); 76 int x=r-l; 77 if(x>0) insert(0,i,x,0); 78 else insert(i,T,-x,0); 79 insert(i+1,i,INF,0); 80 l=r; 81 } 82 insert(n+1,T,l,0); 83 for(int i=1;i<=m;i++) 84 { 85 int s,t,c; 86 scanf("%d %d %d",&s,&t,&c); 87 insert(s,t+1,INF,c); 88 } 89 while(spfa()) mcf(); 90 printf("%d\n",ans); 91 return 0; 92 }

哪天回來補寫一個線性規劃 單純性算法 寫法

1061: [Noi2008]誌願者招募