ACM-ICPC 2018 焦作赛区网络预赛 - F - 最大费用流

题目链接:点击这里

 

解题思路:

如此神奇的构图!

首先建立源点汇点S = 0,T = n + 1,然后对所有i∈[0,n]建立一条(i,i+1)容量为k花费为0的边,使所有的点一开始都形成一个串.

那么对于一个覆盖l->r,用费为w,建立一条(l,r+1)容量为1花费为w的边,这怎么理解呢?

原来在没有覆盖区间的网络里面形成的那一串所有点的流量都是k,也就是说从S出发每次流出一点流量的话，那么每个点后面的边最后都会被经过k次，也就相当于k次覆盖。流过(i,i+1)边的流量为m，那么说明i点最后被覆盖了k-m次，也就是说可以把这个边的流量看做i点的无效流流量。 

但原来的边中没有花费，所以我们就要用有花费的边，替代原来的区间，现在建立了一条(l,r+1)的边流量为1,如果有流量流入这条边,那么相当于[l,r]的点后面的边的一点流量跑到(l.r+1)上面去了。

换句话说就是选取k次，每次选取若干个区间不重叠的使得他们的值和最大，那么选取k次，每个点最多也就被选到k次。我们可以把它考虑成单向的，所以可以转化成跑k次流的最长路径和。

那么接下来就是跑最大费用流就ok了.

 

#include<bits/stdc++.h>
#define inf 0x3f3f3f3f
using namespace std;
typedef long long ll;
const int mx = 1100;
int n,m,K,W,head[mx],tot,S,T;
struct node
{
    int y,flow,c;
    int nxt;
}edge[mx<<1];
int u[mx],v[mx],w[mx],cnt[mx];
void AddEdge(int x,int y,int f,int c)
{
    edge[tot] = {y,f,c,head[x]};
    head[x] = tot++;
    edge[tot] = {x,0,-c,head[y]};
    head[y] = tot++;
}
int dis[mx],pre[mx];
bool vis[mx];
bool spfa()
{
    queue <int> que;
    que.push(S);
    for(int i=0;i<=T;i++) dis[i] = pre[i] = -1;
    dis[S] = 0;
    while(!que.empty()){
        int no = que.front();
        que.pop();
        vis[no] = 0;
        for(int i=head[no];~i;i=edge[i].nxt)
        {
            int u = edge[i].y;
            if(!edge[i].flow) continue;
            if(dis[u]<dis[no]+edge[i].c){
                dis[u] = dis[no] + edge[i].c;
                pre[u] = i;
                if(!vis[u]){
                    que.push(u);
                    vis[u] = 1;
                }
            }
        }
    }
    return ~dis[T];
}
int maxcost()
{
    int ans = 0;
    while(spfa()){
    	int mi = inf;
        for(int i=pre[T];~i;i=pre[edge[i^1].y]) mi = min(edge[i].flow,mi);
        for(int i=pre[T];~i;i=pre[edge[i^1].y])
		{
			edge[i].flow -= mi;
			edge[i^1].flow += mi;    
		} 
        ans += dis[T];
    }
    return ans;
}
int main()
{    
    int t;
    scanf("%d",&t);
    while(t--){
        scanf("%d%d%d",&n,&K,&m);
        memset(head,-1,sizeof(head));
        tot = 0;
        int l = 1,r,top = 0;
        cnt[top++] = 0;
        cnt[top++] = n+1;
        for(int i=1;i<=m;i++){
        	scanf("%d%d%d",v+i,u+i,w+i);
        	cnt[top++] = v[i];
        	cnt[top++] = u[i];
		}
		sort(cnt,cnt+top);
		for(int i=1;i<top;i++) if(cnt[i]!=cnt[i-1]) cnt[l++] = cnt[i];
		for(int i=0;i<l-1;i++) AddEdge(i,i+1,K,0);
		for(int i=1;i<=m;i++)
		{
			int L = lower_bound(cnt,cnt+l,v[i])-cnt;
			int R = lower_bound(cnt,cnt+l,u[i])-cnt;
			AddEdge(L,R+1,1,w[i]);
		}
		S = 0,T = l-1;
		printf("%d\n",maxcost());
    }
    return 0; 
}