先不考虑他让参与人数最多,我们可以S连每个人流量k费用0,每个人连能去的小组,流量1费用手续费,小组连T,拆边建费用
考虑让参与人数最多,每个人连T流量k-1费用为0,意义是在他参加了社团的情况下他可以把加入社团的机会让给别的没有社团的人
#include<iostream>
#include<cstdlib>
#include<cstdio>
#include<cstring>
#include<cmath>
#include<ctime>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<stack>
#include<queue>
#include<map>
#include<set>
#include<bitset>
using namespace std;
#define MAXN 1010
#define MAXM 1000010
#define ll long long
#define INF 1000000000
#define MOD 1000000007
#define eps 1e-8
struct vec{
int to;
int fro;
int v;
int f;
};
vec mp[MAXM];
int tai[MAXN],cnt=1;
int q[MAXN],hd,tl;
bool iq[MAXN];
int dis[MAXN],rd[MAXN];
int s,t;
int answ,ansf;
inline void be(int x,int y,int f,int v){
mp[++cnt].to=y;
mp[cnt].fro=tai[x];
tai[x]=cnt;
mp[cnt].f=f;
mp[cnt].v=v;
}
inline void bse(int x,int y,int f,int v){
be(x,y,f,v);
be(y,x,0,-v);
}
bool spfa(){
int i,x,y;
memset(dis,0x3f,sizeof(dis));
hd=tl=0;
dis[s]=0;
iq[q[tl++]=s]=1;
while(hd!=tl){
iq[x=q[hd++]]=0;
hd%=MAXN;
for(i=tai[x];i;i=mp[i].fro){
y=mp[i].to;
if(dis[y]>dis[x]+mp[i].v&&mp[i].f){
dis[y]=dis[x]+mp[i].v;
rd[y]=i;
if(!iq[y]){
iq[q[tl++]=y]=1;
tl%=MAXN;
}
}
}
}
if(dis[t]>INF){
return 0;
}
int flow=INF;
for(i=t;i!=s;i=mp[rd[i]^1].to){
flow=min(flow,mp[rd[i]].f);
}
ansf+=flow;
for(i=t;i!=s;i=mp[rd[i]^1].to){
answ+=flow*mp[rd[i]].v;
mp[rd[i]].f-=flow;
mp[rd[i]^1].f+=flow;
}
return 1;
}
int n,m,k;
int c[MAXN],f[MAXN];
int mx[MAXN];
char S[MAXN];
int main(){
int i,j;
scanf("%d%d%d",&n,&m,&k);
s=n+m+1;
t=s+1;
for(i=1;i<=m;i++){
scanf("%d",&c[i]);
}
for(i=1;i<=m;i++){
scanf("%d",&f[i]);
}
for(i=1;i<=n;i++){
bse(s,i,k,0);
bse(i,t,k-1,0);
scanf("%s",S+1);
for(j=1;j<=m;j++){
if(S[j]=='1'){
bse(i,n+j,1,-f[j]);
mx[j]++;
}
}
}
for(i=1;i<=m;i++){
for(j=1;j<=mx[i];j++){
bse(n+i,t,1,(2*j-1)*c[i]);
}
}
while(spfa()){
}
printf("%d\n",answ);
return 0;
}
/*
3 3 1
1 2 3
3 2 1
111
111
111
*/