着实被批评了一下,自己的数论确实太烂了. 题意:一条路上,有n个炸弹,给出每个炸弹的位置,一次走一步的概率是p,走两步的概率是1-p.求安全走完的概率. 定义dp[i] = dp[i-1]*p + dp[i-2]*(1-p) 由于路很长,不能递推. n个炸弹就把路分成了n+1段路,计算走完一段路安全的概率就行,然后相乘. | p ,  1 | | dp[i-1],dp[i-2] | * | 1-p,0 |   = | dp[i],dp[i-1] | 求安全走完一段长为L路的概率,就将矩阵乘L次即…

任意门:http://poj.org/problem?id=1830 开关问题 Time Limit: 1000MS Memory Limit: 30000K Total Submissions: 10742 Accepted: 4314 Description 有N个相同的开关,每个开关都与某些开关有着联系,每当你打开或者关闭某个开关的时候,其他的与此开关相关联的开关也会相应地发生变化,即这些相联系的开关的状态如果原来为开就变为关,如果为关就变为开.你的目标是经过若干次开关操作后使得最后N个开…

题目地址:id=2947">POJ 2947 题意:N种物品.M条记录,接写来M行,每行有K.Start,End,表述从星期Start到星期End,做了K件物品.接下来的K个数为物品的编号. 此题注意最后结果要调整到3-9之间. 思路: 非常easy想到高斯消元. 可是是带同余方程式的高斯消元,開始建立关系的时候就要MOD 7 解此类方程式时最后求解的过程要用到扩展gcd的思想,举个样例,假设最后得到的矩阵为:     1  1   4     0  6   4    则6 * X2 %…

偶数个相邻, 以n*m个点为变量, 建立异或方程组然后高斯消元... O((n*m)^3)复杂度看起来好像有点大...但是压一下位的话就是O((n*m)^3 / 64), 常数小, 实际也跑得很快. ------------------------------------------------------------------------------------- #include<cstdio> #include<bitset> #include<cstring>…

矩阵运算: \(A\times B\)叫做\(A\)左乘\(B\),或者\(B\)右乘\(A\). 行列式性质: \(1.\)交换矩阵的两行(列),行列式取相反数. \(2.\)某一行元素都\(\times k\),行列式值也\(\times k\). \(3.\)某一行加到另一行上,行列式值不变. \(4.\)矩阵某两行(列)元素分别成比例,行列式值为\(0\). \(5.A+B=C\Rightarrow|A|+|B|=|C|\). \(6.\)矩阵与转置矩阵行列式相等. 对于方阵而言: \(…

Description 我们称一个由0和1组成的矩阵是和谐的,当且仅当每个元素都有偶数个相邻的1.一个元素相邻的元素包括它本 身,及他上下左右的4个元素(如果存在). 给定矩阵的行数和列数,请计算并输出一个和谐的矩阵.注意:所有元素为0的矩阵是不允许的. Input 输入一行,包含两个空格分隔的整数m和n,分别表示矩阵的行数和列数. Output 输出包含m行,每行n个空格分隔整数(0或1),为所求矩阵.测试数据保证有解. Sample Input 4 4 Sample Output 0 1 0…

题目链接 高斯消元其实是个大模拟qwq 所以就着代码食用 首先我们读入 ;i<=n;++i) ;j<=n+;++j) scanf("%lf",&s[i][j]); 读入肯定没什么问题(不过我在这卡了一分多钟) 然后我们要进行消元操作 所谓消元操作其实就是对于输入的矩阵 比如说 9 3 2 2 1 4 7 3 1 3 4 5 进行一番乱搞,使得第当前枚举的(比如说枚举第i行第i列)s[i][j]系数变成1. 实际上就是整行同除qwq 比如我们除完第一行第一列的之后,矩…

pro:给定N*M的矩阵,以及初始玩家位置. 规定玩家每次会等概率的向左走,向右走,向下走,原地不动,问走到最后一行的期望.保留4位小数. sol:可以列出方程,高斯消元即可,发现是三角矩阵,O(N*M)----元素个数.  也可以用反复逼近答案. 反复做,dp[i][j]=(dp[i][j+1]+dp[i][j-1]+dp[i][j]+dp[i-1][j])/d[j]+1.0  为了使逼近效果更好,我每次先左一次,再右一次. #include<bits/stdc++.h> #define r…

题意:构造一个$n*m$矩阵 使得每个元素和上下左右的xor值=0 题解:设第一行的每个元素值为未知数 可以依次得到每一行的值 然后把最后一行由题意条件 得到$m$个方程 高斯消元解一下 bitset写起来比较方便 #include <stdio.h> #include <iostream> #include <algorithm> #include <bitset> using namespace std; const int MAXN = 45; int…

看题就知道要使用高斯消元求解! 代码如下: #include<iostream> #include<algorithm> #include<iomanip> #include<cmath> #include<cstring> using namespace std; ][],p,ans[]; ]; int pows(int a,int b) { ; while(b){ ) ans=(ans*a)%p; b>>=; a=(a*a)%p;…