分点新技能+网络流

当某个点在属性不同的决策时不同时,我们考虑分点,这不失为一种好的解决问题的方法
CTSC家园
由于人类对自然的疯狂破坏,人们意识到在大约2300年之后,地球不能再居住了,于是在月球上建立了新的绿地,以便在需要时移民。令人意想不到的是,2177年冬由于未知的原因,地球环境发生了连锁崩溃,人类必须在最短的时间内迁往月球。
现有n个太空站处于地球与月球之间(编号1..n),m艘公共交通太空船在其中来回穿梭,每个太空站Si可容纳无限的人,每艘太空船pi只可容纳Hpi人。对于每一艘太空船pi,将周期性地停靠一系列的太空站(Si1,Si2…Sir),如:(1,3,4)表示停靠太空站1 3 4 1 3 4 1 3 4 …。 任一艘太空船从任一个太空站驶往另一个任意的太空站耗时为1。人只能在太空船停靠太空站(或地球、月球)时上船或下船。初始时的人全在地球上,太空船全在初始站(太空船pi处于Si1),目标是让所有的人尽快地全部转移到月球上。

输入格式:
文件第一行为三个正整数 n(太空站个数)、 m(太空船个数)、 k(需要运送的地球上的人的个数),其中 1<=m<=13, 1<=n<=20, 1<=k<=50。
接下来的n行给出了太空船的信息,第i+1行说明太空船pi,此行第一个数表示pi可容纳的人数Hpi,第二个数表示pi停靠一个周期的太空站个数r,1<=r<=n+2, 随后r个数便是停靠的太空站的编号(Si1,Si2,…,Sir), 地球用0表示,月球用-1表示。0时刻时,所有太空船都在初始站,随后开始运行,在时刻1,2,3…等正点时刻各艘太空船停靠相应的太空站,即人只有在0,1,2…等正点时刻才能上下太空船。
输出格式:
文件只有一个数,若问题有解,输出完成全部人员安全转移的时刻,否则输出0。
样例输入:
2 2 1
1 3 0 1 2
1 3 1 2 –1
样例输出:
5
数据范围:
1<=m<=13, 1<=n<=20, 1<=k<=50。
时间限制:
1S
空间限制:
256M

//对于那些不断变化的边,我们用分点的方法来解决,对于那些求那些求最小可得时间的问题
//我们可以进行二分答案,而对于那些可以利用前面状态的,我们可以选择递推(更优)。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<queue>
using namespace std;
#define N 200001
#define INF 2e9
int fir[N],to[N*10],ne[N*10],C[N*10],cnt=1;
int n,m,k;
int h[N],r[N],S[1001][1001],ss,tt;
queue<int>q;
bool inq[N];
int tot,dis[N];
void add(int x,int y,int z)
{
    to[++cnt]=y;ne[cnt]=fir[x];fir[x]=cnt;C[cnt]=z;
    to[++cnt]=x;ne[cnt]=fir[y];fir[y]=cnt;C[cnt]=0;
}

bool BFS(int s,int t)
{
//  for(int i=1;i<=tot;i++) inq[i]=1;
    while(!q.empty()) q.pop();
    for(int i=1;i<=tot;i++) dis[i]=-1;
    q.push(s);dis[s]=0;inq[s]=1;
    while(!q.empty())
    {
        int i=q.front();q.pop();
        for(int p=fir[i];p;p=ne[p])
        if(C[p])
        {
            int V=to[p];
            if(dis[V]==-1)
            {
                dis[V]=dis[i]+1;
                q.push(V);
            }
        }
    }
    return (dis[t]!=-1);
}
int dfs(int now,int t,int fl)
{
    int res=0;
    if(now==t||fl==0) return fl;
    for(int i=fir[now];i;i=ne[i])
    if(C[i])
    {
        int V=to[i];
        if(dis[V]==dis[now]+1&&C[i])
        {
            int f=dfs(V,t,min(fl,C[i]));
            if(f)
            {
                fl-=f;
                C[i]-=f;
                C[i^1]+=f;
                res+=f;
            }
            if(!fl) break;
        }
    }
    return res;
}
int dinic(int s,int t)
{
    int res=0;
    while(BFS(s,t))
    {
        int fl=dfs(s,t,2e9);
        if(!fl) return res;
        res+=fl;
    }
    return res;
}
int s,t;
int ans;
int main()
{
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&h[i],&r[i]);
        for(int j=1;j<=r[i];j++)
        {
            scanf("%d",&S[i][j]);
            S[i][j]+=2;
//          cout<<S[i][j]<<" ";
        }
//      cout<<endl;
    }
    tot=2;
    s=1;t=2;
//  add(s,4,INF);
//  add(3,t,INF);
    n+=2;
    for(int i=1;i<=n;i++) tot++;
    for(ans=1;k>0;ans++)
    {
        if(ans>100)
        {
            ans=2;
            break;
        }
        for(int i=1;i<=n;i++)
        {
            tot++;
            add(tot-n,tot,INF);
            if(i==1) add(tot,t,INF);
//          if(i==1) cout<<tot<<" ";
            if(i==2) add(s,tot,INF);
        }
        if(ans>1)
        for(int i=1;i<=m;i++)
        {
            int dd = (ans % r[i]);if (dd==0) dd=r[i];
            int d = dd - 1;if (d==0) d=r[i];
            dd=S[i][dd];d=S[i][d];
            add((ans-1)*n+d+2,ans*n+dd+2,h[i]);
//          cout<<(ans-1)*n+d+2<<" "<<ans*n+dd+2<<endl;
        }
        k-=dinic(s,t);
    }
    cout<<ans-2;
}
/*
2 2 1                 
1 3 0 1 2
1 3 1 2 -1
*/
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值