洛谷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 }
原文地址:https://www.cnblogs.com/bztMinamoto/p/9509268.html