zoukankan      html  css  js  c++  java
  • 【题解】P3980 [NOI2008]志愿者招募(费用流求线性规划)

    【题解】P3980 [NOI2008]志愿者招募(费用流求线性规划)

    题意

    (m)种志愿者以及(n)天,每种志愿者有固定的服务区间(天),雇佣一个志愿者有不同的费用(c_i)。假设每种志愿者无限多,第(i)天至少要有(a[i])个志愿者,问最小花费

    加入辅助变量(y_i)把问题化为(=a[i])。设(i)种志愿者总共(x[i])个,(b[i][j])表示(i)号志愿者是否存在于(j)天,那么现在我们就是一个线性规划问题了,具体的是:

    [min X C^T\ mathrm{s.t.}\ (B,B_Y)(X,Y)^T= A^T ]

    写出其中一个出来

    [sum_{jin S_i} x_j-y_i= a_i ]

    那么

    [sum_{jin S_i} x_j- a_i-y_i=0 ]

    然而我们发现,(x_j)的存在是一段一段的((B)矩阵的行的1是连续连续的一段),所以我们对所有的等式进行差分。可以发现所有的变量不变量都以(+-)的形式分别出现一次(除了最后一个方程中的变量),这启示我们可以用网络流,因为网络流实际上就是做一个流量平衡的工作。

    但是原来的最后一个方程所有变量都只出现了一次,于是我们把(-( ext{equation }n))也加入方程组中这是可以用网络流解这个问题的关键。

    所以现在,所有的变量都出现了两次,且分别是以+-的形式出现,于是我们可以把每个方程看做一个点,每个变量看做一条边,现在我们是要最小化(sum x_i c_i),那么直接给所有的(x_i)代表的边给上(c_i)的费用跑费用流即可。

    但是网络流只能解除整数解,但是根据(IP)理论,这个线性规划问题存在一组整数解。

    //@winlere
    #include<iostream>
    #include<cstdio>
    #include<cstring>
    #include<algorithm>
    #include<assert.h>
    #include<queue>
    
    using namespace std;  typedef long long ll;
    inline int qr(){
    	int ret=0,f=0,c=getchar();
    	while(!isdigit(c))f|=c==45,c=getchar();
    	while(isdigit(c)) ret=ret*10+c-48,c=getchar();
    	return f?-ret:ret;
    }
    const int maxn=1e4+5;
    const int inf=1e9;
    int head[maxn],last[maxn],L[maxn],R[maxn],C[maxn],A[maxn],sum[maxn],n,m,S,T;
    struct E{int to,nx,w,c;}e[maxn<<5];
    void add(int fr,int to,int w,int c){
    	static int cnt=1;
    	e[++cnt]=(E){to,head[fr],w,c}; head[fr]=cnt;
    	e[++cnt]=(E){fr,head[to],0,-c}; head[to]=cnt;
    }
    
    ll fl[maxn],d[maxn];
    bool in[maxn];
    ll Mincost(){
    	static queue<int> q;
    	ll ret=0;
    	while(1){
    		memset(last,0,sizeof last);
    		memset(d,0x3f,sizeof d);
    		memset(fl,0,sizeof fl);
    		memset(in,0,sizeof in);
    		while(q.size()) q.pop();
    		q.push(S); d[S]=0; fl[S]=inf;
    		while(q.size()){
    			int now=q.front();
    			q.pop(); in[now]=0;
    			if(now==T) continue;
    			for(int t=head[now];t;t=e[t].nx)
    				if(fl[now]&&e[t].w&&d[e[t].to]>d[now]+e[t].c){
    					last[e[t].to]=t;
    					fl[e[t].to]=min(fl[now],0ll+e[t].w);
    					d[e[t].to]=d[now]+e[t].c;
    					if(!in[e[t].to]) q.push(e[t].to),in[e[t].to]=1;
    				}
    		}
    		if(!fl[T]) break;
    		ret+=d[T]*fl[T];
    		for(int t=T;t!=S;t=e[last[t]^1].to)
    			e[last[t]].w-=fl[T],e[last[t]^1].w+=fl[T];
    	}
    	return ret;
    }
    
    int main(){
    	n=qr(); m=qr();
    	for(int t=1;t<=n;++t) A[t]=qr();
    	S=n+2,T=n+3;
    	for(int t=1;t<=m;++t){
    		L[t]=qr(),R[t]=qr(),C[t]=qr(),++sum[L[t]],--sum[R[t]+1];
    		add(R[t]+1,L[t],inf,C[t]);
    	}
    	for(int t=1;t<=n;++t){
    		sum[t]+=sum[t-1],assert(sum[t]>0||A[t]==0);
    		if(A[t]>A[t-1]) add(t,T,A[t]-A[t-1],0);
    		if(A[t]<A[t-1]) add(S,t,A[t-1]-A[t],0);
    		add(t,t+1,inf,0);
    	} add(S,n+1,A[n],0);
    	ll ans=Mincost();
    	cerr<<"ans=";
    	printf("%lld
    ",ans);
    	return 0;
    }
    
    
    

  • 相关阅读:
    [JSOI2007][BZOJ1031] 字符加密Cipher|后缀数组
    leetcode Flatten Binary Tree to Linked List
    leetcode Pascal's Triangle
    leetcode Triangle
    leetcode Valid Palindrome
    leetcode Word Ladder
    leetcode Longest Consecutive Sequence
    leetcode Sum Root to Leaf Numbers
    leetcode Clone Graph
    leetcode Evaluate Reverse Polish Notation
  • 原文地址:https://www.cnblogs.com/winlere/p/12203438.html
Copyright © 2011-2022 走看看