zoukankan      html  css  js  c++  java
  • 4.26 ABC F I hate Matrix Construction 二进制拆位 构造 最大匹配

    LINK:I hate Matrix Construction

    心情如题目名称。

    主要说明一下构造的正确性。

    准确来说这道题困扰我很久。

    容易发现可以拆位构造。

    这样题目中的条件也比较容易使用。

    最后等价于每一行每一列有一个 当前行/列有一个1或者0的限制。

    考虑直接进行构造。

    容易发现这类似于最大匹配 尝试利用最大匹配来做 不过这样的话时间复杂度会爆掉且不好写。

    可以观察到一个性质 一个点只会为一行或者一列提供贡献 如果行列的需求一样那么这个点直接放即最优。

    最后问题变成了 有一些没有放值得位置 要求合理的放值满足一些 行或列的限制 且这些位置最多只会满足行列中的一种的约束。

    关于剩下的条件 是否存在增广路的问题:如果当且列放了值影响到行的最优决策了 首先点不是共用的 且这个行会影响到其他列。

    如果存在一条增广路 但是考虑到 此时的列 和其他的列是相同的 所以当且列也直接存在增广路 所以影响关系是不存在的。

    综上 直接进行构造 因为不存在类似于最大匹配的增广路存在。

    一种比较简单的构造方法:

    先让剩下的全部变成0 对于1找0的存在 看能否更替即可。

    最后注意再check一下 可能变得不合法了。

    写法是借鉴别人的 我也同时意识到 随便写代码复杂度很高。

    const ll MAXN=510;
    int n,cc;
    int s[MAXN],t[MAXN];
    ull v[MAXN],u[MAXN];
    int b[MAXN][MAXN];	
    int r[MAXN][2],c[MAXN][2];
    ull a[MAXN][MAXN];
    inline int solve()
    {
    	rep(0,63,k)
    	{
    		memset(b,-1,sizeof(b));
    		memset(c,0,sizeof(c));
    		memset(r,0,sizeof(r));
    		rep(1,n,i)
    		{
    			int x=u[i]&1;
    			if(s[i]==x)continue;
    			rep(1,n,j)if(b[i][j]!=-1&&b[i][j]!=x)return 0;
    			else b[i][j]=x;
    		}
    		rep(1,n,i)
    		{
    			int x=v[i]&1;
    			if(t[i]==x)continue;
    			rep(1,n,j)if(b[j][i]!=-1&&b[j][i]!=x)return 0;
    			else b[j][i]=x;
    		}
    		rep(1,n,i)rep(1,n,j)
    		{
    			if((u[i]&1)==(v[j]&1))b[i][j]=(u[i]&1);
    			if(b[i][j]==-1)b[i][j]=0;
    			++r[i][b[i][j]];
    			++c[j][b[i][j]];
    		}
    		rep(1,n,i)
    			if(s[i]&&(u[i]&1))
    			{
    				if(r[i][1])continue;
    				rep(1,n,j)
    				{
    					if(!t[j]&&!(v[j]&1)&&c[j][0]>1)
    					{
    						b[i][j]=1;
    						--c[j][0];--r[i][0];
    						++c[j][1];++r[i][1];
    						break;
    					}
    				}
    			}
    		rep(1,n,i)
    		{
    			if(t[i]&&(v[i]&1))
    			{
    				if(c[i][1])continue;
    				rep(1,n,j)
    				{
    					if(!s[j]&&!(u[j]&1)&&r[j][0]>1)
    					{
    						b[j][i]=1;
    						--r[j][0];--c[i][0];
    						++r[j][1];++c[i][1];
    						break;
    					}
    				}
    			}
    		}
    		rep(1,n,i)
    		{
    			if((u[i]&1)&&!s[i]&&r[i][1]!=n)return 0;
    			if(!(u[i]&1)&&s[i]&&r[i][0]!=n)return 0;
    			if(!(v[i]&1)&&t[i]&&c[i][0]!=n)return 0;
    			if((v[i]&1)&&!t[i]&&c[i][1]!=n)return 0;
    			if((u[i]&1)==s[i]&&!s[i]&&!r[i][0])return 0;
    			if((u[i]&1)==s[i]&&s[i]&&!r[i][1])return 0;
    			if((v[i]&1)==t[i]&&t[i]&&!c[i][1])return 0;
    			if((v[i]&1)==t[i]&&!t[i]&&!c[i][0])return 0;
    		}
    		rep(1,n,i)rep(1,n,j)if(b[i][j])a[i][j]=a[i][j]|(1ull<<k);
    		rep(1,n,i)u[i]>>=1,v[i]>>=1;
    	}
    	return 1;
    }
    signed main()
    {
    	freopen("1.in","r",stdin);
    	gt(n);
    	rep(1,n,i)gt(s[i]);
    	rep(1,n,i)gt(t[i]);
    	rep(1,n,i)scanf("%llu",&u[i]);
    	rep(1,n,i)scanf("%llu",&v[i]);
    	if(!solve())puts("-1");
    	else rep(1,n,i)rep(1,n,j)printf("%llu%c",a[i][j],j!=n?' ':'
    ');
    	return 0;
    }
    
  • 相关阅读:
    算法题-数组算法题
    Linux-shell脚本的调试和追踪
    Linux-循环loop
    Linux-条件判断式
    Linux-第一行#!/bin/bash的含义
    Linux-排序命令:sort、wc、uniq
    Linux-选取命令:cut grep
    Linux-shell变量
    Linux-重定向、追加、tee
    Linux三剑客-sed编辑文本
  • 原文地址:https://www.cnblogs.com/chdy/p/12787454.html
Copyright © 2011-2022 走看看