zoukankan      html  css  js  c++  java
  • 【codevs1034】家园——网络流

    由数据范围看是可以用网络流的。。。

    其实还是接近于很裸的网络流,不过是变成了动态建边,图会很大,不过题目数据范围小啊。。。

    枚举时间,自己想一个上界,cyc大佬说是100,那就100吧,到达上界前找到就输出当前时间并退出,超出上界就输出0。

    cyc大佬反复强调tot=1而不是2!!!

    然后就跑啊跑,记得加当前弧优化!

    #include<cstdio>
    #include<iostream>
    #include<cstring>
    using namespace std;
    const int inf=0x3f3f3f3f;
    struct point
    {
        int next,flow,to;
    }e[50010];
    int n,m,k,t,tot=1;
    int p[22],a[22][22],first[4000],cur[4000],d[4000],q[5001];
    void add(int u,int v,int w)
    {
        tot++;e[tot].to=v;e[tot].next=first[u];first[u]=tot;e[tot].flow=w;
        tot++;e[tot].to=u;e[tot].next=first[v];first[v]=tot;e[tot].flow=0;
    }
    bool bfs()
    {
        memset(d,-1,sizeof(d));
        int head=0,tail=1;q[0]=0;d[0]=0;
        while(head!=tail)
        {
            int x=q[head++];if(head>=5000)head=0;
            for(int i=first[x];i;i=e[i].next)
            {
                if(d[e[i].to]==-1&&e[i].flow>0)
                {
                    q[tail++]=e[i].to;
                    if(tail>=5000)tail=0;
                    d[e[i].to]=d[x]+1;
                }
            }
        }
        if(d[t*22+1+n]==-1)return false;
        return true;
    }
    int dfs(int x,int a)
    {
        if(x==t*22+1+n||a==0)return a;
        int flow=0,f;
        for(int& i=cur[x];i;i=e[i].next)
        {
            if(d[e[i].to]==d[x]+1&&(f=dfs(e[i].to,min(a,e[i].flow)))>0)
            {
                e[i].flow-=f;
                e[i^1].flow+=f;
                a-=f;
                flow+=f;
                if(a==0)break;
            }
        }
        return flow;
    }
    int main()
    {
        scanf("%d %d %d",&n,&m,&k);
        for(int i=1;i<=m;i++)
        {
            scanf("%d %d",&p[i],&a[i][0]);
            for(int j=1;j<=a[i][0];j++)
            {
                scanf("%d",&a[i][j]);
                if(a[i][j]==-1)a[i][j]=n+1;
            }
        }
        for(t=1;t<=100;t++)//从1开始枚举时间
        {
            for(int i=0;i<=n+1;i++)
            add((t-1)*22+i,t*22+i,inf);//建边
            for(int i=1;i<=m;i++)
            {
                int u=(t-1)*22+a[i][(t-1)%a[i][0]+1];
                int v=t*22+a[i][t%a[i][0]+1];
                add(u,v,p[i]);//建边
            }
            while(bfs())//增广
            {
                for(int i=0;i<=t*22+n+1;i++)cur[i]=first[i];
                k-=dfs(0,inf);
            }
            if(k<=0){printf("%d",t);return 0;}
        }
        printf("0");
        return 0;
    }
  • 相关阅读:
    更准确的mysql全文索引
    range
    牛顿冷却定律 使用
    Servo: The Embeddable Browser Engine
    Node.js V0.12新特性之性能优化
    Lodash,你正在使用的JavaScript库
    Python on Android
    Microsoft HoloLens 技术解谜(下)
    Microsoft HoloLens 技术解谜(上)
    市售体感设备横评
  • 原文地址:https://www.cnblogs.com/JKAI/p/6842827.html
Copyright © 2011-2022 走看看