ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Ybtoj #891. 「高斯消元」生日礼物

2022-02-06 21:03:11  阅读:233  来源: 互联网

标签:891 int Ybtoj long Ny Fl RI 高斯消 define


题面传送门
如果设\(f_{i,j}\)为当前位置要不要取反,那么显然可以列出\(nm\)个方程。
直接解方程是\(O(\frac{n^3m^3}{w})\)的显然过不去。
考虑每个方程的项数很少,所以可以主元法。具体的,令前两行和第一列为变量,则每个位置\((i,j)\)都可以用\((i-2,j-1)\)位置的等式表达出来,这样每个位置都可以用这\(2m+n-2\)个数表示出来。
然后还有最后两行以及最后一列的方程没用。直接解看自由元的数量即可。答案即为\(2\)的幂次。时间复杂度\(O(\frac{(n+m)^3}{w})\)
code:

#include<bits/stdc++.h>
#define I inline
#define max(a,b) ((a)>(b)?(a):(b))
#define min(a,b) ((a)<(b)?(a):(b))
#define abs(x) ((x)>0?(x):-(x))
#define re register
#define RI re int
#define ll long long
#define db double
#define lb long db
#define N 600
#define M N*N+5
#define mod 123456789
#define Mod (mod-1)
#define eps (1e-9)
#define U unsigned int
#define it iterator
#define Gc() getchar() 
#define Me(x,y) memset(x,y,sizeof(x))
#define d(x,y) (m*(x-1)+(y))
#define R(n) (rand()*rand()%(n)+1)
#define Pc(x) putchar(x)
#define LB lower_bound
#define UB upper_bound  
using namespace std;
int n,m,k,cnt,Nx,Ny,Fl[N*3+5];bitset<N*3+5>F[N*N+5],Q[N*3+5],C;ll Ans=1;
int xp[8]={1,1,-1,-1,2,2,-2,-2};int yp[8]={2,-2,2,-2,1,-1,1,-1};
I void GA(int x,int y){C=F[d(x,y)];for(RI i=0;i<8;i++) {Nx=x+xp[i];Ny=y+yp[i];if(Nx<1||Ny<1|Nx>n||Ny>m) continue;C^=F[d(Nx,Ny)];}}
I void Ins(){RI i;for(i=k;i;i--){if(!C[i]) continue;if(!Fl[i]) {Fl[i]=1,Q[i]=C;cnt++;return;}C^=Q[i];}}
int main(){
	freopen("present.in","r",stdin);freopen("present.out","w",stdout);
	RI i,j,h;scanf("%d%d",&n,&m);k=2*m+n-2;if(n==1) {puts("1");return 0;} for(i=1;i<=m;i++) F[d(1,i)][i]=1,F[d(2,i)][i+m]=1;for(i=3;i<=n;i++) F[d(i,1)][i-2+2*m]=1;
	for(i=3;i<=n;i++)for(j=2;j<=m;j++)GA(i-2,j-1),F[d(i,j)]=C;for(i=1;i<=n-2;i++)GA(i,m),Ins();for(i=1;i<=m;i++) GA(n-1,i),Ins(),GA(n,i),Ins();for(i=cnt+1;i<=k;i++) Ans=Ans*2%mod;printf("%lld\n",Ans);
}

标签:891,int,Ybtoj,long,Ny,Fl,RI,高斯消,define
来源: https://www.cnblogs.com/275307894a/p/15866330.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有