2020 icpc 济南

Eter`nal

于 2021-12-02 20:49:37 发布

阅读量271

点赞数

分类专栏： icpc 文章标签：题解

本文链接：https://blog.csdn.net/m0_51488693/article/details/121686146

版权

icpc 专栏收录该内容

7 篇文章 0 订阅

订阅专栏

Matrix Equation 高斯消元

大意：给定 01 矩阵 A,C

求满足 $A\times B=B\odot C$ 01 矩阵 B 的个数。

$\times$ : $D_{i,j}=\sum_{k=1}^n{A_{i,k}B_{k,j}}$

$\odot$ : $D_{i,j}=A_{i,j}B_{i,j}$

思路：通过观察上述式子我们发现 $D_{i,j}$ 叉乘是 k 在变化，列是不变的，对于B的各个列之间是相互独立的。所以我们可以对 B 的每一列单独拉出来考虑。对于 B 的第 j 列记为 b

根据式子我们可以列出：

$A_{1,1}b_1+A_{1,2}b_2+......A_{1,n}bn)\mod2=A_{1,j}b_1$

$A_{2,1}b_1+A_{2,2}b_2+......A_{2,n}bn)\mod2=A_{2,j}b_2$

…

$A_{n,1}b_1+A_{n,2}b_2+......A_{n,n}bn)\mod2=A_{n,j}b_n$

有 mod 2 的存在，所以也就等价于异或线性方程组：

$A_{1,1}b_1\oplus A_{1,2}b_2\oplus......A_{1,n}bn=A_{1,j}b_1$

$A_{2,1}b_1\oplus A_{2,2}b_2\oplus......A_{2,n}bn=A_{2,j}b_2$

…

$A_{n,1}b_1\oplus A_{n,2}b_2\oplus......A_{n,n}bn=A_{n,j}b_n$

继续化简：

$(A_{1,1}\oplus C_{1,j})b_1\oplus A_{1,2}b_2\oplus......A_{1,n}bn=0$

$A_{2,1}b_1\oplus (A_{2,2}\oplus C_{2,j})b_2\oplus......A_{2,n}bn=0$

…

$A_{n,1}b_1\oplus A_{n,2}b_2\oplus......(A_{n,n}\oplus C_{n,j})b_n=0$

然后用高斯消元求自由元的个数就行了。自由元的个数为 cnt ，该列对最终结果的贡献就是 $2^{cnt}$

乘法原理把每一列乘起来就行了。

代码如下：

#include <bits/stdc++.h>
#define int long long 
#define rep(i,bbb,eee) for(int i=bbb;i<=eee;i++)
#define frep(i,bbb,eee) for(int i=bbb;i>=eee;i--)
#define mem(a,b) memset(a,b,sizeof(a))
#define cout(xxx) cout<<fixed<<setprecision(xxx)
#define inf 0x3f3f3f3f
#define pb push_back
#define AC signed
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
const int N=1010,M=998244353;

int qmi(int a,int b)
{
	int res=1;
	while(b)
	{
		if(b&1)res=res*a%M;
		a=a*a%M;
		b>>=1;
	}
	return res%M;
}
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
inline void print(int x)
{   
   if(x<0){putchar('-');x=-x;}
   if(x>9) print(x/10);
   putchar(x%10+'0');
}
//____________________________________________//
int n;
int A[N][N],C[N][N],a[N][N];
int gauss()
{
	int r,c;
	for(r=0,c=0;c<n;c++)
	{
		int t=-1;
		for(int i=r;i<n;i++)
			if(a[i][c])t=i;

		if(t==-1)continue;

		for(int i=c;i<=n;i++)swap(a[t][i],a[r][i]);

		for(int i=r+1;i<n;i++)
			if(a[i][c])
				for(int j=n;j>=c;j--)
					a[i][j]^=a[r][j];

		r++;
	}
	return n-r;
}
void solve()
{
	cin>>n;
	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			cin>>A[i][j];

	for(int i=0;i<n;i++)
		for(int j=0;j<n;j++)
			cin>>C[i][j];

	int ans=1;
	for(int j=0;j<n;j++)
	{
		for(int i=0;i<n;i++)
			for(int k=0;k<n;k++)
				a[i][k]=A[i][k];

		for(int i=0;i<n;i++)a[i][i]^=C[i][j];
		int t=gauss();
		ans=ans*qmi(2,t)%M;
	}
	cout<<ans<<"\n";
}
AC main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int _=1;
	//cin>>_;
	while(_--)solve();
	return 0;
}

Tree Constructer 构造+二分图染色

大意：若u的点权为 x ，v 的点权为 y 。当且仅当 $x|y=2^{60}-1$ 时u,v 之间有条边。点的个数<=100，现给出一个棵树，要求根据树中边的关系构造出一组合法点权。

思路：构造题重点就是就是想到一种简单的构造方式。根据题目，我们要找到是 a-b-c

$a|b=b|c=2^{60}-1$ 即 a、c 是一列，b是另一列，类似于染色，黑白色交替。

如果我们令点数目较少的一类为白色点的话，白色点的数目最大不超过50。要保证白色点的之间或运算不为 $2^{60}-1$ , 我们可以令白色点最高0，且第 i 号白色点的第 i 位为0，这样确保白色点之间是合法的，对于黑色点，最高位是1，周围所有的白色点的编号与之对应的位是1。这样黑白之间也是合法的。因为白色点不超过50个，除了最高位的1，其他的是1的位最高是50，那么更高的位必然全部都是0，黑色点之间也是合法的。

代码如下：

#include <bits/stdc++.h>
#define int long long 
#define rep(i,bbb,eee) for(int i=bbb;i<=eee;i++)
#define frep(i,bbb,eee) for(int i=bbb;i>=eee;i--)
#define mem(a,b) memset(a,b,sizeof(a))
#define cout(xxx) cout<<fixed<<setprecision(xxx)
#define inf 0x3f3f3f3f
#define pb push_back
#define AC signed
#define fi first
#define se second
using namespace std;
typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
const int N=110,M=998244353;

int qmi(int a,int b)
{
	int res=1;
	while(b)
	{
		if(b&1)res=res*a%M;
		a=a*a%M;
		b>>=1;
	}
	return res%M;
}
inline int read()
{
    int x=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)&&ch!='-')ch=getchar();
    if(ch=='-')f=-1,ch=getchar();
    while(isdigit(ch))x=x*10+ch-'0',ch=getchar();
    return f*x;
}
inline void print(int x)
{   
   if(x<0){putchar('-');x=-x;}
   if(x>9) print(x/10);
   putchar(x%10+'0');
}
//____________________________________________//
int n,ans[N],id[N];
vector<int> g[N],w,b;
void dfs(int u,int fa,int c)
{
	if(c)w.pb(u);
	else b.pb(u);
	for(auto v:g[u])
	{
		if(v==fa)continue;
		dfs(v,u,c^1);
	}
}
void solve()
{
	cin>>n;
	rep(i,1,n)id[i]=-1;
	rep(i,1,n-1)
	{
		int x,y;
		cin>>x>>y;
		g[x].pb(y);
		g[y].pb(x);
	}
	dfs(1,0,0);
	if(w.size()>b.size())swap(w,b);
	int t=(1LL<<59)-1;
	for(int i=0;i<w.size();i++)
	{
		int u=w[i];
		id[u]=i;		
		ans[u]=t^(1LL<<i);
	}
	t=1LL<<59;
	for(int i=0;i<b.size();i++)
	{
		int u=b[i];
        ans[u]=t;
		for(auto v:g[u])
		{
        
			if(id[v]!=-1)
			{
				ans[u]|=(1LL<<id[v]);
			}
		}
	}
	rep(i,1,n)cout<<ans[i]<<" ";
}
AC main()
{
	ios::sync_with_stdio(false);cin.tie(0);
	int _=1;
	//cin>>_;
	while(_--)solve();
	return 0;
}

Eter`nal

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
0
评论
2020 icpc 济南

Matrix Equation 高斯消元大意：给定 01 矩阵 A,C求满足 A×B=B⊙CA\times B=B\odot CA×B=B⊙C 01 矩阵 B 的个数。×\times× : Di,j=∑k=1nAi,kBk,jD_{i,j}=\sum_{k=1}^n{A_{i,k}B_{k,j}}Di,j=∑k=1nAi,kBk,j⊙\odot⊙: Di,j=Ai,jBi,jD_{i,j}=A_{i,j}B_{i,j}Di,j=Ai,jBi,j思路：通过观察上述式子我们发现Di.
复制链接

扫一扫