zoukankan      html  css  js  c++  java
  • 【[网络流24题]方格取数问题】

    发现这个图还挺不好连得

    因为并不是非常经典的最大权闭合子图了

    试图拆点也发现没有什么效果

    但是转念一想发现最大权闭合子图之所以能跑最小割就是因为我们强行分组了

    这里是不是也可以强行分组啊

    于是我们可以强行分组,把所有点分成黑点和白点,黑点向(T)连边,(S)向白点连边,边权是这个点的权值

    之后白点向他能控制的黑点连边就行了

    代码

    #include<algorithm>
    #include<iostream>
    #include<cstring>
    #include<cstdio>
    #include<queue>
    #define maxn 100005
    #define re register
    #define LL long long
    #define inf 999999999
    #define max(a,b) ((a)>(b)?(a):(b))
    #define min(a,b) ((a)<(b)?(a):(b))
    const int dx[]={0,0,1,-1};
    const int dy[]={1,-1,0,0};
    inline int read()
    {
    	char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
    	while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();
    	return x;
    }
    struct E{int v,nxt,w,f;}e[maxn<<3];
    int d[maxn],head[maxn],cur[maxn];
    int a[101][101],b[101][101],c[101][101];
    int n,m,num=1,S,T,tot,ans;
    inline void add(int x,int y,int z) {e[++num].v=y;e[num].w=z;e[num].nxt=head[x];head[x]=num;}
    inline void C(int x,int y,int z) {add(x,y,z),add(y,x,0);}
    inline int BFS()
    {
    	std::queue<int> q;
    	for(re int i=S;i<=T;i++) cur[i]=head[i],d[i]=0;
    	d[S]=1,q.push(S);
    	while(!q.empty())
    	{
    		int k=q.front();q.pop();
    		for(re int i=head[k];i;i=e[i].nxt)
    		if(!d[e[i].v]&&e[i].w>e[i].f) d[e[i].v]=d[k]+1,q.push(e[i].v);
    	}
    	return d[T];
    }
    int dfs(int x,int now)
    {
    	if(x==T||!now) return now;
    	int flow=0,ff;
    	for(re int& i=cur[x];i;i=e[i].nxt)
    	if(d[e[i].v]==d[x]+1)
    	{
    		ff=dfs(e[i].v,min(now,e[i].w-e[i].f));
    		if(ff<=0) continue;
    		now-=ff,flow+=ff;
    		e[i].f+=ff,e[i^1].f-=ff;
    		if(!now) break;
    	}
    	return flow;
    }
    int main()
    {
    	n=read(),m=read();
    	for(re int i=1;i<=n;i++)
    		for(re int j=1;j<=m;j++)
    			a[i][j]=read(),b[i][j]=++tot,ans+=a[i][j];
    	for(re int i=1;i<=n;i++)
    		for(re int j=1;j<=m;j++) c[i][j]=(i+j)%2;
    	T=++tot;
    	for(re int i=1;i<=n;i++)
    		for(re int j=1;j<=m;j++)
    			if(c[i][j]) C(S,b[i][j],a[i][j]);else C(b[i][j],T,a[i][j]);
    	for(re int i=1;i<=n;i++)
    		for(re int j=1;j<=m;j++)
    			for(re int k=0;c[i][j]&&k<4;k++)
    			{
    				int xx=i+dx[k],yy=j+dy[k];
    				if(xx>=1&&yy>=1&&xx<=n&&yy<=m) if(!c[xx][yy]) C(b[i][j],b[xx][yy],inf);
    			}
    	while(BFS()) ans-=dfs(S,inf);
    	printf("%d
    ",ans);
    	return 0;
    }
    
  • 相关阅读:
    hexo及next主题修改
    LeetCode#476 Number Complement
    Html学习笔记(二) 简单标签
    Haproxy的应用
    STM32 一个初始化EXTI的例子
    sql语句优化原则
    在Docker中运行gocd
    Gnome Ubuntu16安装Nvidia显卡396驱动,CUDA9.2以及cudnn9.2
    吴裕雄--天生自然数据结构:十大经典排序算法——希尔排序
    吴裕雄--天生自然数据结构:十大经典排序算法——插入排序
  • 原文地址:https://www.cnblogs.com/asuldb/p/10290698.html
Copyright © 2011-2022 走看看