星际转移问题

分析(引用BYvoid 大牛的分析):

分层图网络流问题,枚举答案,构造网络流判定。

建模方法:
首先判断从地球到月球是否存在一条路线,如果不存在那么无解,否则把每个太空站按照每天拆分成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;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值