线性代数五之高斯消元——[SDOI2010]外星千足虫，[HNOI2013]游走，[HNOI2011]XOR和路径，[hdu 4035]Maze

最新推荐文章于 2023-03-12 15:18:42 发布
最新推荐文章于 2023-03-12 15:18:42 发布
阅读量234 收藏
点赞数 1
分类专栏： 数学 # 高斯消元 文章标签： 数学 高斯消元
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Emm_Titan/article/details/119154604
版权

多类型高斯消元杂题

[SDOI2010]外星千足虫

description

solution

高斯消元的模板题

虽然感觉问了个最早确定所有虫的时间戳，但并没有什么用

在高斯消元过程中，使用的最大行编号就是最早确定时间

code

#include <cstdio>
#include <bitset>
#include <iostream>
using namespace std;
#define maxn 1005
int n, m;
bitset < maxn > a[maxn << 1];

int main() {
	scanf( "%d %d", &n, &m );
	for( int i = 1, x;i <= m;i ++ ) {
		for( int j = 1;j <= n;j ++ )
			scanf( "%1d", &x ), a[i][j] = x;
		scanf( "%d", &x ), a[i][n + 1] = x;
	}
	int ans = 0;
	for( int i = 1;i <= n;i ++ ) {
		int row = i;
		while( row <= m && ! a[row][i] ) row ++;
		if( row == m + 1 ) return ! printf( "Cannot Determine\n" );
		if( i ^ row ) swap( a[row], a[i] );
		ans = max( ans, row );
		for( int j = 1;j <= m;j ++ )
			if( i == j || ! a[j][i] ) continue;
			else a[j] ^= a[i];
	}
	printf( "%d\n", ans );
	for( int i = 1;i <= n;i ++ )
		if( a[i][n + 1] ) printf( "?y7M#\n" );
		else printf( "Earth\n" );
	return 0;
}

[HNOI2013]游走

description

solution

将路径拆分成每条边期望经过次数乘以边权的求和

贪心的把最大权值放在期望经过次数最少的边上

但是边的级别是空间时间不足以承受的

事实上，边的期望至于边连接的两个端点有关

每个点到与之相连边的概率一样，期望一样

所以到某条特定边的期望就是经过该点期望除以该点连接的总边数

一条边被经过的期望则是两个端点到这条边的期望和

于是乎求边的期望就转化为求点的期望

而点的期望至于相邻点有关

设 E i : i E_i:i Ei​:i点期望， n u m i : i num_i:i numi​:i点总边数， x 1 , x 2 , . . . , x k x_1,x_2,...,x_k x1​,x2​,...,xk​与 i i i相连，则有 E i = ∑ j = 1 k E j n u m j E_i=\sum_{j=1}^k\frac{E_j}{num_j} Ei​=∑j=1k​numj​Ej​​

每个点都能列出一个这样的方程，高斯消元解决

特别地，游走是从点 1 1 1开始的，所以期望要 + 1 +1 +1，游走到点 n n n就不会继续进行了，所以计算期望时不能纳入考虑

code

#include <cstdio>
#include <vector>
#include <iostream>
#include <algorithm>
using namespace std;
#define maxn 505
#define eps 1e-6
#define maxm 125005
vector < int > G[maxn]; 
int n, m;
int num[maxn], from[maxm], to[maxm];
double E[maxm];
double x[maxn][maxn];

double Fabs( double x ) {
	return x < 0 ? -x : x;
}

int main() {
	scanf( "%d %d", &n, &m );
	for( int i = 1, u, v;i <= m;i ++ ) {
		scanf( "%d %d", &u, &v );
		G[u].push_back( v );
		G[v].push_back( u );
		from[i] = u, to[i] = v;
		num[u] ++, num[v] ++;
	}
	x[1][n] = 1;
	for( int i = 1;i < n;i ++ ) {
		x[i][i] = 1;
		for( auto j : G[i] )
			if( j != n )
				x[i][j] = -1.0 / num[j];
	}
	for( int i = 1;i < n;i ++ ) {
		int k = i;
		for( int j = i + 1;j < n;j ++ )
			if( Fabs( x[j][i] ) > Fabs( x[k][i] ) )
				k = j;
		if( i ^ k ) swap( x[i], x[k] );
		for( int j = n;j >= i;j -- )
			x[i][j] /= x[i][i];
		for( int j = 1;j < n;j ++ )
			if( i ^ j )
				for( int k = n;k >= i;k -- )
					x[j][k] -= x[j][i] * x[i][k];
	}
	for( int i = 1;i <= m;i ++ ) {
		if( from[i] != n )
			E[i] += x[from[i]][n] / num[from[i]];
		if( to[i] != n )
			E[i] += x[to[i]][n] / num[to[i]];
	}
	sort( E + 1, E + m + 1 );
	double ans = 0;
	for( int i = 1;i <= m;i ++ )
		ans += E[i] * ( m - i + 1 );
	printf( "%.3f", ans );
	return 0;
}

[HNOI2011]XOR和路径

description

solution

一般看到充满二进制意味的异或都要考虑拆位

此题也不例外，每一个单独计算期望

设 f i : f_{i}: fi​: 表示 i → n i\rightarrow n i→n路径该位为 1 1 1的概率，则 1 − f i 1-f_i 1−fi​表示 i → n i\rightarrow n i→n路径为 0 0 0的概率， d i : i d_i:i di​:i的出度

∀ ( u , v ) ∈ E d g e f u = 1 d u ( ∑ w ( u , v ) = 0 f v + ∑ w ( u , v ) = 1 ( 1 − f v ) ) ⇔ ∀ ( u , v ) ∈ E d g e d u × f u = ( ∑ w ( u , v ) = 0 f v + ∑ w ( u , v ) = 1 ( 1 − f v ) ) \forall_{(u,v)\in Edge}f_u=\frac{1}{d_u}\Big(\sum_{w(u,v)=0}f_v+\sum_{w(u,v)=1}(1-f_v)\Big)\Leftrightarrow \forall_{(u,v)\in Edge}d_u\times f_u=\Big(\sum_{w(u,v)=0}f_v+\sum_{w(u,v)=1}(1-f_v)\Big) ∀(u,v)∈Edge​fu​=du​1​(∑w(u,v)=0​fv​+∑w(u,v)=1​(1−fv​))⇔∀(u,v)∈Edge​du​×fu​=(∑w(u,v)=0​fv​+∑w(u,v)=1​(1−fv​))
在这里插入图片描述

换言之， u → n u\rightarrow n u→n的路径想要为 1 : 1: 1:如果 ( u → v ) (u\rightarrow v) (u→v)这条边为 1 1 1，必须 v → n v\rightarrow n v→n的路径为 0 0 0，否则为 1 1 1

方程化为

d u f u − ∑ w ( u , v ) = 0 f v + ∑ w ( u , v ) = 1 f u = ∑ w ( u , v ) = 1 1 d_uf_u-\sum_{w(u,v)=0}f_v+\sum_{w(u,v)=1}f_u=\sum_{w(u,v)=1}1 du​fu​−∑w(u,v)=0​fv​+∑w(u,v)=1​fu​=∑w(u,v)=1​1

高斯消元， a n s = ∑ i 2 i f i ( 1 ) ans=\sum_i2^if_i(1) ans=∑i​2ifi​(1)

code

#include <cstdio>
#include <vector>
#include <cstring>
#include <iostream>
using namespace std;
#define eps 1e-8
#define maxn 105
vector < pair < int, int > > G[maxn];
int n, m;
int d[maxn];
double ans[maxn];
double a[maxn][maxn];

double Fabs( double x ) {
	return x < 0 ? -x : x;
}

void build( int x ) {
	a[n][n] = 1;
	for( int i = 1;i < n;i ++ ) {
		a[i][i] = d[i];
		for( auto j : G[i] ) {
			if( j.second & x ) a[i][j.first] ++, a[i][n + 1] ++;
			else a[i][j.first] --;		
		}
	}
}

void gauss() {
	for( int i = 1;i <= n;i ++ ) {
		int k = i;
		for( int j = i;j <= n;j ++ )
			if( Fabs( a[j][i] ) > Fabs( a[k][i] ) )
				k = j;
		if( k ^ i ) swap( a[i], a[k] );
		for( int j = i + 1;j <= n;j ++ )
			if( Fabs( a[j][i] ) < eps ) continue;
			else {
				double t = a[j][i] / a[i][i];
				for( k = i;k <= n + 1;k ++ )
					a[j][k] -= a[i][k] * t;
				a[j][i] = 0;
			}
	}
	for( int i = n;i;i -- ) {
		for( int j = i + 1;j <= n;j ++ )
			a[i][n + 1] -= a[i][j] * ans[j];
		ans[i] = a[i][n + 1] / a[i][i];
	}
	memset( a, 0, sizeof( a ) );
}

int main() {
	scanf( "%d %d", &n, &m );
	int maxx = 0;
	for( int i = 1, u, v, w;i <= m;i ++ ) {
		scanf( "%d %d %d", &u, &v, &w );
		G[u].push_back( make_pair( v, w ) );
		d[u] ++;
		if( u ^ v ) {
			G[v].push_back( make_pair( u, w ) );
			d[v] ++;
		}
		maxx = max( maxx, w );
	}
	double ret = 0;
	for( int i = 1;i <= maxx;i <<= 1 ) {
		build( i );
		gauss();
		ret += ans[1] * i;
	}
	printf( "%.3f\n", ret );
	return 0;
}

Maze(树上高斯消元)

problem

solution

令 E i : i E_i:i Ei​:i节点逃出期望经过的边数， e d g e i : i edge_i:i edgei​:i相连边数
∀ i , i ∈ l e a f E i = k i × E 1 + e i × 0 + ( 1 − k i − e i ) × ( E f a i + 1 ) = k i × E 1 + ( 1 − k i − e i ) × E f a i + ( 1 − k i − e i ) \forall_{i,i\in leaf}E_i=k_i\times E_1+e_i\times 0+(1-k_i-e_i)\times (E_{fa_i}+1)\\=k_i\times E_1+(1-k_i-e_i)\times E_{fa_i}+(1-k_i-e_i) ∀i,i∈leaf​Ei​=ki​×E1​+ei​×0+(1−ki​−ei​)×(Efai​​+1)=ki​×E1​+(1−ki​−ei​)×Efai​​+(1−ki​−ei​)

∀ i , i ∉ l e a f E i = k i × E 1 + e i × 0 + ( 1 − k i − e i ) × 1 e d g e i × ( ∑ j ∈ s o n i ( E j + 1 ) + E f a i + 1 ) = k i × E 1 + 1 − k i − e i e d g e i × E f a i + 1 − k i − e i m × ∑ j ∈ s o n i E j + ( 1 − k i − e i ) \forall_{i,i\notin leaf}E_{i}=k_i\times E_1+e_i\times 0+(1-k_i-e_i)\times\frac{1}{edge_i}\times\Big(\sum_{j\in son_i}(E_j+1)+E_{fa_i}+1\Big)\\=k_i\times E_1+\frac{1-k_i-e_i}{edge_i}\times E_{fa_i}+\frac{1-k_i-e_i}{m}\times\sum_{j\in son_i}E_j+(1-k_i-e_i) ∀i,i∈/​leaf​Ei​=ki​×E1​+ei​×0+(1−ki​−ei​)×edgei​1​×(j∈soni​∑​(Ej​+1)+Efai​​+1)=ki​×E1​+edgei​1−ki​−ei​​×Efai​​+m1−ki​−ei​​×j∈soni​∑​Ej​+(1−ki​−ei​)

转移涉及父亲儿子且彼此依赖，而我们非常想得到不交叉的线性递推关系

想办法变成只跟父亲有关的转移，形式的设 E i = A i × E 1 + B i × E f a i + C i E_i=A_i\times E_1+B_i\times E_{fa_i}+C_i Ei​=Ai​×E1​+Bi​×Efai​​+Ci​

类比解方程组，列出以下式子

∀ i , i ∉ l e a f \forall_{i,i\notin leaf} ∀i,i∈/​leaf​
∑ j ∈ s o n i E j = ∑ j ∈ s o n i A j × E 1 + B j × E f a j + C j = ∑ j ∈ s o n i A j × E 1 + B j × E i + C j \sum_{j\in son_i}E_j=\sum_{j\in son_i}A_j\times E_1+B_j\times E_{fa_j}+C_j=\sum_{j\in son_i}A_j\times E_1+B_j\times E_i+C_j\\ j∈soni​∑​Ej​=j∈soni​∑​Aj​×E1​+Bj​×Efaj​​+Cj​=j∈soni​∑​Aj​×E1​+Bj​×Ei​+Cj​

E i = k i E 1 + 1 − k i − e i e d g e i E f a i + 1 − k i − e i m ∑ j ∈ s o n i E j + ( 1 − k i − e i ) = k i E 1 + 1 − k i − e i e d g e i E f a i + 1 − k i − e i e d g e i ∑ j ∈ s o n i ( A j E 1 + B j E i + C j ) + ( 1 − k i − e i ) E_i=k_iE_1+\frac{1-k_i-e_i}{edge_i}E_{fa_i}+\frac{1-k_i-e_i}{m}\sum_{j\in son_i}E_j+(1-k_i-e_i)\\=k_iE_1+\frac{1-k_i-e_i}{edge_i}E_{fa_i}+\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}(A_jE_1+B_jE_i+C_j)+(1-k_i-e_i) Ei​=ki​E1​+edgei​1−ki​−ei​​Efai​​+m1−ki​−ei​​j∈soni​∑​Ej​+(1−ki​−ei​)=ki​E1​+edgei​1−ki​−ei​​Efai​​+edgei​1−ki​−ei​​j∈soni​∑​(Aj​E1​+Bj​Ei​+Cj​)+(1−ki​−ei​)
在这里插入图片描述

⇔ \Leftrightarrow ⇔
( 1 − 1 − k i − e i e d g e i ∑ j ∈ s o n i B j ) E i = (1-\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}B_j)E_i= (1−edgei​1−ki​−ei​​j∈soni​∑​Bj​)Ei​= ( k i + 1 − k i − e i e d g e i ∑ j ∈ s o n i A j ) E 1 + 1 − k i − e i e d g e i E f a i + ( 1 − k i − e i ) + 1 − k i − e i e d g e i ∑ j ∈ s o n i C j (k_i+\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}A_j)E_1+\frac{1-k_i-e_i}{edge_i}E_{fa_i}+(1-k_i-e_i)+\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}C_j (ki​+edgei​1−ki​−ei​​j∈soni​∑​Aj​)E1​+edgei​1−ki​−ei​​Efai​​+(1−ki​−ei​)+edgei​1−ki​−ei​​j∈soni​∑​Cj​
⇒ A i = k i + 1 − k i − e i e d g e i ∑ j ∈ s o n i A j ( 1 − 1 − k i − e i e d g e i ∑ j ∈ s o n i B j ) \Rightarrow A_i=\frac{k_i+\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}A_j}{(1-\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}B_j)} ⇒Ai​=(1−edgei​1−ki​−ei​​∑j∈soni​​Bj​)ki​+edgei​1−ki​−ei​​∑j∈soni​​Aj​​ B i = 1 − k i − e i e d g e i ( 1 − 1 − k i − e i e d g e i ∑ j ∈ s o n i B j ) B_i=\frac{\frac{1-k_i-e_i}{edge_i}}{(1-\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}B_j)} Bi​=(1−edgei​1−ki​−ei​​∑j∈soni​​Bj​)edgei​1−ki​−ei​​​ C i = ( 1 − k i − e i ) + 1 − k i − e i e d g e i ∑ j ∈ s o n i C j ( 1 − 1 − k i − e i e d g e i ∑ j ∈ s o n i B j ) C_i=\frac{(1-k_i-e_i)+\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}C_j}{(1-\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}B_j)} Ci​=(1−edgei​1−ki​−ei​​∑j∈soni​​Bj​)(1−ki​−ei​)+edgei​1−ki​−ei​​∑j∈soni​​Cj​​
在这里插入图片描述

∀ i , i ∈ l e a f \forall_{i,i\in leaf} ∀i,i∈leaf​
E i = k i × E 1 + ( 1 − k i − e i ) × E f a i + ( 1 − k i − e i ) = A i × E 1 + B i × E f a i + C i ⇒ { A i = k i B i = 1 − k i − e i C i = 1 − k i − e i E_i=k_i\times E_1+(1-k_i-e_i)\times E_{fa_i}+(1-k_i-e_i)=A_i\times E_1+B_i\times E_{fa_i}+C_i\\ \Rightarrow \begin{cases}A_i=k_i\\B_i=1-k_i-e_i\\C_i=1-k_i-e_i\end{cases} Ei​=ki​×E1​+(1−ki​−ei​)×Efai​​+(1−ki​−ei​)=Ai​×E1​+Bi​×Efai​​+Ci​⇒⎩⎪⎨⎪⎧​Ai​=ki​Bi​=1−ki​−ei​Ci​=1−ki​−ei​​
从叶子节点开始倒着往上递推（也被称之为树上高斯消元）

树上高斯消元

把有关父亲和儿子的递推式通过 A i x + B i y + C i A_ix+B_iy+C_i Ai​x+Bi​y+Ci​的形式转化为只由一方单向线性递推

然后正着/倒着递推消元

直到算出 A 1 , B 1 , C 1 A_1,B_1,C_1 A1​,B1​,C1​

E 1 = A 1 × E 1 + B 1 × 0 + C 1 ⇒ E 1 = C 1 1 − A 1 E_1=A_1\times E_1+B_1\times 0+C_1\Rightarrow E_1=\frac{C_1}{1-A_1} E1​=A1​×E1​+B1​×0+C1​⇒E1​=1−A1​C1​​

当 A 1 A_1 A1​无限趋近于 1 1 1时，无解

如果 i i i的 A i , B i , C i A_i,B_i,C_i Ai​,Bi​,Ci​的分母形式 1 − 1 − k i − e i e d g e i ∑ j ∈ s o n i B j 1-\frac{1-k_i-e_i}{edge_i}\sum_{j\in son_i}B_j 1−edgei​1−ki​−ei​​∑j∈soni​​Bj​趋近于 0 0 0，也是无解

注意精度问题，eps=1e-8都不行，因为涉及 / 100 /100 /100

code

#include <cmath>
#include <cstdio>
#include <vector>
using namespace std;
#define eps 1e-10
#define maxn 10005
vector < int > G[maxn];
int T, n;
double k[maxn], e[maxn], A[maxn], B[maxn], C[maxn];

bool dfs( int i, int fa ) {
	if( G[i].size() == 1 && G[i][0] == fa ) {
		A[i] = k[i], B[i] = C[i] = 1 - k[i] - e[i];
		return 1;
	}
	A[i] = k[i];
	B[i] = ( 1 - k[i] - e[i] ) / G[i].size();
	C[i] = 1 - k[i] - e[i];
	double t = 0;
	for( auto j : G[i] ) {
		if( j == fa ) continue;
		else if( ! dfs( j, i ) ) return 0;
		A[i] += A[j] * B[i];
		C[i] += C[j] * B[i];
		t += B[j] * B[i];
	}
	if( fabs( 1 - t ) < eps ) return 0;
	else {
		A[i] /= ( 1 - t );
		B[i] /= ( 1 - t );
		C[i] /= ( 1 - t );
		return 1;
	}
}

int main() {
	scanf( "%d", &T );
	for( int Case = 1;Case <= T;Case ++ ) {
		scanf( "%d", &n );
		for( int i = 1;i <= n;i ++ )
			G[i].clear();
		for( int i = 1, u, v;i < n;i ++ ) {
			scanf( "%d %d", &u, &v );
			G[u].push_back( v );
			G[v].push_back( u );  
		}
		for( int i = 1;i <= n;i ++ ) {
			scanf( "%lf %lf", &k[i], &e[i] );
			k[i] /= 100, e[i] /= 100;
		}
		if( dfs( 1, 0 ) && fabs( 1 - A[1] ) > eps )
			printf( "Case %d: %f\n", Case, C[1] / ( 1 - A[1] ) );
		else
			printf( "Case %d: impossible\n", Case );
	}
	return 0;
}
ikrvxt
关注
专栏目录
bzoj [Sdoi2010]外星千足虫
Elijahqi
04-08 232
http://www.elijahqi.win/archives/2976 题目描述 公元2089年6月4日，在经历了17年零3个月的漫长旅行后，“格纳格鲁一号”载人火箭返回舱终于安全着陆。此枚火箭由美国国家航空航天局（NASA）研制发射，行经火星、金星、土卫六、木卫二、谷神星、“张衡星”等23颗太阳系星球，并最终在小行星“杰森星”探寻到了地外生命。宇航员在“杰森星”地表岩层下45.70米位...
高斯消元小结
m0_49959202的博客
09-19 306
文章目录线性空间1.算法分析1.1 高斯消元1.2 线性基2.模板2.1 高斯消元2.1.1 解累加方程2.1.1.1 整数系数2.1.1.2 浮点数系数2.1.2 解异或方程2.2 线性基3.例题3.1 高斯消元3.2 线性基 线性空间 1.算法分析 1.1 高斯消元 模拟线性代数的运算 1.2 线性基      线性基的本质是用一个数x的尽可能高位的1来代表x     线性基是向量空间的一组基，通常可以解决有关异或的一些题目。     通俗一点的讲法就是由一个集合构造出来的另一个集合，它有以下几个性质：
【bzoj1923】[Sdoi2010]外星千足虫 高斯消元解xor方程组
DQS的树洞
03-17 812
DescriptionInput第一行是两个正整数 N, M。 接下来 M行，按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包含一个“01”串和一个数字，用一个空格隔开。“01”串按位依次表示每只虫 子是否被放入机器：如果第 i 个字符是“0”则代表编号为 i 的虫子未被放入，“1” 则代表已被放入。后面跟的数字是统计的昆虫足数 mod 2 的结果。 由于 NASA的实验机器精确
[SDOI2010]外星千足虫（高斯消元）
wang_zitao的博客
10-29 273
洛谷2447 模2相当于异或运算，将高斯消元的板子改一改就行了。 注意到n的范围到1000，于是可以考虑用bitset优化。 消元过程中若发现某一列上每行均为0，可直接判断无解，否则记录下为1的在当前行之后最小的一行，对每一列求最大值，即为所需方程个数。 #include&amp;amp;lt;bitset&amp;amp;gt; #include&amp;amp;lt;cstdio&amp;amp;gt; #include&amp;amp;lt;cstdlib&am
P3211-[HNOI2011]XOR和路径【高斯消元】
QuantAsk
03-03 165
正题 题目链接:https://www.luogu.com.cn/problem/P3211 题目大意 一个nnn个点mmm条边的无向图，从111到nnn随机游走。求期望路径异或和。 2≤n≤100,1≤m≤1042\leq n\leq 100,1\leq m\leq 10^42≤n≤100,1≤m≤104 解题思路 因为是异或的期望，很难直接处理，所以考虑按位考虑每一位是111的概率。 然后nnn很小就是一个很显然的高斯消元了。设fif_ifi​表示i∼ni\sim ni∼n是111的概率。 fx=1
[HNOI2011]XOR和路径——高斯消元、随机游走模型
偶耶的博客
07-22 321
就我这未优化的大常数代码能达到rank1？？
millipede-elisp:打印美丽的千足虫
05-17
在实际使用时，用户需要将`millipede.el`添加到Emacs的加载路径中，并通过调用相应的Elisp函数来启动千足虫。此外，还可以通过修改配置文件来自定义千足虫的颜色、速度、节段数量等属性，以满足个人喜好。 总的来说...
artxgj#twitter-chinese-text#千足虫2
07-25
千足虫Search mdbg for definitionTweets containing 千足虫2021-12-19 03:13:05 UTC ~ 联合早报
millipede-haskell:打印美丽的千足虫
05-09
千足虫 :bug: 打印美丽的千足虫 用法 # prints a millipede of species 0 $ millipede # prints a millipede of species n, scientists have discovered 1309 species of millipedes to date $ millipede <n>
四川省成都石室中学2013届高三生物下学期“三诊”模拟考试试题新人教版
08-19
- 共同进化是指不同物种间相互影响并共同演化的过程，此蜣螂与千足虫间的互动即属此例。 - 种群基因频率的改变是物种进化的基本特征。 - 生殖隔离的缺乏意味着这类蜣螂仍能与普通蜣螂交配，未形成新物种。 8. ...
[SDOI2010]外星千足虫
weixin_30847939的博客
12-22 85
传送门 题目要求解线性方程组，不过这次是要在模2的意义下做。也就是相当于把变量互相异或。 本题用\(O(n^3)\)的高斯消元过不了，不过在互相异或的情况下，消元其实就是两个方程组直接互相异或，用bitset优化一下即可。 至于如何求最少使用的方程，因为每次都是找一个最小的，使得当前位为1的方程，所以这个就直接符合使用最少的方程，无需特殊处理。 #include<cstdio> #in...
P2447 [SDOI2010]外星千足虫
qq_35975367的博客
08-04 228
P2447 [SDOI2010]外星千足虫 题意： 有n个未知数 给你一个m行n+1列的式子，对于每行，1到n列为这个n个未知数的系数，第n+1列为该行式子的和mod2，问n个未知数是否有唯一解，并输出，并输出最少需要前k个式子就可以得到唯一解。有多个解输出Cannot Determine 题解： 很明显高斯消元异或方程组，没错这算是个模板题 这个k怎么求？在高斯消元求解异或方程组的过程中是有交换行的，我们只需要在交换行的过程中记录答案就行 复杂度是O(13∗n2∗m)O(\frac{1}{3}*n^2*m
P2447 [SDOI2010] 外星千足虫（异或高斯消元 + bitset运用）
最新发布
zhi6fui的博客
03-12 287
现在你面前摆有1…N 编号的 N 只千足虫，你的任务是鉴定每只虫子所属的星球，但不允许亲自去数它们的足。Charles 每次会在这 N 只千足虫中选定若干只放入“昆虫点足机”（the Insect Feet Counter, IFC）中，“点足机”会自动统计出其内所有昆虫足数之和。Charles 会将这个和数 \bmodmod 22 的结果反馈给你，同时告诉你一开始放入机器中的是哪几只虫子。他的这种统计操作总共进行 M。
Bzoj 1923: [Sdoi2010]外星千足虫(bitset优化高斯消元)
_zidaoziyan的博客
09-13 675
传送门: 题意:给你一些虫子,每个有奇数条足或者偶数条足,有n个虫子,进行m次试验, 每次试验有一个字符串,后面跟一个字符,字符串中位置为0表示这个虫子没有 被放入,为1则相反,最后这几个虫子的足之和对2取模为后面跟着的字符, 求最后每个虫子的奇偶或者不确定,最早在第几次试验的时候可以确定 思路:高斯消元 #include<bitset> #include<algorithm>
[BZOJ1923][SDOI2010]外星千足虫
abcyan1235的博客
04-10 104
bzoj luogu 题意 给你\(m\)个\(n\)元异或方程，方程是一个一个给的，问你给到第几个的时候出解，或者是在得到\(m\)个方程后仍不能确定解。保证方程不会前后矛盾，即不会出现无解。 sol 我今天突然发现异或方程组的高斯消元和线性基就是一个东西。 这道题其实就是问你什么时候你可以得到\(n\)个线性无关方程。 所以你可以把每个方程组类似线性基那样插入，每次取异或。 一旦得到...
[BZOJ1923][Sdoi2010]外星千足虫（高斯消元）
zyf2000
01-12 667
题目描述传送门题解高斯消元求解异或方程组 和普通消元的思路差不多，同样是找到一个可用的交换，然后把原来的加减运算全部改为异或 如果a(i,i)无法交换成0的话，说明这里有可能出现无解或多解，那么回代的时候判断一下剩下的常数项和b(i)是否相等就可以了 不过这道题保证一定有解，所以出现上面的情况就一定有多解，直接返回就行了 注意交换的时候记录一下最多用到的方程k 安利一个讲解：解线性方程组—
bzoj 1923: [Sdoi2010]外星千足虫 （高斯消元）
clover_hxy的博客
04-23 955
1923: [Sdoi2010]外星千足虫 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 753  Solved: 474 [Submit][Status][Discuss] Description Input 第一行是两个正整数 N, M。 接下来 M行，按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包
【无码专区11】异或2（结论 / 推式子 + 哈希hash + 大整数高精度 加减乘除重载考察）
热门推荐
11-13 11万+
直接写吐了，真的nm难调啊
【无码专区2】序列划分（数学）
10-19 10万+
天啦，只看20位！！
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值