传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=3661
思路:
网络流,但这不是标算,bzoj为了大概是为了让网络流过才把时限开到100s
建两个源点,S向S‘连容量为K的边,这个是为了限制总流量的
每天的每只兔子拆成两个点,之间连1的边进行点容量限制,表示这只兔子今天只算一只(不限制后面会出问题)
每只兔子向下一天的自己连1的边,表示不用换人
每天设立两个中转点,之间连L的边,拆成两个也是为了限制流量
每只兔子向中转点连边,中转点再向下一天的每只兔子连边,表示可以换L只兔子
点数1280000+,边数不超过5000000
然后就没有然后了
正解与污误解的分界线
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
const int N=820,maxn=1500010,maxm=5000010,inf=1061109567;
using namespace std;
int n,m,K,L,map[N][N],rab[N][N][2],cnt,lim[N][2],Edge[N][N];char ch[N];
struct Flow{
int pre[maxm],now[maxn],son[maxm],val[maxm],tot,dis[maxn],S,SS,T,q[maxm+10],head,tail;
void add(int a,int b,int c){pre[++tot]=now[a],now[a]=tot,son[tot]=b,val[tot]=c;}
void ins(int a,int b,int c){add(a,b,c),add(b,a,0);}
void init(){memset(now,0,sizeof(now)),tot=1,S=maxn-3,SS=maxn-2,T=maxn-1;}
bool bfs(){
memset(dis,-1,sizeof(dis));
q[tail=1]=S,dis[S]=head=0;
while (head!=tail){
if (++head>maxn) head=1;
int x=q[head];
for (int y=now[x];y;y=pre[y])
if (dis[son[y]]==-1&&val[y]){
if (++tail>maxn) tail=1;
dis[son[y]]=dis[x]+1,q[tail]=son[y];
}
}
return dis[T]>0;
}
int find(int x,int low){
if (x==T) return low;
int y,res=0;
for (y=now[x];y;y=pre[y]){
if (dis[son[y]]!=dis[x]+1||!val[y]) continue;
int tmp=find(son[y],min(low,val[y]));
res+=tmp,low-=tmp,val[y]-=tmp,val[y^1]+=tmp;
if (!low) break;
}
if (!y) dis[x]=-1;
return res;
}
void build(){
ins(S,SS,K);
for (int i=1;i<=m;i++)
for (int j=1;j<=n;j++){
ins(rab[i][j][0],rab[i][j][1],1);//点容量限制
Edge[i][j]=tot-1;
}
for (int i=1;i<=n;i++) if (map[1][i]) ins(SS,rab[1][i][0],1);
for (int i=2;i<=m;i++)
for (int j=1;j<=n;j++)
if (map[i][j]) ins(rab[i-1][j][1],rab[i][j][0],1);//向下一天兔子连边
for (int i=1;i<m;i++)
for (int j=1;j<=n;j++)
if (map[i][j]) ins(rab[i][j][1],lim[i][0],1);
for (int i=1;i<m;i++) ins(lim[i][0],lim[i][1],L);
for (int i=1;i<m;i++)
for (int j=1;j<=n;j++)
if (map[i+1][j]) ins(lim[i][1],rab[i+1][j][0],1);
for (int i=1;i<=n;i++) ins(rab[m][i][1],T,1);
}
void work(){
int res=0;
while (bfs()) res+=find(S,inf);
if (res!=K){puts("-1");return;}
for (int i=1;i<=m;i++,puts(""))
for (int j=1;j<=n;j++)
if (!val[Edge[i][j]]) printf("%d ",j);
}
}F;
int main(){
scanf("%d%d%d%d",&n,&m,&K,&L);F.init();
for (int i=1;i<=n;i++){
scanf("%s",ch+1);
for (int j=1;j<=m;j++) map[j][i]=ch[j]-'0',rab[j][i][0]=++cnt,rab[j][i][1]=++cnt;
}
for (int i=1;i<=m;i++) lim[i][0]=++cnt,lim[i][1]=++cnt;
F.build(),F.work();
return 0;
}