您的位置:首页 > 其它

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

2017-06-05 19:48 246 查看
给一个很大的区间,然后给一些区间的权重为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;
}
内容来自用户分享和网络整理,不保证内容的准确性,如有侵权内容,可联系管理员处理 点击这里给我发消息
标签: