洛谷 P4043 [AHOI2014/JSOI2014]支线剧情 有上下界最小费用可行流

题目背景
宅男JYY非常喜欢玩RPG游戏，比如仙剑，轩辕剑等等。不过JYY喜欢的并不是战斗场景，而是类似电视剧一般的充满恩怨情仇的剧情。这些游戏往往都有很多的支线剧情，现在JYY想花费最少的时间看完所有的支线剧情。

题目描述
JYY现在所玩的RPG游戏中，一共有N个剧情点，由1到N编号，第i个剧情点可以根据JYY的不同的选择，而经过不同的支线剧情，前往Ki种不同的新的剧情点。当然如果为0，则说明i号剧情点是游戏的一个结局了。
JYY观看一个支线剧情需要一定的时间。JYY一开始处在1号剧情点，也就是游戏的开始。显然任何一个剧情点都是从1号剧情点可达的。此外，随着游戏的进行，剧情是不可逆的。所以游戏保证从任意剧情点出发，都不能再回到这个剧情点。由于JYY过度使用修改器，导致游戏的“存档”和“读档”功能损坏了，
所以JYY要想回到之前的剧情点，唯一的方法就是退出当前游戏，并开始新的游戏，也就是回到1号剧情点。JYY可以在任何时刻退出游戏并重新开始。不断开始新的游戏重复观看已经看过的剧情是很痛苦，JYY希望花费最少的时间，看完所有不同的支线剧情。

输入输出格式

输入格式：
输入格式： 输入一行包含一个正整数N。
接下来N行，第i行为i号剧情点的信息；
第一个整数为，接下来个整数对，Bij和Tij，表示从剧情点i可以前往剧情点，并且观看这段支线剧情需要花费的时间。

输出格式：
输出一行包含一个整数，表示JYY看完所有支线剧情所需要的最少时间。

输入输出样例

输入样例#1：
6
2 2 1 3 2
2 4 3 5 4
2 5 5 6 6
0
0
0
输出样例#1：
24

说明
JYY需要重新开始3次游戏，加上一开始的一次游戏，4次游戏的进程是

样例解释

1->2->4，1->2->5，1->3->5和1->3->6。
对于100%的数据满足N<=300,0<=Ki<=50,1<=Tij<=300,Sigma(Ki)<=5000

分析：
显然就是一个DAG的最小路径覆盖，但是边有费用，就变成了最小费用路径覆盖（注意路径可能多了，这里只追求最小代价）。
和上下界最小流一样做就好了。记录$deg$表示入流的下界-出流的下界，然后如果$deg>0$，那么SS连向这个点，否则这个点连向TT。
然后直接在这个网络跑最小费用最大流，然后加上所有下界*这条边的费用，这些代价是没有算进去的。上下界最小流的$(T,S)$边包含了这些流量。
跑完之后不需要汇点退流，因为这样算的是最小费用最小流，一开始提到了不一定要路径最少。

代码：

// luogu-judger-enable-o2
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
#include <queue>

const int maxn=307;
const int maxe=2e4+7;
const int inf=0x3f3f3f3f;

using namespace std;

int n,m,x,y,ans,cnt,s,t;
int deg[maxn],ls[maxn],dis[maxn],vis[maxn],pre[maxn];

struct edge{
    int x,y,w,c,op,next;
}g[maxe];

queue <int> q;

void add(int x,int y,int w,int c)
{
    g[++cnt]=(edge){x,y,w,c,cnt+1,ls[x]};
    ls[x]=cnt;
    g[++cnt]=(edge){y,x,0,-c,cnt-1,ls[y]};
    ls[y]=cnt;
}

bool spfa()
{
    while (!q.empty()) q.pop();
    memset(dis,inf,sizeof(dis));
    memset(vis,0,sizeof(vis));
    dis[s]=0;
    q.push(s);
    vis[s]=1;
    while (!q.empty())
    {
        int x=q.front();
        q.pop();
        for (int i=ls[x];i>0;i=g[i].next)
        {
            int y=g[i].y;
            if ((g[i].w) && (dis[y]>dis[x]+g[i].c))
            {
                dis[y]=dis[x]+g[i].c;
                pre[y]=i;
                if (!vis[y])
                {
                    vis[y]=1;
                    q.push(y);
                }
            }
        }
        vis[x]=0;
    }
    return (dis[t]!=inf);
}

void mcf()
{
    int i=t,x=1e9;
    while (i!=s)
    {
        x=min(x,g[pre[i]].w);
        i=g[pre[i]].x;
    }
    i=t;
    while (i!=s)
    {
        g[pre[i]].w-=x;
        g[g[pre[i]].op].w+=x;
        ans+=g[pre[i]].c*x;
        i=g[pre[i]].x;
    }
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        scanf("%d",&m);
        for (int j=1;j<=m;j++)
        {
            scanf("%d%d",&x,&y);
            add(i,x,inf,y);
            ans+=y;
            deg[i]--;
            deg[x]++;
        }
    }
    int S=1,T=n+2,SS=n+3,TT=n+4;
    for (int i=1;i<=n;i++) add(i,T,inf,0);
    for (int i=1;i<=n;i++)
    {
        if (deg[i]>0) add(SS,i,deg[i],0);
        else if (deg[i]<0) add(i,TT,-deg[i],0);
    }
    add(T,S,inf,0);
    s=SS,t=TT;
    while (spfa()) mcf();
    printf("%d\n",ans);
}