http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3229

題目大意：射命丸文要給幻想鄉的居民照相，共照n天m個人，每天射命丸文照相數不多於d個，且一個人n天一共被拍的照片不能少於g個，且每天可照的人有限制，且這些人今天照的相片必須在[l,r]以內，求是否有可行解，如果有則輸出最多照片數，並且輸出每天每個可以被照的人的被照的照片數。

————————————————————————

https://www.cnblogs.com/kane0526/archive/2013/04/05/3001108.html

按照這篇博客的順序先刷了這道題。

（話說射命丸文怎麽求成屌絲了……）

建圖很顯然，源點到每天連一條[0,d]，每天到每天能連的人建[l,r]，每個人到匯點連[g,INF]。

現在就是有源匯上下界網絡流的題了。

按照無源匯上下界網絡流的想法（可參考我的上一篇博客ZOJ2314）是很容易建的，有源匯的話我們就手動把它變成無源匯的即可。

那麽我們從匯點向源點連一條[0,INF]不就可以了嗎？

接下來就是我們建立一個超級源匯點，跑一遍無源匯上下界網絡流即可。

但！是！還！沒！完！我們跑完的流量實際上只是流向超級匯點的流量，它只是流向匯點的流量的一部分，在殘余網絡裏還可能存在一些被卡在管子裏的流量無法流向超級匯點但可以流向匯點。

所以我們機智的刪掉了超級源匯點和一些用於無源匯上下界網絡流的邊，再跑一邊從源點到匯點的最大流即可。

至於每個人每天的照片，那自然是每天到每個人的反邊的容量+下界了！

#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
const int N=1501;
const int M=740001;
const int INF=1e9;
inline int
 
 read(){
    int X=0,w=0;char ch=0;
    while(!isdigit(ch)){w|=ch==‘-‘;ch=getchar();}
    while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
    return w?-X:X;
}
struct node{
    int nxt;
    int to;
    int w;
}edge[M];
int head[N],du[N],id[366][1001],low[366][1001],cnt=-1;
int S,T;
void add(int u,int v,int w){
    cnt++;
    edge[cnt].to=v;
    edge[cnt].w=w;
    edge[cnt].nxt=head[u];
    head[u]=cnt;
}
int lev[N],cur[N],dui[N];
bool bfs(int m){
    int r=0;
    for(int i=1;i<=m;i++){
    lev[i]=-1;
    cur[i]=head[i];
    }
    dui[0]=S,lev[S]=0;
    int u,v;
    for(int l=0;l<=r;l++){
    u=dui[l];
    for(int e=head[u];e!=-1;e=edge[e].nxt){
        v=edge[e].to;
        if(edge[e].w>0&&lev[v]==-1){ 
        lev[v]=lev[u]+1;
        r++;
        dui[r]=v; 
        if(v==T)return 1; 
        }
    }
    }
    return 0;
}
int dinic(int u,int flow,int m){
    if(u==m)return flow;
    int res=0,delta;
    for(int &e=cur[u];e!=-1;e=edge[e].nxt){
    int v=edge[e].to;
    if(edge[e].w>0&&lev[u]<lev[v]){ 
        delta=dinic(v,min(edge[e].w,flow-res),m); 
        if(delta>0){
        edge[e].w-=delta;
        edge[e^1].w+=delta;
        res+=delta;
        if(res==flow)break; 
        }
    }
    }
    if(res!=flow)lev[u]=-1;
    return res;
}
inline void init(){
    memset(head,-1,sizeof(head));
    memset(du,0,sizeof(du));
    memset(id,0,sizeof(id));
    cnt=-1;
    return;
}
int main(){
    int n,m;
    while(scanf("%d%d",&n,&m)!=EOF&&n){
    int st=m+n+1,ed=m+n+2;
    init();
    for(int i=1;i<=m;i++){
        int g=read();
        add(i,ed,INF-g);
        add(ed,i,0);
        du[i]-=g;
        du[ed]+=g;
    }
    for(int i=m+1;i<=m+n;i++){
        int c=read(),d=read();
        add(st,i,d);
        add(i,st,0);
        for(int j=1;j<=c;j++){
        int t=read()+1,l=read(),r=read();
        add(i,t,r-l);
        add(t,i,0);
        du[i]-=l;
        du[t]+=l;
        id[i-m][t]=cnt;
        low[i-m][t]=l;
        }
    }
    int h1=head[st],h2=head[ed];
    add(ed,st,INF);add(st,ed,0);
    S=ed+1;T=S+1;
    int ans=0,full=0;
    for(int i=1;i<=m+n+2;i++){
        if(du[i]>0){
        add(S,i,du[i]);
        add(i,S,0);
        full+=du[i];
        }else if(du[i]<0){
        add(i,T,-du[i]);
        add(T,i,0);
        }
    }
    while(bfs(T)==1)ans+=dinic(S,INF,T);
    if(ans!=full)puts("-1");
    else{
        head[S]=head[T]=-1;
        head[st]=h1;head[ed]=h2;
        S=st;T=ed;
        while(bfs(ed)==1)ans+=dinic(S,INF,T);
        printf("%d\n",ans);
        for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            if(id[i][j]){
            printf("%d\n",edge[id[i][j]].w+low[i][j]);
            }
        }
        }
    }
    puts("");
    }
    return 0;
}

ZOJ3229：Shoot the Bullet——題解