[bzoj4823][Cqoi2017]老C的方块

来自FallDream的博客，未经允许，请勿转载，谢谢。

---

挺有意思的一道题....

看完题面比较明确是最小割，考虑怎么建图 想了比较久

突破口应该是题目中那张奇怪的图

观察这个奇怪的图和方块，很容易发现每个图案，其实都是每个分割线周围各一个1*2的块拼在一起的

假如中间有个分割线，那么就是这样，两边分别三选一

然后可以依此根据块的位置不同建出图，我大概画了一部分，大家可以参考

这张图中行从下到上，列从左到右

然后注意把点拆成入点和出点，中间连费用的边 最小割即可。

可以哈希来查点

#include<iostream>
#include<cstdio>
#include<map>
#include<cstring>
#define S 0
#define MN 100000
#define INF 2000000000
#define ll long long
using namespace std;
inline int read()
{
    int x =  , f = ; char ch = getchar();
    while(ch < '' || ch > ''){ if(ch == '-') f = -;  ch = getchar();}
    while(ch >= '' && ch <= ''){x = x *  + ch - '';ch = getchar();}
    return x * f;
}

map<ll,int> mp;
struct edge{int to,next,w;}e[];
int n,m,c,T,head[MN*+],C[MN*+],q[MN*+],top=,d[MN*+],cnt=,ans=;
struct Block{int x,y,w;}b[MN+];
inline void ins(int f,int t,int w)
{
    if(!t) return;
    e[++cnt]=(edge){t,head[f],w};head[f]=cnt;
    e[++cnt]=(edge){f,head[t],};head[t]=cnt;
}
inline int num(int x,int y)
{
    if(x>n||y>m||x<||y<) return ;
    int res=mp[1LL*x*MN+y];
    return res?res+c:;
}

int dfs(int x,int f)
{
    if(x==T) return f;
    int used=;
    for(int&i=C[x];i;i=e[i].next)
        if(e[i].w&&d[e[i].to]==d[x]+)
        {
            int w=dfs(e[i].to,min(f-used,e[i].w));
            used+=w;e[i].w-=w;e[i^].w+=w;
            if(used==f) return f;
        }
    return d[x]=-,used;
}

bool bfs()
{
    memset(d,,sizeof(int)*(T+));int i,j;
    for(d[q[top=i=]=S]=;i<=top;++i)
        for(int j=C[q[i]]=head[q[i]];j;j=e[j].next)
            if(e[j].w&&!d[e[j].to])
                d[q[++top]=e[j].to]=d[q[i]]+;
    return d[T];
}

int main()
{
    m=read();n=read();c=read();T=c*+;
    for(int i=;i<=c;++i) b[i].y=read(),b[i].x=read(),b[i].w=read(),mp[1LL*b[i].x*MN+b[i].y]=i;
    for(int i=;i<=c;++i)
    {
        ins(i+c,i,b[i].w);
        if(b[i].x&)
        {
            int cas=b[i].y%;
            if(cas==)
            {
                ins(S,i+c,INF);
                ins(i,num(b[i].x+,b[i].y),INF);
                ins(i,num(b[i].x-,b[i].y),INF);
                ins(i,num(b[i].x,b[i].y+),INF);
            }
            if(cas==) ins(i,num(b[i].x,b[i].y+),INF);
            if(cas==)
            {
                ins(i,num(b[i].x+,b[i].y),INF);
                ins(i,num(b[i].x-,b[i].y),INF);
                ins(i,num(b[i].x,b[i].y+),INF);
            }
            if(cas==) ins(i,T,INF);
        }
        else
        {
            int cas=b[i].y%;
            if(cas==) ins(i,num(b[i].x,b[i].y-),INF);
            if(cas==)
            {
                ins(S,i+c,INF);
                ins(i,num(b[i].x,b[i].y-),INF);
                ins(i,num(b[i].x+,b[i].y),INF);
                ins(i,num(b[i].x-,b[i].y),INF);
            }
            if(cas==) ins(i,T,INF);
            if(cas==)
            {
                ins(i,num(b[i].x,b[i].y-),INF);
                ins(i,num(b[i].x-,b[i].y),INF);
                ins(i,num(b[i].x+,b[i].y),INF);
            }
        }
    }
    while(bfs()) ans+=dfs(S,INF);
    printf("%d\n",ans);
    return ;
}