Description 一个黑白网格,点一次会改变这个以及与其连通的其他方格的颜色,求最少点击次数使得所有全部变成黑色. Sol 高斯消元解异或方程组. 先建立一个方程组. \(x_i\) 表示这个点是否被用过. 因为第二次使用同一个点,这个点的贡献就被消除了,所以每个点只会被用 0/1 次. \( ...
分类：
其他好文 时间：
2016-11-29 23:11:26
阅读次数：
208
科学计算 最小二乘leastsq 非线性方程组求解 from scipy.optimize import fsolve 可以传入雅可比行列式的方式来提高运算速度，如result=fsolve(f,[1,1,1],fprime=j) 插值 from scipy import interpolate 积 ...
分类：
其他好文 时间：
2016-11-21 07:46:14
阅读次数：
335
用 GSL 求解超定方程组及矩阵的奇异值分解(SVD)最近在学习高动态图像（HDR）合成的算法，其中需要求解一个超定方程组，因此花了点时间研究了一下如何用 GSL 来解决这个问题。GSL 里是有个最小...
分类：
其他好文 时间：
2016-11-20 14:06:14
阅读次数：
346
一、鸡兔同笼 鸡和兔子35只，腿一共有94条，求出鸡和兔子的数量各是多少？（鸡的数目是23，兔子的数目是12） $n = 0; for($ i=1;$i<35;$i++){ $n ++; //n=34; $ j=35-i; if(i*2+j*4==94){ printf("鸡的数目是%d只, 兔子的 ...
分类：
其他好文 时间：
2016-11-19 07:43:11
阅读次数：
168
【题目链接】 http://acm.hdu.edu.cn/showproblem.php?pid=5446 【题目大意】 给出一个合数M的每一个质因子，同时给出n,m，求C(n,m)%M。 【题解】 首先我们可以用Lucas定理求出对答案对每个质因子的模，然后我们发现只要求解这个同余方程组就可以得到 ...
分类：
其他好文 时间：
2016-11-17 18:33:19
阅读次数：
251
题意： 有 解法：中国剩余定理。定义为有 k 对关系：P % ai = bi，其中 ai 两两之间互质。（而两两之间不互质就是把原来的关系式化为：P = bi (mod ai) → ai * x + bi = P，用拓展欧几里德求解同余方程组了。） 而 ai 两两互质时，可知道 a2*a3*...* ...
分类：
其他好文 时间：
2016-11-15 23:52:30
阅读次数：
314
题目：求在小于等于N的正整数中有多少个X满足：X mod a[0] = b[0], X mod a[1] = b[1], X mod a[2] = b[2], …, X mod a[i] = b[i], … (0 < a[i] <= 10)。 解法：先同上题一样用拓展欧几里德求出同余方程组的最后一个 ...
分类：
其他好文 时间：
2016-11-15 10:21:45
阅读次数：
131
Description 求解模线性方程组, \(m_i\) 不互质. Sol 扩展欧几里得+中国剩余定理. 首先两两合并跟上篇博文一样. 每次通解就是每次增加两个数的最小公倍数,这对取模任意一个数都是0. 伪代码如下 Code ...
分类：
其他好文 时间：
2016-11-12 09:33:59
阅读次数：
177
