Loading

随笔分类 -  微积分，线性代数，群论

POJ3735 Training little cats 矩阵快速幂 + 矩阵快速幂优化
摘要：有n只猫，每只猫初始石刻没糖果。 进行m次重复的k次操作。 g i 给第i只猫一个糖果 s i j 交换第i只猫和第j只猫的糖果 c i 吃掉第i只猫的所有糖果 m <=1e9 明显是用矩阵快速幂处理。 难点在于建立转移矩阵。 注意到只需对单位矩阵操作。对于n只猫建立 n+1 * 1 的矩阵，多的一 阅读全文
posted @ 2020-07-21 14:04 MQFLLY 阅读(146) 评论(0) 推荐(0)
POJ 3185 The Water Bowls 翻转计数/异或高斯消元
摘要：给定20个桶的状态，可以翻转桶把相邻的状态变换。问最少的变换次数。 此类问题通常有以下解决方案 1.暴力枚举 2.高斯消元 + 暴力枚举 3.找规律性质 这题和上一题POJ1830 很相似，可以用高斯消元异或方程组枚举自由元做。 但是比较麻烦； 注意到每个桶至多翻转一次。 从右往左的第一个1必然要通 阅读全文
posted @ 2020-07-17 10:46 MQFLLY 阅读(177) 评论(0) 推荐(0)
POJ1830 开关问题 高斯消元 异或
摘要：设xixi表示第ii个开关的操作情况,那么xi=1xi=1为按了这个开关,否则xi=0xi=0表示这个开关并没有按下去,那么我们接着统计ai,jai,j表示第i个开关和第j个开关的联系情况,ai,j=1ai,j=1表示按下j会影响i的状态,否则就是不会影响,即ai,j=0ai,j=0,特殊地一个点就 阅读全文
posted @ 2020-07-17 09:07 MQFLLY 阅读(123) 评论(0) 推荐(0)
高斯消元模板
摘要：double a[maxn][maxn]; int n; void Gauss() { for (int i = 1; i <= n; i++) { int Max = i; for (int j = i + 1; j <= n; j++) { if (fabs(a[j][i]) > fabs(a[ 阅读全文
posted @ 2020-07-17 09:03 MQFLLY 阅读(155) 评论(0) 推荐(0)
HDU - 5950 Recursive sequence
摘要：题意：f[1]=a,f[2]=b,f[i]=f[i-1]+2*f[i-2]+i^4(I>=3),求f[n]%mod 明显的矩阵快速幂。 这是我一开始的矩阵 这是结果的 还是太年轻了 struct Mat { ll m[7][7]; }; Mat E; Mat P; void init() { for 阅读全文
posted @ 2020-07-04 18:37 MQFLLY 阅读(141) 评论(0) 推荐(0)
HDU 4686 Arc of Dream 矩阵快速幂
摘要：若能根据关系 推出转移矩阵 就是简单题 然而这需要多做题的经验 ll Ax, Ay, A0, B0 ,Bx,By; struct Mat { ll m[5][5]; }; Mat E; Mat P; void init() { for (int i = 0; i < 5; i++) E.m[i][i 阅读全文
posted @ 2020-07-03 15:54 MQFLLY 阅读(136) 评论(0) 推荐(0)
模板-高斯消元
摘要：高斯-约旦消元法 int n; double a[maxn][maxn]; void Gauss() { scanf("%d", &n); for (int i = 1; i <=n; i++) { for (int j = 1; j <=n + 1; j++) { scanf("%lf", &a[ 阅读全文
posted @ 2020-03-12 12:33 MQFLLY 阅读(122) 评论(0) 推荐(0)