分层图网络流问题,枚举答案,构造网络流判定。
建模方法:
首先判断从地球到月球是否存在一条路线,如果不存在那么无解,否则把每个太空站按照每天拆分成d 个点,<i,d>表示第i
个站第d 天。建立附加源S 汇T,顺序枚举答案Day。
1、对于第Day 天,从S 到<0,Day>连接一条容量为无穷大的有向边。
2、从<-1,Day>到T 连接一条容量为无穷大的有向边。
3、对于第i 个太空船,设第Day-1天在a 处,第Day 天在b 处,从<a,Day-1>到<b,Day>连接一条容量为该太空船容量的有向边。
4、对于第i 个太空站,从<i,Day-1>到<i,Day>连接一条容量为无穷大的有向边。
5、求当前网络最大流,如果最大流量大于等于地球上人数K,停止枚举,当前Day 值就是答案。
建模分析:
我们把网络优化问题转化为枚举答案+可行性判定问题。枚举天数,按天数把图分层,因为乘船每坐一站天数都要增加一,把
太空船航线抽象成图中的一条边,跨图的两层。由于太空船容量有限,边上也要加上容量限制。除了坐船以外,人还可以在某个
空间站等待下一班太空船的到来,所以每个点要与下一层同一点连接一条容量为无穷的边。这样在层限制的图上求出的网络最大
流就是在当前天数以内能够从地面到月球的最多的人数,该人数随天数递增不递减,存在单调性。所以可以枚举答案或二分答案,
用网络流判定。网络流判定问题更适合枚举答案,而不是二分,因为新增一些点和边只需要在原有的基础上增广,不必重新求网
络流。
#include <iostream>//纪念写了一天的代码,有点小问题,实在是懒得写了,就是判断是否连通的时候,我直接就用<365,正确的应该是运用广搜来判断是否连通
#include <cstdio>
#include <cstring>
using namespace std;
const int oo=1e4;
const int mm=1111111;
const int mn=9999;
int node,src,dest,edge;
int ver[mm],flow[mm],next[mm];
int head[mn],work[mn],dis[mn],q[mn];
int n,m;
int P(int vi ,int ti){
return ti * (n + 2) + vi + 1;
}
void prepare(int _node,int _src,int _dest)
{
node=_node,src=_src,dest=_dest;
for(int i=0;i<node;i++) head[i]=-1;
edge=0;
}
void addedge(int u,int v,int c)
{
ver[edge]=v;flow[edge]=c;next[edge]=head[u];head[u]=edge++;
ver[edge]=u;flow[edge]=0;next[edge]=head[v];head[v]=edge++;
}
bool Dinic_bfs()
{
int i,u,v,l,r=0;
for(i=0;i<node;i++)
dis[i]=-1;
dis[q[r++]=src]=0;
for(l=0;l<r;l++)
for(i=head[u=q[l]];i>=0;i=next[i])
if(flow[i]&&dis[v=ver[i]]<0)
{
dis[q[r++]=v]=dis[u]+1;
if(v==dest) return 1;
}
return 0;
}
int Dinic_dfs(int u,int exp)
{
if(u==dest) return exp;
for(int &i=work[u],v,tmp;i>=0;i=next[i])
{
if(flow[i]&&dis[v=ver[i]]==dis[u]+1&&(tmp=Dinic_dfs(v,min(exp,flow[i])))>0)
{
flow[i]-=tmp;
flow[i^1]+=tmp;
// cout<<"u="<<u<<" v="<<v<<" tmp="<<tmp<<endl;
return tmp;
}
}
return 0;
}
int Dinic_flow()
{
int i,ret=0,delta;
while(Dinic_bfs())
{
for(i=0;i<node;i++)
work[i]=head[i];
while(delta=Dinic_dfs(src,oo)) {ret+=delta;}
}
return ret;
}
int main()
{
int k,i,j,tmp,ans;
int hp[110],r[110],rr[110][110];
while(~scanf("%d%d%d",&n,&m,&k))
{
for(i=0;i<m;i++)
{
scanf("%d",&hp[i]);
scanf("%d",&r[i]);
for(j=0;j<r[i];j++){
scanf("%d",&rr[i][j]);
if(rr[i][j]==-1) rr[i][j]=n+1;
}
}
int S=oo-3;
int T=S+1;
prepare(mn,S,T);
ans=0;
addedge(S,P(0,0),oo);
addedge(P(n+1,0),T,oo);
for(i=1,ans=0;ans<k&&i<365;i++)
{
addedge(S,P(0,i),oo);addedge(P(n+1,i),T,oo);
for(j=0;j<n+2;j++)
addedge(P(j,i-1),P(j,i),oo);
for(j=0;j<m;j++)
{
int u=rr[j][(i-1)%r[j]];
int v=rr[j][i%r[j]];
addedge(P(u,i-1),P(v,i),hp[j]);
// cout<<"u="<<P(u,i-1)<<" v="<<P(v,i)<<" hp="<<hp[j]<<endl;
}
ans+=Dinic_flow();//cout<<"ans="<<ans<<endl;
}
if(i<365)
printf("%d\n",i-1);
else
printf("0\n");
}
return 0;
}