POJ 3680 <离散化+最小费用最大流模版>

给一个很大的区间，然后给一些区间的权重为wi，要求整个区间的所有点被覆盖不超过K个区间，问应该选哪些区间使得权重最大。

区间很大但是数量不多，那么就离散化就好。

int getid(int x){
	return lower_bound(q2.begin(),q2.end(),x)-q2.begin()+1;//返回下标+1 
}

得到下标。

sort(q2.begin(),q2.end());
        q2.erase(unique(q2.begin(),q2.end()),q2.end());

可转化成一个图，区间i->i+1的费用为0，容量无穷，现在就成了找最大费用最大流。

可以把权值变成负数，然后刚开始这条边就满流，最小费用就成了负数最小，结果加个负号就成了最大的费用。

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <algorithm>
#include <queue>
#define V 1100
#define E 1000100
#define inf 99999999
using namespace std;
int vis[V]; //V为点的数量 
int dist[V];
int pre[V];

struct Edge{
    int u,v,c,cost,next;
}edge[E];
int head[V],cnt;

void init(){
    cnt=0;
    memset(head,-1,sizeof(head));
}
void addedge(int u,int v,int c,int cost)
{
    edge[cnt].u=u;edge[cnt].v=v;edge[cnt].cost=cost;
    edge[cnt].c=c;edge[cnt].next=head[u];head[u]=cnt++;

    edge[cnt].u=v;edge[cnt].v=u;edge[cnt].cost=-cost;
    edge[cnt].c=0;edge[cnt].next=head[v];head[v]=cnt++;
}

bool spfa(int begin,int end){
    int u,v;
    queue<int> q;
    for(int i=0;i<=end+2;i++){
        pre[i]=-1;
        vis[i]=0;
        dist[i]=inf;
    }
    vis[begin]=1;
    dist[begin]=0;
    q.push(begin);
    while(!q.empty()){
        u=q.front();
        q.pop();
       // cout << "find " << endl;
        vis[u]=0;
        for(int i=head[u];i!=-1;i=edge[i].next){
            if(edge[i].c>0){
                v=edge[i].v;
                if(dist[v]>dist[u]+edge[i].cost){
                    dist[v]=dist[u]+edge[i].cost;
                    pre[v]=i;
                    if(!vis[v]){
                        vis[v]=true;
                        q.push(v);
                    }
                }
            }
        }
    }
    return dist[end]!=inf;
}

int MCMF(int begin,int end){
    int ans=0,flow;
    int flow_sum=0;
    while(spfa(begin,end)){
        flow=inf;
        //cout << "一次" << endl;
        for(int i=pre[end];i!=-1;i=pre[edge[i].u])
            if(edge[i].c<flow)
                flow=edge[i].c;
        for(int i=pre[end];i!=-1;i=pre[edge[i].u]){
            edge[i].c-=flow;
            edge[i^1].c+=flow;
        }
        ans+=dist[end];
        flow_sum += flow;
    }
    //cout << flow_sum << endl;
    return ans;
}
struct ttt{
	int l,r,w;
};
ttt q1[205];
vector<int>q2;
int getid(int x){
	return lower_bound(q2.begin(),q2.end(),x)-q2.begin()+1;//返回下标+1 
}
int main()
{
   // freopen("in.txt","r",stdin);
    int n,m,i,a,b,c,t1,t2,t3;
    int T;
    cin >> T;
    while(T--){
        scanf("%d %d",&n,&m);
		init();
		q2.clear();
        for(i=1;i<=n;i++){
        	scanf("%d %d %d",&t1,&t2,&t3);
        	q1[i].l=t1;q1[i].r=t2;q1[i].w=t3;
			q2.push_back(t1);
			q2.push_back(t2);
        }
        sort(q2.begin(),q2.end());
        q2.erase(unique(q2.begin(),q2.end()),q2.end());
		addedge(0,1,m,0);//下标0-1 
		addedge(q2.size(),q2.size()+1,m,0);
		for(i=1;i<q2.size();i++){
			addedge(i,i+1,inf,0);//离散化后这里应该为下标， 
		}
		for(i=1;i<=n;i++){
			int u=getid(q1[i].l);
			int v=getid(q1[i].r);
			addedge(u,v,1,-q1[i].w);
			edge[cnt-2].c=1; //不能为++，因为初始可能不为1 
		}
		printf("%d\n",-MCMF(0,q2.size()+1));
    }
    return 0;
}