1. 程式人生 > >bzoj1061 [Noi2008]志願者招募(網路流解決線性規劃問題)

bzoj1061 [Noi2008]志願者招募(網路流解決線性規劃問題)

bzoj1061 [Noi2008]志願者招募

題意:
一個專案需要N 天才能完成,其中第i 天至少需要Ai 個人。一共有M 類志願者可以招募。其中第i 類可以從第Si 天工作到第Ti 天,招募費用是每人Ci 元。布布希望用盡量少的費用招募足夠的志願者,求最優的招募方案的費用。

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

就是構造出說每個變數一正一負各出現一次的式子。
然後對於>=的處理新增輔助變數。
對於正負變數看成流入流出,滿足不等式 …..=0就是流量平衡。

程式碼:

#include<algorithm>
#include<iostream> #include<cstring> #include<cstdio> #include<queue> using namespace std; const int N=505; const int M=805; const int inf=0x3f3f3f3f; queue<int> Q; int n,m,head[N],to[4*M],nxt[4*M],w[4*M],id,U[M],V[M],W[M],S,T,num=1,dep[N]; bool vis[N]; void build(int u,int
v,int ww) { num++; to[num]=v; nxt[num]=head[u]; head[u]=num; w[num]=ww; num++; to[num]=u; nxt[num]=head[v]; head[v]=num; w[num]=0; } bool bfs() { memset(vis,0,sizeof(vis)); Q.push(S); vis[S]=1; dep[S]=1; while(!Q.empty()) { int u=Q.front(); Q.pop(); for
(int i=head[u];i;i=nxt[i]) { int v=to[i]; if(w[i]<=0||vis[v]) continue; vis[v]=1; dep[v]=dep[u]+1; Q.push(v); } } return vis[T]; } int dfs(int u,int d) { if(!d||u==T) return d; int ret=0; for(int i=head[u];i;i=nxt[i]) { int v=to[i]; if(w[i]<=0||dep[v]!=dep[u]+1) continue; int flow=dfs(v,min(d,w[i])); ret+=flow; d-=flow; w[i]-=flow; w[i^1]+=flow; if(!d) break; } if(!ret) dep[u]=-1; return ret; } int main() { scanf("%d%d%d",&n,&m,&id); for(int i=1;i<=m;i++) scanf("%d%d%d",&U[i],&V[i],&W[i]); for(int i=1;i<=m;i++) { if(i==id) S=U[i],T=V[i]; else if(W[i]<=W[id]) build(U[i],V[i],W[id]-W[i]+1),build(V[i],U[i],W[id]-W[i]+1); } int ret=0; while(bfs()) ret+=dfs(S,inf); printf("%d\n",ret); return 0; }