bzoj3442 学习小组

题目链接：bzoj3442
题目大意：
共有n个学生，m个学习小组，每个学生只愿意参加其中的一些学习小组，但校领导规定一个学生最多参加k个学习小组。每个学生参加学习小组都要交一定的手续费，不同的学习小组有不同的手续费。然而若有a个学生参加第i个学习小组，那么就要给这个学习小组组织者奖励Ci*a^2元。在参与学生（而不是每个学习小组的人数总和）尽量多的情况下，求财务处最少要支出多少钱（若为负数，则输出负数）（支出=总奖励费-总手续费）。

题解：
最大流
如果忽略参与学生尽量多这个限制的话，我还是会做的QwQ

• $S\rightarrow i$，从源点到学生连一条流量为$k$，费用为0的边
• 根据兴趣连学生到兴趣小组的边，流量为1费用为0
• 兴趣小组到T连边，看到这种套路Ci*a^2就是要拆边，$a^2=1+3+5+\dots$。噢要记得给拆出来的每条边都减上一个学生要交的手续费。

但是还要求参与学生尽量多、、
其实，它翻译一下就是：如果学生有兴趣的话，至少要参加一个小组，其余的不做硬性规定。那么构图只要对每个学生再加一条$i\rightarrow T$流量为k-1费用为0的边就可以了。

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<iostream>
#include<algorithm>
using namespace std;
#define maxn 210
#define maxm 150010

const int inf=1e9;
struct node
{
    int y,c,f,next,ot;
}a[maxm*2];int len,first[maxn];
void ins(int x,int y,int c,int f)
{
    len++;int n1=len;a[len].y=y;a[len].c=c;
    a[len].f=f;a[len].next=first[x];first[x]=len;
    len++;int n2=len;a[len].y=x;a[len].c=-c;
    a[len].f=0;a[len].next=first[y];first[y]=len;
    a[n1].ot=n2;a[n2].ot=n1;
}
queue<int> q;bool vis[maxn];
int S,T,pre[maxn],d[maxn],l[maxn],flow[maxn];
int mymin(int x,int y){return (x<y)?x:y;}
int bfs()
{
    for (int i=1;i<=T;i++) {pre[i]=-1;d[i]=inf;}
    vis[S]=true;d[S]=0;flow[S]=inf;q.push(S);
    while (!q.empty())
    {
        int x=q.front();q.pop();
        for (int k=first[x];k!=-1;k=a[k].next)
        {
            int y=a[k].y;
            if (a[k].f<=0) continue;
            if (pre[y]==-1 || d[y]>d[x]+a[k].c)
            {
                d[y]=d[x]+a[k].c;
                pre[y]=x;l[y]=k;
                flow[y]=mymin(flow[x],a[k].f);
                if (!vis[y])
                {
                    vis[y]=true;
                    q.push(y);
                }
            }
        }vis[x]=false;
    }
    if (pre[T]==-1) return 0;
    return flow[T];
}
int MCMF()
{
    int delta,ret=0;
    while (delta=bfs())
    {
        ret+=delta*d[T];
        int x=T;
        while (x!=S)
        {
            a[l[x]].f-=delta;
            a[a[l[x]].ot].f+=delta;
            x=pre[x];
        }
    }
    return ret;
}
char s[110];int c[110];
int main()
{
    //freopen("a.in","r",stdin);
    //freopen("a.out","w",stdout);
    int n,m,k,x,i,j;
    scanf("%d%d%d",&n,&m,&k);
    len=0;S=n+m+1;T=S+1;
    memset(first,-1,sizeof(first));
    for (i=1;i<=n;i++)
    {
        ins(S,i,0,k);
        ins(i,T,0,k-1);
    }
    for (i=1;i<=m;i++) scanf("%d",&c[i]);
    for (i=1;i<=m;i++) 
    {
        scanf("%d",&x);
        for (j=1;j<=100;j++)
         ins(i+n,T,c[i]*(2*j-1)-x,1);
    }
    for (i=1;i<=n;i++)
    {
        scanf("%s",s+1);
        for (j=1;j<=m;j++)
         if (s[j]=='1') ins(i,j+n,0,1);
    }
    printf("%d\n",MCMF());
    return 0;
}