題目鏈接
題意：有N天，M個妹紙，接下來是一行共M個數，表示M個妹紙要求你在N天內總共給他們拍攝至少Gi個照片。然後有N天，每天有個Ci和Di，表示今天有Ci個妹紙要拍攝，但是今天最多拍攝Di張照片，然後是Ci個妹紙，第一個是妹紙的編號，0～M-1，然後是L～R，表示這個妹紙要求你最少拍L張照片，最多拍R張。
解
這道題，由於我ISAP的退出的頂端高度沒有寫好，竟然WA了幾十次，加上一整個晚上的自閉。
終於，我給我的頂端加了1，過了。這裏ISAP的高度必須是結點數+1，而又有結點數是從0開始的，所以，我們必須將更新更新爲0～node+2，爲什麼是“node + 2”，因爲node + 2，是第一個不可行的高度。必須初始化到node + 2。
關於ISAP的方面就講那麼多了，接下去是關於這道題的解題思路了。
首先，每個妹紙至少總共拍Gi張照片，也就是確定了下限是Gi，所以在這裏要先處理下限；
再則，每天至少0張，至多Di張照片，我們得知了上限；
然後，每個妹紙每天拍Li～Ri張照片，說明的是有上下限。
於是，根據上述幾則條例，我們可以確定如何建圖，然後就是直接跑一個最大流即可，最大流*2是否等於|入流 - 出流|總和，是判斷有無解的標準。
有解的時候，我們刪除超級源匯點，將原來的源點和匯點變成現在的源匯點，再跑一次最大流，即可得知每條邊的流量，加上下限，也就是每條邊的實際流了。
#include <iostream>
#include <cstdio>
#include <cmath>
#include <string>
#include <cstring>
#include <algorithm>
#include <limits>
#include <vector>
#include <stack>
#include <queue>
#include <set>
#include <map>
#include <bitset>
#include <unordered_map>
#include <unordered_set>
#define lowbit(x) ( x&(-x) )
#define pi 3.141592653589793
#define e 2.718281828459045
#define INF 0x7f7f7f7f
#define HalF (l + r)>>1
#define lsn rt<<1
#define rsn rt<<1|1
#define Lson lsn, l, mid
#define Rson rsn, mid+1, r
#define QL Lson, ql, qr
#define QR Rson, ql, qr
#define myself rt, l, r
using namespace std;
typedef unsigned long long ull;
typedef unsigned int uit;
typedef long long ll;
const int maxN = 1505, maxM = 8e5 + 7;
int N, M, head[maxN], cnt, du[maxN], _Index, line[maxM], fl[maxM], sum;
struct Eddge
{
int nex, to, flow;
Eddge(int a=-1, int b=0, int c=0):nex(a), to(b), flow(c) {}
}edge[maxM];
inline void addEddge(int u, int v, int f)
{
edge[cnt] = Eddge(head[u], v, f);
head[u] = cnt++;
}
inline void _add(int u, int v, int f) { addEddge(u, v, f); addEddge(v, u, 0); }
struct ISAP
{
int S, T, gap[maxN], cur[maxN], deep[maxN], que[maxN], ql, qr;
inline void init()
{
for(int i=0; i<=T + 2; i++)
{
gap[i] = deep[i] = 0;
cur[i] = head[i];
}
++gap[deep[T] = 1];
que[ql = qr = 1] = T;
while(ql <= qr)
{
int u = que[ql ++];
for(int i=head[u], v; ~i; i=edge[i].nex)
{
v = edge[i].to;
if(!deep[v]) { ++gap[deep[v] = deep[u] + 1]; que[++qr] = v; }
}
}
}
inline int aug(int u, int Flow)
{
if(u == T) return Flow;
int flow = 0;
for(int &i = cur[u], v; ~i; i=edge[i].nex)
{
v = edge[i].to;
if(deep[u] == deep[v] + 1)
{
int tmp = aug(v, min(edge[i].flow, Flow));
flow += tmp; Flow -= tmp; edge[i].flow -= tmp; edge[i ^ 1].flow += tmp;
if(!Flow) return flow;
}
}
if(!(--gap[deep[u]])) deep[S] = T + 2;
++gap[++deep[u]]; cur[u] = head[u];
return flow;
}
inline int Max_Flow()
{
init();
int ret = aug(S, INF);
while(deep[S] <= T + 1) ret += aug(S, INF);
return ret;
}
} mf;
inline void init()
{
cnt = 0; mf.S = N + M + 2; mf.T = N + M + 3; _Index = 0; sum = 0;
for(int i=0; i<=N + M + 3; i++) { head[i] = -1; du[i] = 0; }
}
int main()
{
while(scanf("%d%d", &N, &M) != EOF)
{
init();
int s = 0, t = N + M + 1;
for(int i=1, Gi; i<=M; i++)
{
scanf("%d", &Gi);
_add(i + N, t, INF);
du[i + N] -= Gi; du[t] += Gi;
}
for(int day = 1, Ci, Di; day <= N; day++)
{
scanf("%d%d", &Ci, &Di);
_add(s, day, Di);
for(int i=1, id, l, r; i<=Ci; i++)
{
scanf("%d%d%d", &id, &l, &r); id += N + 1;
_add(day, id, r - l);
du[day] -= l; du[id] += l;
line[++_Index] = cnt - 1; fl[_Index] = l;
}
}
for(int i=0; i<=t; i++)
{
if(du[i] > 0)
{
_add(mf.S, i, du[i]);
sum += du[i];
}
else if(du[i] < 0)
{
_add(i, mf.T, -du[i]);
sum -= du[i];
}
}
_add(t, s, INF);
if(mf.Max_Flow() * 2 != sum) { printf("-1\n\n"); continue; }
head[mf.S] = head[mf.T] = -1;
mf.S = s; mf.T = t;
printf("%d\n", mf.Max_Flow());
for(int i=1; i<=_Index; i++)
{
printf("%d\n", fl[i] + edge[line[i]].flow);
}
printf("\n");
}
return 0;
}
