poj3680 Intervals (费用流)

建图（(x,y,c,l)表示x到y，费用c，流量l）

(S,1,0,K)

(i,i+1,0,K) 这个边上的流量，表示i还可以被覆盖的次数

(N,T,0,K)

(i,j,w,1)对于权值为w的区间[i,j]

然后跑最大费用最大流

因为没有负权值，所以肯定尽量跑满

 #include<cstring>
 #include<cstdio>
 #include<algorithm>
 #include<queue>
 #define CLR(a,x) memset(a,x,sizeof(a))
 #define MP make_pair
 using namespace std;
 typedef long long ll;
 typedef unsigned long long ull;
 typedef pair<int,int> pa;
 const int maxn=,maxp=,inf=1e9;

 inline ll rd(){
     ll x=;char c=getchar();int neg=;
     while(c<''||c>''){if(c=='-') neg=-;c=getchar();}
     while(c>=''&&c<='') x=x*+c-'',c=getchar();
     return x*neg;
 }

 int l[maxn],r[maxn],w[maxn],N,M,K,tmp[maxp],S=,T=;
 struct Edge{
     int b,l,ne,c;
 }eg[maxp*];
 int egh[maxp],ect=;
 int dis[maxp],fae[maxp];
 bool flag[maxp];

 inline void adeg(int a,int b,int c,int l){
     // printf("~%d %d %d %d\n",a,b,c,l);
     eg[++ect].b=b,eg[ect].l=l,eg[ect].c=c,eg[ect].ne=egh[a],egh[a]=ect;
     eg[++ect].b=a,eg[ect].l=,eg[ect].c=-c,eg[ect].ne=egh[b],egh[b]=ect;
 }

 queue<int> q;
 inline bool spfa(){
     CLR(dis,-);dis[S]=;
     CLR(fae,);q.push(S);
     while(!q.empty()){
         int p=q.front();q.pop();
         // printf("~%d %d\n",p,dis[S]);
         flag[p]=;
         for(int i=egh[p];i;i=eg[i].ne){
             int b=eg[i].b;if(!eg[i].l) continue;
             // printf("!!%d %d %d\n",b,eg[i].c,eg[i].l);
             if(dis[b]<dis[p]+eg[i].c){
                 dis[b]=dis[p]+eg[i].c;
                 fae[b]=i;
                 if(!flag[b]) q.push(b),flag[b]=;
             }
         }
     }return dis[T]>=;
 }

 int main(){
     //freopen("","r",stdin);
     int i,j,k;
     for(int t=rd();t;t--){
         N=rd(),K=rd();
         for(i=;i<=N;i++){
             tmp[i]=l[i]=rd(),tmp[i+N]=r[i]=rd(),w[i]=rd();
         }sort(tmp+,tmp+N+N+);
         int M=unique(tmp+,tmp+N+N+)-tmp-;
         CLR(egh,);ect=;
         for(i=;i<=N;i++){
             l[i]=lower_bound(tmp+,tmp+M+,l[i])-tmp;
             r[i]=lower_bound(tmp+,tmp+M+,r[i])-tmp;
             adeg(l[i],r[i],w[i],);
         }
         for(i=;i<M;i++){
             adeg(i,i+,,inf);
         }adeg(M,T,,K);adeg(S,,,K);
         int ans=;
         while(spfa()){
             int mi=inf;
             for(i=T;i;i=eg[fae[i]^].b){
                 if(fae[i]) mi=min(mi,eg[fae[i]].l);
             }
             for(i=T;i;i=eg[fae[i]^].b){
                 eg[fae[i]].l-=mi,eg[fae[i]^].l+=mi;
             }
             ans+=mi*dis[T];
         }
         printf("%d\n",ans);
     }
     return ;
 }