分点新技能+网络流

来源:互联网 发布:im短域名 编辑:程序博客网 时间:2024/05/22 07:15

当某个点在属性不同的决策时不同时,我们考虑分点,这不失为一种好的解决问题的方法
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 2e9int 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 21 3 1 2 -1*/
原创粉丝点击