[CTSC1999]家园

题目链接

由于数据范围较小,可以对每一天进行考虑。

考虑每一天时,对每个太空站都建个新的点代表当前天的太空站,地球和月球也新建。

从源点向地球连容量为inf的边,月球向汇点连容量为inf的边。

然后每个太空站可以从上一天留下\(inf\)个人,所以从前一天连一条容量为\(inf\)的边。

考虑每个太空船,可以从前一天所在的点转移\(h_i\)个人过来,所以连一条容量为\(h_i\)的边。

然后跑最大流,最大流为当前天可送达的人数。

无解可以用并查集判,也可以判如果天数大于500,即一个一个的送,经过所有点都送不完的话,就无解。

常数巨大的代码

#include<bits/stdc++.h>
using namespace std;
 
inline void read(int &x) {
    x=0;int f=1;char ch=getchar();
    for(;!isdigit(ch);ch=getchar()) if(ch=='-') f=-f;
    for(;isdigit(ch);ch=getchar()) x=(x<<1)+(x<<3)+(ch^48);x*=f;
}
 
inline void print(int x) {
    if(x<0) putchar('-'),x=-x;
    if(!x) return ;print(x/10),putchar((x%10)^48);
}
inline void write(int x) {if(!x) putchar('0');else print(x);putchar('\n');}
 
#define maxn 1000050

int n,m,s,t,max_flow,tot=1,k;
int head[maxn],dis[maxn],vis[maxn],a[200][200],num[200],v[200];
struct edge{int to,nxt,w;}e[maxn<<1];

void add(int u,int v,int w) {e[++tot]=(edge){v,head[u],w},head[u]=tot;}
void ins(int u,int v,int w) {add(u,v,w),add(v,u,0);}

int bfs() {
    memset(vis,0,sizeof vis);
    memset(dis,63,sizeof dis);//write(dis[0]);
    queue<int > q;q.push(s);vis[s]=1;dis[s]=0;
    while(!q.empty()) {
        int now=q.front();q.pop();vis[now]=0;
        for(int i=head[now];i;i=e[i].nxt)
            if(e[i].w>0&&dis[e[i].to]>dis[now]+1) {
                dis[e[i].to]=dis[now]+1;
                if(!vis[e[i].to]) vis[e[i].to]=1,q.push(e[i].to);
            }
    }return dis[t]<1e9;
}

int dfs(int u,int flow) {
    if(u==t) return flow;
    for(int i=head[u];i;i=e[i].nxt)
        if(e[i].w>0&&dis[e[i].to]==dis[u]+1) {
            int f=dfs(e[i].to,min(flow,e[i].w));
            if(f>0) {e[i].w-=f,e[i^1].w+=f;return f;}
        }
    dis[u]=-1;return 0;
}

int dinic() {
    max_flow=0;
    while(bfs()) {
        int flow;
        while((flow=dfs(s,1e9))) max_flow+=flow;
    }return max_flow;
}

int main() {
    s=0,t=1e4;
    read(n),read(m),read(k);
    for(int i=1;i<=m;i++) {
        read(v[i]),read(num[i]);
        for(int j=0;j<num[i];j++) read(a[i][j]),a[i][j]+=2;
    }
    int ans=0,res=0;
    while(1) {
        ins(s,ans*(n+2)+2,1e9);
        ins(ans*(n+2)+1,t,1e9);
        if(ans) {
            for(int i=1;i<=n;i++) ins((ans-1)*(n+2)+i+2,ans*(n+2)+i+2,1e9);
            for(int i=1;i<=m;i++) ins((ans-1)*(n+2)+a[i][(ans-1)%num[i]],ans*(n+2)+a[i][ans%num[i]],v[i]);
        }
        res+=dinic();if(res>=k) return write(ans),0;
        ans++;if(ans==500) return puts("0"),0;
    }
    return 0;
}

posted @ 2018-11-07 15:09  Hyscere  阅读(222)  评论(0编辑  收藏  举报