1305: [CQOI2009]dance跳舞

题目链接

题目大意：给定n个男生和n个女生，一些互相喜欢而一些不，举行几次舞会，每次舞会要配成n对，不能有相同的组合出现，每个人只能与不喜欢的人跳k次舞，求最多举行几次舞会

题解：将一个人x拆成两个点，分别为lx和dx，表示喜欢/不喜欢

对男孩

$连(lx,dx,k),表示消耗了一次跟不喜欢的人的跳舞机会$
$连(st,lx)，容量待定$

对女孩

$连(dx,lx,k),同上$
$连(lx,ed)，容量待定$

对于两个人x,y

若x喜欢y，$连(lx,ly,1)$
否则，$连(dx,dy,1)$

当然可以二分，但是这里用了动态加点

当答案增大1，就给连接源/汇的边扩容1，直至无法满流，输出ans-1即可

我的收获：网络流枚举答案强啊

#include <iostream>
#include <cstdio>
#include <cstring>
#include <queue>
#include <algorithm>
using namespace std;

const int M=234577;
#define T(x) (x)*n*m
#define INF 0x7fffffff
#define pii pair<int,int>
const int dx[]={0,1,-1,0,0};
const int dy[]={0,0,0,1,-1};

int n,m,t,now,tim,cnt,tot,st,ed,door;
int id[25][25],dis[25][25];
int num[M],d[M],head[M],last[M];
bool Exit;
char s[25][25];

queue<pii> q;

struct edge{int to,c,nex;}e[M*10];

void add(int u,int v,int w){e[t]=(edge){v,w,head[u]};last[u]=head[u]=t++;}

void insert(int x,int y,int z){add(x,y,z),add(y,x,0);};

int dfs(int x,int in)
{
    if(x==ed) return in;
    int ans=0,f;
    for(int i=last[x];i!=-1;last[x]=i=e[i].nex)
    {
        int v=e[i].to;
        if(e[i].c&&d[v]==d[x]-1){
            f=dfs(v,min(in-ans,e[i].c));
            ans+=f;e[i].c-=f,e[i^1].c+=f;
            if(Exit||ans==in) return ans;
        }
    }
    if(--num[d[x]]==0) Exit=1;
    d[x]++,num[d[x]]++,last[x]=head[x];
    return ans;
}

int isap(int tim)
{
    memset(d,0,sizeof(d));
    Exit=0;
    int flow=0;
    while(!Exit) flow+=dfs(st,INF);
    return flow;
}

void build()
{
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(s[i][j]=='D')
            {
                if(tim) insert(id[i][j]+T(tim-1),id[i][j]+T(tim),INF);
                insert(id[i][j]+T(tim),ed,1);
            }
}

void getdis(int x,int y)
{
    memset(dis,0,sizeof(dis));
    q.push(make_pair(x,y));
    while(!q.empty()){
        int nx=q.front().first,ny=q.front().second;q.pop();
        for(int k=1;k<=4;k++){
            int fx=nx+dx[k],fy=ny+dy[k];
            if(dis[fx][fy]||fx<1||fx>n||fy<1||fy>m||s[fx][fy]!='.') continue;
            dis[fx][fy]=dis[nx][ny]+1;
            insert(id[fx][fy],id[x][y]+T(dis[fx][fy]-1),1);
            q.push(make_pair(fx,fy));
        }
    } 
}

void work()
{
    while(now<tot)
    {
        if(tim>=n*m){puts("impossible");return ;}
        build();
        tim++;now+=isap(tim);
    }
    cout<<tim<<endl;
}

void init()
{
    memset(head,-1,sizeof(head));
    memset(last,-1,sizeof(last));
    cin>>n>>m;
    st=0,ed=T(n*m)+1;
    for(int i=1;i<=n;i++) scanf("%s",s[i]+1); 
    for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) id[i][j]=++cnt;//！！！！ 
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++){
            if(s[i][j]=='.') tot++,insert(st,id[i][j],1);
            if(s[i][j]=='D') getdis(i,j);
        }
}

int main()
{
    init();
    work();
    return 0;
}