bzoj 4276 [ONTAK2015]Bajtman i Okrągły Robin

http://www.elijahqi.win/archives/3244
Description

有n个强盗，其中第i个强盗会在[a[i],a[i]+1]，[a[i]+1,a[i]+2]，…，[b[i]-1,b[i]]这么多段长度为1时间中选出一个时间进行抢劫，并计划抢走c[i]元。作为保安，你在每一段长度为1的时间内最多只能制止一个强盗，那么你最多可以挽回多少损失呢？

Input

第一行包含一个正整数n(1<=n<=5000)，表示强盗的个数。
接下来n行，每行包含三个正整数a[i],b[i],ci，依次描述每一个强盗。

Output

输出一个整数，即可以挽回的损失的最大值。

Sample Input

4
1 4 40
2 4 10
2 3 30
1 3 20
Sample Output

90
HINT

Source

By Claris

每个人流量1 费用c[i] 每个人都能覆盖一段区间 所以线段树优化建图 暴力费用流即可

#include<queue>
#include<cstdio>
#include<cctype>
#include<cstring>
#include<algorithm>
#define lc (x<<1)
#define rc (x<<1|1)
using namespace std;
inline char gc(){
    static char now[1<<16],*S,*T;
    if (T==S){T=(S=now)+fread(now,1,1<<16,stdin);if(T==S) return EOF;}
    return *S++;
}
inline int read(){
    int x=0,f=1;char ch=gc();
    while(!isdigit(ch)) {if(ch=='-') f=-1;ch=gc();}
    while(isdigit(ch)) x=x*10+ch-'0',ch=gc();
    return x*f;
}
const int inf=0x3f3f3f3f;
const int N=1e4;
struct node{
    int x,y,next,z,c;
}data[N*100];
bool flag[N<<2];int num=1,a[N],b[N],c[N],n;
int dis[N<<2],T,pre[N<<2],path[N<<2],h[N<<2],rr;
inline void insert1(int x,int y,int z,int c){
    data[++num].y=y;data[num].next=h[x];h[x]=num;data[num].z=z;data[num].c=c;data[num].x=x;
    data[++num].y=x;data[num].next=h[y];h[y]=num;data[num].z=0;data[num].c=-c;data[num].x=y;
}
inline bool spfa(){
    queue<int>q;q.push(0);for (int i=0;i<=T;++i) dis[i]=-inf;
    dis[0]=0;flag[0]=1;
    while(!q.empty()){
        int x=q.front();q.pop();flag[x]=0;
        for (int i=h[x];i;i=data[i].next){
            int y=data[i].y,z=data[i].z,c=data[i].c;
            if (dis[x]+c>dis[y]&&z){
                dis[y]=dis[x]+c;pre[y]=x;path[y]=i;
                if (!flag[y]) flag[y]=1,q.push(y);
            }
        }
    }return dis[T]!=-inf;
}
inline void build(int x,int l,int r){
    if (l==r) {insert1(x,T,1,0);return;}int mid=l+r>>1;
    insert1(x,lc,inf,0);insert1(x,rc,inf,0);
    build(lc,l,mid);build(rc,mid+1,r);
}
inline void modify(int x,int l,int r,int l1,int r1,int id){
    if(l1<=l&&r1>=r) {insert1((rr<<2)+id,x,1,0);return;}int mid=l+r>>1;
    if (l1<=mid) modify(lc,l,mid,l1,r1,id);
    if (r1>mid) modify(rc,mid+1,r,l1,r1,id);
}
int main(){
    freopen("bzoj4276.in","r",stdin);
    n=read();rr=0;
    for (int i=1;i<=n;++i) a[i]=read(),b[i]=read()-1,c[i]=read(),rr=max(rr,b[i]);
    T=(rr<<2)+n+1;build(1,1,rr);
    for (int i=1;i<=n;++i){
        insert1(0,(rr<<2)+i,1,c[i]);
        modify(1,1,rr,a[i],b[i],i);
    }int ans=0;
    while(spfa()){
        int minn=inf,now=T;
        while(now) minn=min(minn,data[path[now]].z),now=pre[now];ans+=dis[T]*minn;now=T;
        while(now) data[path[now]].z-=minn,data[path[now]^1].z+=minn,now=pre[now];
    }printf("%d\n",ans);
    return 0;
}