• 洛谷P2754 [CTSC1999]家园(最大流)


    传送门

    这题思路太强了……大佬们怎么想到的……我这菜鸡根本想不出来……

    先判断是否能到达,对每一艘飞船能到的地方用并查集合并一下,最后判断一下是否连通

    然后考虑几天怎么判断,我们可以枚举。

    每一个点表示“第几个空间站在第几天”这个状态,那么枚举的答案每加一,就要新建所有空间站的点

    源点向每一个地球连一条容量$inf$的边,每个空间站向下一秒的空间站连容量$inf$的边,表示可以转移到下一秒

    然后飞船在哪一个星球可以判断,那么从这一个时间的出发点向下一个时间的到达点连边,容量为飞船载人数

    只要网络流的流量大于等于$k$就说明可以了

    然后实际上不用每次建图,只要每一次新建边,在残量网络上跑dinic就行了

     1 //minamoto
     2 #include<iostream>
     3 #include<cstdio>
     4 #include<cstring>
     5 #include<queue>
     6 #define inf 0x3f3f3f3f
     7 using namespace std;
     8 #define getc() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++)
     9 char buf[1<<21],*p1=buf,*p2=buf;
    10 inline int read(){
    11     #define num ch-'0'
    12     char ch;bool flag=0;int res;
    13     while(!isdigit(ch=getc()))
    14     (ch=='-')&&(flag=true);
    15     for(res=num;isdigit(ch=getc());res=res*10+num);
    16     (flag)&&(res=-res);
    17     #undef num
    18     return res;
    19 }
    20 const int N=105,M=1000005;
    21 int head[M],Next[M],ver[M],edge[M],dep[M],tot=1;
    22 int fa[N],p[N],g[N][N],num[N];
    23 int n,m,s,t,k,ans,mx;
    24 queue<int> q;
    25 int ff(int x){return fa[x]==x?x:fa[x]=ff(fa[x]);}
    26 void unite(int x,int y){
    27     x=ff(x),y=ff(y);
    28     if(x!=y) fa[x]=y;
    29 }
    30 inline void add(int u,int v,int e){
    31     ver[++tot]=v,Next[tot]=head[u],head[u]=tot,edge[tot]=e;
    32     ver[++tot]=u,Next[tot]=head[v],head[v]=tot,edge[tot]=0;
    33 }
    34 bool bfs(){
    35     while(!q.empty()) q.pop();
    36     for(int i=1;i<=ans*(n+1);++i) dep[i]=-1;
    37     q.push(s),dep[s]=0,dep[t]=-1;
    38     while(!q.empty()){
    39         int u=q.front();q.pop();
    40         for(int i=head[u];i;i=Next[i]){
    41             int v=ver[i];
    42             if(dep[v]<0&&edge[i]){
    43                 dep[v]=dep[u]+1,q.push(v);
    44                 if(v==t) return true;
    45             }
    46         }
    47     }
    48     return false;
    49 }
    50 int dfs(int u,int limit){
    51     if(!limit||u==t) return limit;
    52     int flow=0,f;
    53     for(int i=head[u];i;i=Next[i]){
    54         int v=ver[i];
    55         if(dep[v]==dep[u]+1&&(f=dfs(v,min(limit,edge[i])))){
    56             flow+=f,limit-=f;
    57             edge[i]-=f,edge[i^1]+=f;
    58             if(!limit) break;
    59         }
    60     }
    61     return flow;
    62 }
    63 int main(){
    64     n=read(),m=read(),k=read();
    65     s=0,t=M-2;
    66     for(int i=1;i<=n+2;++i) fa[i]=i;
    67     for(int i=1;i<=m;++i){
    68         p[i]=read(),num[i]=read();
    69         for(int j=0;j<num[i];++j){
    70             g[i][j]=read();
    71             if(g[i][j]==0) g[i][j]=n+1;
    72             if(g[i][j]==-1) g[i][j]=n+2;
    73             if(j!=0) unite(g[i][j],g[i][j-1]);
    74         }
    75     }
    76     if(ff(n+1)!=ff(n+2)) return puts("0"),0;
    77     for(ans=1;;++ans){
    78         add(s,(ans-1)*(n+1)+n+1,inf);
    79         for(int i=1;i<=m;++i){
    80             int x=(ans-1)%num[i],y=ans%num[i];
    81             if(g[i][x]==n+2) x=t;
    82             else x=(ans-1)*(n+1)+g[i][x];
    83             if(g[i][y]==n+2) y=t;
    84             else y=ans*(n+1)+g[i][y];
    85             add(x,y,p[i]);
    86         }
    87         while(bfs()) mx+=dfs(s,inf);
    88         if(mx>=k) return printf("%d
    ",ans),0;
    89         for(int i=1;i<=n+1;++i) add((ans-1)*(n+1)+i,ans*(n+1)+i,inf);
    90     }
    91     return 0;
    92 }
  • 相关阅读:
    ExplorerControls的显示问题
    VS选项中没有C#相关设置选项?
    根据点提取栅格值
    原来我的代码暴露在外面!
    .NET反编译之manager,base.AutoScaleMode修复
    无法加载Dll”ArcGISVersion.dll”:0x8007007E
    ArcEngine9.3迁移至ArcObject10.1
    地籍宗地出图(二)
    地籍宗地出图(一)
    CAD设置图层透明显示
  • 原文地址:https://www.cnblogs.com/bztMinamoto/p/9509268.html
Copyright © 2020-2023  润新知