poj3680Intervals【最大权不相交路径】

题意：和nefu495最长k可重区间集问题【最大权不相交路径】网络流24题
一样，只是权值不是线段长度变成了给定值。

做法：一样的……

【建模方法】
方法1
按左端点排序所有区间，把每个区间拆分看做两个顶点<i.a><i.b>，建立附加源S汇T，以及附加顶点S'。
1、连接S到S'一条容量为K，费用为0的有向边。
2、从S'到每个<i.a>连接一条容量为1，费用为0的有向边。
3、从每个<i.b>到T连接一条容量为1，费用为0的有向边。
4、从每个顶点<i.a>到<i.b>连接一条容量为1，费用为区间长度的有向边。
5、对于每个区间i，与它右边的不相交的所有区间j各连一条容量为1，费用为0的有向边。
求最大费用最大流，最大费用流值就是最长k可重区间集的长度。
方法2
离散化所有区间的端点，把每个端点看做一个顶点，建立附加源S汇T。
1、从S到顶点1（最左边顶点）连接一条容量为K，费用为0的有向边。
2、从顶点2N（最右边顶点）到T连接一条容量为K，费用为0的有向边。
3、从顶点i到顶点i+1(i+1<=2N)，连接一条容量为无穷大，费用为0的有向边。
4、对于每个区间[a,b]，从a对应的顶点i到b对应的顶点j连接一条容量为1，费用为区间长度的有向边。
求最大费用最大流，最大费用流值就是最长k可重区间集的长度。
【建模分析】
这个问题可以看做是求K条权之和最大的不想交路径，每条路径为一些不相交的区间序列。由于是最大费用流，两条路径之间一定有一些区间相交，可以看做事相交部分重复了2次，
而K条路经就是最多重复了K次。最简单的想法就是把区间排序后，不相交的区间之间连接一条边，由于每个区间只能用一次，所以要拆点，点内限制流量。如果我们改变一下思路，
把端点作为网络中的顶点，区间恰恰是特定一些端点之间的边，这样建模的复杂度更小。方法1的边数是O(N^2)的，而方法2的边数是O(N)的，可以解决更大规模的问题。

这次用的是方法二，之前一直不理解方法二是如何控制流量的QAQ画个图就懂了



反倒是方法一又不懂了卧槽……(反正方法二更好）

#include<cstdio>
#include<algorithm>
using namespace std;
const int mm=333300;
const int mn=1111;
const int oo=1000000000;
int node,src,dest,edge;
int reach[mm],flow[mm],cost[mm],next[mm];
int head[mn],dis[mn],q[mn],p[mn],a[mn],x[mn],y[mn];
bool vis[mn];
inline void prepare(int _node,int _src,int _dest)
{
node=_node,src=_src,dest=_dest;
for(int i=0;i<node;++i)head[i]=-1;
edge=0;
}
inline void addedge(int u,int v,int f,int c)
{
reach[edge]=v,flow[edge]=f,cost[edge]=c,next[edge]=head[u],head[u]=edge++;
reach[edge]=u,flow[edge]=0,cost[edge]=-c,next[edge]=head[v],head[v]=edge++;
}
bool spfa()
{
int i,u,v,l,r=0,tmp;
for(i=0;i<node;++i)dis[i]=oo;
dis[q[r++]=src]=0;
p[src]=p[dest]=-1;
for(l=0;l!=r;(++l==mn)?l=0:1)
for(i=head[u=q[l]],vis[u]=0;i>=0;i=next[i])
if(flow[i]&&dis[v=reach[i]]>(tmp=dis[u]+cost[i]))
{
dis[v]=tmp;
p[v]=i^1;
if(vis[v])continue;
vis[q[r++]=v]=1;
if(r==mn)r=0;
}
return p[dest]>=0;
}
int SpfaFlow()
{
int i,ret=0,delta;
while(spfa())
{
for(i=p[dest],delta=oo;i>=0;i=p[reach[i]])
if(delta>flow[i^1])delta=flow[i^1];
for(i=p[dest];i>=0;i=p[reach[i]])
flow[i]+=delta,flow[i^1]-=delta;
ret-=delta*dis[dest];
}
return ret;
}
int find(int x,int m)
{
int l=0,r=m;
while(l<r)
{
m=(l+r)>>1;
if(a[m]==x)return m;
if(a[m]>x)r=m-1;
else l=m+1;
}
return l;
}
int val[mn];
int main()
{
// freopen("input.txt","r",stdin);
// freopen("output.txt","w",stdout);
int i,u,v,n,k,m,t;
scanf("%d",&t);
while(t--)
{
scanf("%d%d",&n,&k);
for(m=i=0;i<n;++i)scanf("%d%d%d",&x[i],&y[i],&val[i]),a[m++]=x[i],a[m++]=y[i];
sort(a,a+m);
for(v=0,u=1;u<m;++u)
if(a[u]>a[v])a[++v]=a[u];
m=v+1;
prepare(m+2,0,m+1);
addedge(src,1,k,0);
addedge(m,dest,k,0);
for(u=1;u<m;++u)addedge(u,u+1,oo,0);
for(i=0;i<n;++i)
addedge(find(x[i],m-1)+1,find(y[i],m-1)+1,1,-val[i]);
printf("%d\n",SpfaFlow());
}
return 0;
}