线性规划||网络流(费用流):COGS 288. [NOI2008] 志愿者招募

时间:2022-09-22 20:01:05

[NOI2008] 志愿者招募

输入文件:employee.in   输出文件:employee.out   简单对比
时间限制:2 s  
内存限制:512 MB

【问题描述】

申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管。布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募一批短期志愿者。经过估算,这个项目需要N 天才能完成,其中第i 天至少需要Ai 个人。
布布通过了解得知,一共有M 类志愿者可以招募。其中第i 类可以从第Si 天工作到第Ti 天,招募费用是每人Ci
元。新官上任三把火,为了出色地完成自己的工作,布布希望用尽量少的费用招募足够的志愿者,但这并不是他的特长!于是布布找到了你,希望你帮他设计一种最
优的招募方案。

【输入格式】

输入文件的第一行包含两个整数N, M,表示完成项目的天数和可以招募的志愿者的种类。
接下来的一行中包含N 个非负整数,表示每天至少需要的志愿者人数。
接下来的M 行中每行包含三个整数Si, Ti, Ci,含义如上文所述。为了方便起见,我们可以认为每类志愿者的数量都是无限多的。

【输出格式】
输入文件中仅包含一个整数,表示你所设计的最优方案的总费用。

【输入样例】
3 3
2 3 4
1 2 2
2 3 5
3 3 2
【输出样例】
14

【样例说明】
招募3 名第一类志愿者和4 名第三类志愿者。
【数据规模和约定】
30%的数据中,1 ≤ N, M ≤ 10,1 ≤ Ai ≤ 10;
100%的数据中,1 ≤ N ≤ 1000,1 ≤ M ≤ 10000,题目中其他所涉及的数据均
不超过2^31-1。

  https://www.byvoid.com/zhs/blog/noi-2008-employee

 #include <iostream>
#include <cstring>
#include <cstdio>
#include <queue>
using namespace std;
const int INF=;
const int maxn=;
const int maxm=;
int cnt=,fir[maxn],to[maxm],nxt[maxm],cap[maxm],val[maxm];
void addedge(int a,int b,int c,int v){
nxt[++cnt]=fir[a];to[cnt]=b;cap[cnt]=c;val[cnt]=v;fir[a]=cnt;
}
int N,M,S,T;
int dis[maxn],path[maxn],tot[maxn],vis[maxn]; queue<int>q;
int MCMF(){
int ret=;
while(true){
memset(dis,,sizeof(dis));dis[S]=;
q.push(S);vis[S]=;
while(!q.empty()){
int node=q.front();q.pop();vis[node]=;
for(int i=fir[node];i;i=nxt[i])
if(cap[i]&&dis[to[i]]>dis[node]+val[i]){
dis[to[i]]=dis[node]+val[i];
path[to[i]]=i;
if(!vis[to[i]]){
vis[to[i]]=;
q.push(to[i]);
}
}
}
if(dis[T]==)
break; int p=T,f=INF;
while(p!=S){
f=min(f,cap[path[p]]);
p=to[path[p]^];
}
ret+=dis[T]*f;p=T;
while(p!=S){
cap[path[p]]-=f;
cap[path[p]^]+=f;
p=to[path[p]^];
}
}
return ret;
} int main(){
freopen("employee.in","r",stdin);
freopen("employee.out","w",stdout);
scanf("%d%d",&N,&M);
S=;T=N+;
for(int i=;i<=N;i++)
scanf("%d",&tot[i]);
for(int i=,a,b,c;i<=M;i++){
scanf("%d%d%d",&a,&b,&c);
addedge(a,b+,INF,c);
addedge(b+,a,,-c);
}
for(int i=;i<=N+;i++){
int c=tot[i]-tot[i-];
if(c>){
addedge(S,i,c,);
addedge(i,S,,);
}
if(c<){
addedge(i,T,-c,);
addedge(T,i,,);
}
if(i>){
addedge(i,i-,INF,);
addedge(i-,i,,);
}
}
printf("%d\n",MCMF());
return ;
}

  然后又用线性规划。

  额,松弛型先天不足,数组开不了,无法通过,再学一下吧。

 #include <iostream>
#include <cstring>
#include <cstdio>
#include <cmath>
using namespace std;
const int maxr=;
const int maxc=; const double eps=1e-;
int n,m,nxt[maxc];
int N[maxr],B[maxr];
double a[maxr][maxr],v;
double b[maxr],c[maxr]; int SGN(double x){
return (x>eps)-(x<-eps);
} void Init(){
B[]=N[]=;v=0.0;
for(int i=;i<=n;i++)N[++N[]]=i;
for(int i=;i<=m;i++)B[++B[]]=i+n;
} void Pivot(int l,int e){
b[e]=b[l]/a[l][e];
a[e][l]=1.0/a[l][e];
for(int i=;i<=N[];i++)
if(N[i]!=e)a[e][N[i]]=a[l][N[i]]/a[l][e]; int pre=;
for(int i=;i<=N[];i++)
if(N[i]!=e&&SGN(a[e][N[i]])!=)
{nxt[pre]=i;pre=i;}
nxt[pre]=; for(int i=;i<=B[];i++)
if(B[i]!=l){
b[B[i]]-=a[B[i]][e]*b[e];
a[B[i]][l]=-a[B[i]][e]*a[e][l];
for(int j=nxt[];j;j=nxt[j])
if(N[j]!=e)a[B[i]][N[j]]-=a[B[i]][e]*a[e][N[j]];
} v+=c[e]*b[e];
c[l]=-c[e]*a[e][l];
for(int i=;i<=N[];i++)
if(N[i]!=e)
c[N[i]]-=c[e]*a[e][N[i]];
for(int i=;i<=N[];i++)if(N[i]==e)N[i]=l;
for(int i=;i<=B[];i++)if(B[i]==l)B[i]=e;
} void Simplex(){
while(true){
int e=maxr,l=maxr;
for(int i=;i<=N[];i++)
if(SGN(c[N[i]])>&&e>N[i])e=N[i]; if(e==maxr)break; double lam=-;
for(int i=;i<=B[];i++)
if(SGN(a[B[i]][e])>){
double tmp=b[B[i]]/a[B[i]][e];
if(lam==-||SGN(lam-tmp)>||SGN(lam-tmp)==&&l>B[i])
{lam=tmp;l=B[i];}
} Pivot(l,e);
}
} int main(){
#ifndef ONLINE_JUDGE
freopen("employee.in","r",stdin);
freopen("employee.out","w",stdout);
#endif
scanf("%d%d",&n,&m);
for(int i=,t;i<=n;i++){
scanf("%d",&t);
c[i]=t;
} for(int i=,l,r,t;i<=m;i++){
scanf("%d%d%d",&l,&r,&t);
for(int j=l;j<=r;j++)
a[i+n][j]=;
b[i+n]=t;
} Init();
Simplex(); printf("%d\n",(int)(v+0.5));
return ;
}

  这个程序可以AC,但解法并不具有共性啊……

 #include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int maxr=;
const int maxc=; int n,m,nxt[maxc];
int a[maxr][maxc]; void Pivot(int l,int e){
int pre=maxc-;
for(int i=;i<=n;i++)
if(a[l][i]!=){nxt[pre]=i;pre=i;}
nxt[pre]=-; for(int i=,t;i<=m;i++)
if(i!=l&&(t=a[i][e])){
a[i][e]=;
for(int j=nxt[maxc-];j!=-;j=nxt[j])
a[i][j]+=t*a[l][j];
}
} void Simplex(){
while(true){
int e=,l=;
for(int i=;i<=n;i++)
if(a[][i]>){e=i;break;}
if(e==)break;
for(int i=;i<=m;i++)
if(a[i][e]<&&(!l||a[l][]>a[i][]))
{l=i;} Pivot(l,e);
}
} int main(){
#ifndef ONLINE_JUDGE
freopen("employee.in","r",stdin);
freopen("employee.out","w",stdout);
#endif
scanf("%d%d",&n,&m);
for(int i=;i<=n;i++)
scanf("%d",&a[][i]);
for(int i=,l,r,t;i<=m;i++){
scanf("%d%d%d",&l,&r,&t);
for(int j=l;j<=r;j++)
a[i][j]=-;
a[i][]=t;
}
Simplex();
printf("%d\n",a[][]);
return ;
}

  这里的线性规划式子都是原式的对偶线性规划式。

线性规划||网络流(费用流):COGS 288. [NOI2008] 志愿者招募的更多相关文章

  1. 【费用流】BZOJ1061&colon; &lbrack;Noi2008&rsqb;志愿者招募(这题超好)

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 5291  Solved: 3173[Submit][Stat ...

  2. 线性规划费用流解法&lpar;Bzoj1061&colon; &lbrack;Noi2008&rsqb;志愿者招募&rpar;

    题面 传送门 Sol 线性规划费用流解法用与求解未知数为非负数的问题 这道题可以列出一堆形如 \(x[i]+x[j]+x[k]+...>=a[p]\) 的不等式 我们强行给每个式子减去一个东西, ...

  3. 【费用流】BZOJ1061&lbrack;NOI2008&rsqb;-志愿者招募

    [题目大意] 一个项目需要n天完成,其中第i天至少需要Ai个人.共有m类人可以招募,其中第i类可以从第Si天做到第Ti天,每人的招募费用为Ci元.求最小招募费用. [思路] byvoid神犇的建图详解 ...

  4. &lbrack;NOI2008&rsqb; 志愿者招募&lbrack;流量平衡&rsqb;

    288. [NOI2008] 志愿者招募 ★★★★   输入文件:employee.in   输出文件:employee.out   简单对比时间限制:2 s   内存限制:512 MB [问题描述] ...

  5. 【bzoj1061】&lbrack;NOI2008&rsqb;志愿者招募 线性规划与费用流

    题目描述 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难题:为即将启动的奥运新项目招募一批短期志愿者.经过估算,这个项目需要N 天才能完成,其中第i ...

  6. 网络流解线性规划问题 BZOJ1061&colon; &lbrack;Noi2008&rsqb;志愿者招募

    线性规划定义: 在给定有限的资源和竞争约束情况下,很多问题都可以表述为最大化或最小化某个目标.如果可以把目标指定为某些变量的线性函数,而且如果可以将资源约束指定为这些变量的等式或不等式,则得到了一个线 ...

  7. 从&lbrack;NOI2008志愿者招募&rsqb;浅谈线性规划在网络流构图上的巧用

    首先来看一下题..http://www.lydsy.com/JudgeOnline/problem.php?id=1061 1061: [Noi2008]志愿者招募 Description 申奥成功后 ...

  8. 【BZOJ 1061】 1061&colon; &lbrack;Noi2008&rsqb;志愿者招募 (线性规划与网络流)&ast;&ast;

    1061: [Noi2008]志愿者招募 Description 申奥成功后,布布经过不懈努力,终于成为奥组委下属公司人力资源部门的主管.布布刚上任就遇到了一个难 题:为即将启动的奥运新项目招募一批短 ...

  9. 【费用流】NOI2008志愿者招募

    1061: [Noi2008]志愿者招募 Time Limit: 20 Sec  Memory Limit: 162 MBSubmit: 5171  Solved: 3089[Submit][Stat ...

随机推荐

  1. 排序算法总结(四)快速排序【QUICK SORT】

    感觉自己这几篇都是主要参考的Wikipedia上的,快排就更加是了....wiki上的快排挺清晰并且容易理解的,需要注意的地方我也添加上了注释,大家可以直接看代码.需要注意的是,wikipedia上快 ...

  2. IOS图像拉伸解决方案

    UIButton实现背景拉伸,即图片两端不拉伸中间拉伸的办法有如下两种: 第一种方法很简单而且使用性更广.做法就是直接拉伸想要setBackgroundImage的image,代码如下: UIImag ...

  3. 杂谈SharpDx中的WIC组件——我们需要WIC的图片编码功能么?

    在前文 SharpDX之Direct2D教程II——加载位图文件和保存位图文件 中,发现在VB2010中不能很好的运用SharpDx中的WIC组件进行图片的编码工作.可能是我的设置问题,也可能是Sha ...

  4. poj1094

    Sorting It All Out Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 29539   Accepted: 10 ...

  5. VS辅助工具

    tfs签入策略修改工具Team Foundation Server 2015 Power Tools(tfpt) 添加签入禁止策略 \*/bin \*/obj \*/Release \*/Debug ...

  6. 00002、div的文字垂直居中与背景的渐变

    1.div可以放多行的文字,当显示文字较少时,文字可垂直居中 text-align: center; display: table-cell; vertical-align: middle; text ...

  7. java jackson 忽略不存在的属性字段 和 按照属性名转json

    @JsonAutoDetect(fieldVisibility = Visibility.ANY, getterVisibility = Visibility.NONE, isGetterVisibi ...

  8. 最小生成树——Prim算法和Kruskal算法

    洛谷P3366 最小生成树板子题 这篇博客介绍两个算法:Prim算法和Kruskal算法,两个算法各有优劣 一般来说当图比较稀疏的时候,Kruskal算法比较快 而当图很密集,Prim算法就大显身手了 ...

  9. 【JVM】深度分析Java的ClassLoader机制(源码级别)

    原文:深度分析Java的ClassLoader机制(源码级别) 为了更好的理解类的加载机制,我们来深入研究一下ClassLoader和他的loadClass()方法. 源码分析 public abst ...

  10. WebSocketTest 异步通讯,实时返回数据

    using System;using System.Net;using System.Net.Sockets;using System.Text;using System.Threading;usin ...