BZOJ 1061 志愿者招募 单纯形 OR 费用流

题意:某个项目需要N天才能完成 其中第i天至少需要Ai个人,一共有M类志愿者.其中第i类可以从第Si天工作到第Ti天,招募费用是每人Ci元.

n<=1000,m<=1e4.问招募足够的志愿者需要的最小费用?

利用要满足的等式 -> 流量平衡关系建图

设第i类招了x[i]人
则对样例有不等式
x[1]>=a[1]
x[1]+x[2]>=a[2]
x[2]+x[3]>=a[3]
若第i天超出的人数为d[i],
则x[1]-d[1]=a[1]
x[1]+x[2]-d[2]=a[2]
x[2]+x[3]-d[3]=a[3]
每个不等式化成等式.令a[0]=a[n+1]=0 后对这n+2个等式相减.
x[1]-d[1]=a[1]
x[2]-d[2]+d[1]=a[2]-a[1]
x[3]-x[1]-d[3]+d[2]=a[3]-a[2]
-x[2]-x[3]+d[3]=-a[3]
得到的这n+1个等式中 x[i],d[i]正好在一个等式p中为正,另外一个等式q中为负数.
回想网络流中的流量平衡关系 令为正代表流入,为负代表流出.则q->p有容量为inf边(容量无限制).流量为未知的x[i].d[i].


因为第i个点有条边流量为a[i]-a[i-1](常数),
若该常数为正 则令源点到i容量为a[i]-a[i-1].
若该常数为负 则令i到汇点容量为-(a[i]-a[i-1]).

源点到i和i到汇点的边一定为满流边.则只要求出该图的任意一个最大流都满足上面的等式(流量平衡)
题目要求segma(x[i]*c[i])最小 则对x[i]边加上权值c[i] 求出最小费用最大流即可.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int inf=0x3f3f3f3f;
const int N=2e6+5;
int n,m,cnt,head[N],s,t;
struct node{
	int u,to,next,flow,cost;
}e[N]; 
void add(int u,int v,int flow,int cost)
{
	e[cnt].u=u,e[cnt].to=v,e[cnt].flow=flow,e[cnt].cost=cost;
	e[cnt].next=head[u];
	head[u]=cnt++;
}
void Add(int u,int v,int flow,int cost)
{
	add(u,v,flow,cost);
	add(v,u,0,-cost);
}
int a[N],pre[N],f[N],inq[N],c[N];
int SPFA(int s,int t)
{
	memset(pre,-1,sizeof(pre));
	queue<int> q;
	for(int i=0;i<=t+1;i++)
		c[i]=inf,inq[i]=0,f[i]=0;
	q.push(s),c[s]=0,inq[s]=1,f[s]=inf;
	while(!q.empty())
	{
		int u=q.front();
		q.pop(),inq[u]=0;
		for(int i=head[u];i!=-1;i=e[i].next)
		{
			int v=e[i].to,w=e[i].cost;
			
			if(e[i].flow>0&&c[v]>c[u]+w)
			{
				c[v]=c[u]+w;
				f[v]=min(f[u],e[i].flow);
				pre[v]=i;
				if(!inq[v])
					inq[v]=1,q.push(v);
			}
		}
	}
	return f[t];
}
int MCMF(int s,int t)
{
	int ans=0,i,tmp;
	while(tmp=SPFA(s,t))
	{
		for(int i=pre[t];i!=-1;i=pre[e[i].u])
		{
			ans+=tmp*e[i].cost;
			e[i].flow-=tmp;
			e[i^1].flow+=tmp;
		}
	}
	return ans;
}
int main()
{
	scanf("%d%d",&n,&m);
	memset(head,-1,sizeof(head));
	cnt=0;
	s=0,t=n+2;
	for(int i=1;i<=n;i++)	scanf("%d",&a[i]);
	int u,v,w;
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d%d",&u,&v,&w);
		Add(u,v+1,inf,w);
	}
	for(int i=1;i<=n+1;i++)
	{
		int tmp=a[i]-a[i-1];
		if(tmp>=0)
			Add(s,i,tmp,0);
		else
			Add(i,t,-tmp,0);
		if(i>1)
			Add(i,i-1,inf,0);
	}
	printf("%d\n",MCMF(s,t));
	return 0;
}

单纯形的对偶问题.

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const double eps=1e-7;
const int inf=0x3f3f3f3f,N=1005,M=1e4+5;
int n,m;
double a[M][N],b[M],c[N],v;
inline int read(){
    char c=getchar();int x=0,f=1;
    while(c<'0'||c>'9'){if(c=='-')f=-1; c=getchar();}
    while(c>='0'&&c<='9'){x=x*10+c-'0'; c=getchar();}
    return x*f;
}
void pivot(int l,int e){
    b[l]/=a[l][e];
    for(int j=1;j<=n;j++) if(j!=e) a[l][j]/=a[l][e];
    a[l][e]=1/a[l][e];
    
    for(int i=1;i<=m;i++) if(i!=l&&fabs(a[i][e])>0){
        b[i]-=a[i][e]*b[l];
        for(int j=1;j<=n;j++) if(j!=e) a[i][j]-=a[i][e]*a[l][j];
        a[i][e]=-a[i][e]*a[l][e];
    }
    
    v+=c[e]*b[l];
    for(int j=1;j<=n;j++) if(j!=e) c[j]-=c[e]*a[l][j];
    c[e]=-c[e]*a[l][e];
    
    //swap(B[l],N[e])
}
double simplex()
{
	while(true)
	{
		int e=0,l=0;
		for(e=1;e<=n;e++)	if(c[e]>eps) break;
		if(e==n+1) return v;
		double mn=inf;
		for(int i=1;i<=m;i++)
		{
			if(a[i][e]>eps&&mn>b[i]/a[i][e])
			{
				mn=b[i]/a[i][e];
				l=i;	
			}	
		} 
		if(mn==inf)	return inf;
		pivot(l,e);
	}
}
int main()
{
	n=read(),m=read();
	for(int i=1;i<=n;i++)
		c[i]=read();
	for(int i=1;i<=m;i++)
	{
		int s=read(),t=read();
		for(int j=s;j<=t;j++)
			a[i][j]=1;
		b[i]=read();
	}
	printf("%d\n",(int)(simplex()+0.5));
	return 0;
}