【Codevs1422】【网络流】河城荷取

东方幻想乡系列省选模拟赛Stage4。这些题目现在只在Codevs上有出现了别的地方都消失了。

题目描述 Description

在幻想乡，河城荷取是擅长高科技工业的河童。荷取的得意之作除了光学迷彩外，还有震动整个幻想乡的巨型人形『非想天则』。不过由于人形太过巨大，所以为它充能是一件很麻烦的事。人形一共有N个电能池，编号1..N。其中前L个电能池(即编号为1..L的电能池)连接着外部充能接口，而编号为N的电能池连接着动力炉核心。在N个蓄能池之间有M条单向管道，每条管道有一个激活代价cost和电能传输极限limit。当激活度达到某个值时，所以激活代价小于等于这个值的管道都会被激活，但是每一条管道只能够最多传送limit个单位的电能。外部接口到电能池和电能池到动力炉核心的管道传输没有限制并且激活代价为0。现在荷取想往动力炉核心输入至少K个单位的电能，求需要的最小激活度。

输入描述 Input Description

第1行：4个正整数N,M,L, K

　　第2..M行：4个整数，u,v,limit,cost，表示一条由u到v的管道，传输极限limit，激活代价为cost

输出描述 Output Description

第1行：1个整数，表示最小激活代价

样例输入 Sample Input

6 5 3 3

1 4 2 4

2 4 3 5

3 5 4 2

4 6 2 3

5 6 3 4

样例输出 Sample Output

4

数据范围及提示 Data Size & Hint

数据范围

　　对于30%的数据：1 ≤\le L ≤\le N ≤\le 100，0 ≤\le M ≤\le 2,000，1 ≤\le cost ≤\le 10,000

　　对于60%的数据：1 ≤\le L ≤\le N ≤\le 1,000，0 ≤\le M ≤\le 20,000，1 ≤\le cost ≤\le 10,000

　　对于100%的数据：1 ≤\le L ≤\le N ≤\le 2,000，0 ≤\leM ≤\le 80,000，1 ≤\le cost ≤\le 1,000,000

　　对于100%的数据：1 ≤\le limit ≤\le 1,000

　提示

　　样例解释：

　　当激活度为4时，除了(2,4)外其他管道都能够使用。此时能够输入恰好4个单位电能。具体如下：

　　(1,4) 输送2个单位电力

　　(4,6) 输送2个单位电力

　　(3,5) 输送2个单位电力

　　(5,6) 输送2个单位电力

　　

　　注意：

　　保证任意(u,v)都只出现一次。

我还真是傻叉这种题还写这么久。

二分费用，然后一次次做最大流，如果满足最大流值大于k就符合条件把ans改成当前二分的mid。

要注意的是每二分出一个新的mid需要重构图。

我重构图时候top没有清零导致无限RE真是呵呵。。。

我还真是傻叉。

[code]#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
#define MAXN 2010
#define MAXM 200000
#define MAXINT 0x7fffffff
#define MAXX 1000000
using namespace std;
int dis[MAXN];
int num[MAXN];
int s,t;
int top;
int n,m,l,k;
int u[MAXM],v[MAXM],limit[MAXM],cost[MAXM];
struct edge
{
    int f;
    int c;
    int to;
    edge *rev;
    edge *next;
}e[MAXM],*prev[MAXN];
void insert(int u,int v,int c)
{
    e[++top].to=v;
    e[top].c=c;
    e[top].next=prev[u];
    prev[u]=&e[top];
}
void init()
{
    int que[MAXN],head=0,tail=0;
    for (int i=0;i<=n+1;i++) dis[i]=MAXINT,num[i]=0;
    que[tail++]=t;
    dis[t]=0;
    num[0]=1;
    while (head!=tail)
    {
        int x=que[head++];
        for (edge *i=prev[x];i;i=i->next)
        {
            if ((i->rev)&&(i->rev->c==0)||(dis[i->to]<MAXINT)){}
            else
            {
                dis[i->to]=dis[x]+1;
                ++num[dis[i->to]];
                que[tail++]=i->to;
            }
        }
    }
}
int maxflow()
{
    int st=s,ret=0;
    edge *E[MAXN],*rep[MAXN];
    for (int i=0;i<=n+1;i++) E[i]=prev[i];
    while (dis[s]<n+2)
    {
        if (st==t)
        {
            int delta=MAXINT;
            for (int i=s;i!=t;i=E[i]->to) delta=min(delta,E[i]->c);
            for (int i=s;i!=t;i=E[i]->to) 
            {
                E[i]->c-=delta;
                E[i]->f+=delta;
                E[i]->rev->c+=delta;
                E[i]->rev->f-=delta;
            }
            ret+=delta;
            st=s;
        }
        edge *i;
        for (i=E[st];i;i=i->next)
            if (i->c>0&&dis[st]==dis[i->to]+1) break;
        if (i)
        {
            E[st]=i;
            rep[i->to]=i->rev;
            st=i->to;
        }
        else
        {
            if ((--num[dis[st]])==0) break;
            E[st]=prev[st];
            int x=n+2;
            for (edge *t=prev[st];t;t=t->next)
                if (t->c>0) x=min(x,dis[t->to]);
            dis[st]=x+1;
            ++num[dis[st]];
            if (st!=s) st=rep[st]->to;
        }
    }
    return ret;
}
void rebuild(int x)
{
    top=0;
    memset(e,0,sizeof(e));
    memset(prev,0,sizeof(prev));
    for (int i=1;i<=l;i++) insert(0,i,MAXINT),insert(i,0,0),prev[0]->rev=prev[i],prev[i]->rev=prev[0];
    for (int i=1;i<=m;i++)
        if (x>=cost[i]) insert(u[i],v[i],limit[i]),insert(v[i],u[i],0),prev[u[i]]->rev=prev[v[i]],prev[v[i]]->rev=prev[u[i]];
    insert(n,n+1,MAXINT);insert(n+1,n,0);prev
->rev=prev[n+1];prev[n+1]->rev=prev
;
}
int main()
{
    scanf("%d%d%d%d",&n,&m,&l,&k);
    s=0;
    t=n+1;
    for (int i=1;i<=m;i++) scanf("%d%d%d%d",&u[i],&v[i],&limit[i],&cost[i]);
    int l=0,r=MAXX;
    int minn=MAXINT;
    while (l<=r)
    {
        int mid=(l+r)>>1;
        rebuild(mid);
        init();
        int temp=maxflow();
        if (temp>=k)
        {
            minn=mid;
            r=mid-1;
        }
        else
            l=mid+1;
    }
    printf("%d\n",minn);
}