http://poj.org/problem?id=2891 解线性模方程组。 比较坑爹,数据比较大,很容易溢出。 1 Program poj2891; 2 3 var m:int64; 4 5 a,r:array[1..30000000]of int64; 6 7 ans,x,y,lcm:int6
分类:
其他好文 时间:
2016-03-02 21:43:25
阅读次数:
149
一、线性最小二乘法 1、基本思路 令,其r(x)是事先选定的一组线性无关的函数。ak是待定系数。然后拟合的准则就是使得yi与f(xi)的距离的平方和最小,称之为最小二乘准则 2、系数的确定 ,要使距离的平方和最小,那只要取得,使得取到极值,就可以解除待定系数ak,记 然后线性方程组为,所以当R列满秩
分类:
编程语言 时间:
2016-02-29 00:17:27
阅读次数:
263
线性代数,原理不讲了。。。 1 /* 用于求整数解得方程组. */ 2 3 #include <iostream> 4 #include <string> 5 #include <cmath> 6 using namespace std; 7 8 const int maxn = 105; 9 10
分类:
其他好文 时间:
2016-02-17 18:55:21
阅读次数:
185
https://en.wikipedia.org/wiki/Gaussian_elimination Gauss–Seidel method - Wikipedia, the free encyclopedia.html
分类:
其他好文 时间:
2016-02-11 11:09:16
阅读次数:
169
最著名的十大公式 No.10 圆的周长公式(The Length of the Circumference of a Circle) No.9 傅立叶变换(The Fourier Transform) No.8 德布罗意方程组(The de Broglie Relations) No.7 1+1=2
分类:
其他好文 时间:
2016-02-05 19:13:37
阅读次数:
317
http://acm.hdu.edu.cn/showproblem.php?pid=1573 求小于等于N的正整数中有多少个X满足: X mod a0 = b0 X mod a1 = b1 …… X mod ai = bi (0<ai<=10) 输入:第一行为一个正整数T,表示有T组测试数据。每组测
分类:
其他好文 时间:
2016-02-02 21:29:35
阅读次数:
223
http://acm.hdu.edu.cn/showproblem.php?pid=3579 题解:同余方程组的裸题。注意输出是最小的正整数,不包括0。 #include<cstdio> #include<cstdlib> #include<cstring> #include<iostream> u
分类:
其他好文 时间:
2016-02-02 21:29:15
阅读次数:
169
同余方程组 例题1:pku2891Strange Way to Express Integers 中国剩余定理求的同余方程组mod 的数是两两互素的。然而本题(一般情况,也包括两两互素的情况,所以中国剩余定理成为了“时代的眼泪”)mod的数可能不是互素,所以要转换一下再求。 P=b1(mod a1)
分类:
其他好文 时间:
2016-02-02 14:09:54
阅读次数:
136
今天接到scy的压缩包,开始做数论专题。那今天就总结一下拓展欧几里得求解不定方程和同余方程组。 首先我们复习一下欧几里得算法: 1 int gcd(int a,int b){ 2 if(b==0) return a; 3 return gcd(b,a%b);4 } 拓展欧几里得算法: 推导过程: 给
分类:
其他好文 时间:
2016-02-01 23:51:53
阅读次数:
281
http://poj.org/problem?id=2891题意:求解一个数x使得 x%8 = 7,x%11 = 9; 若x存在,输出最小整数解。否则输出-1;ps:思路:这不是简单的中国剩余定理问题,由于输入的ai不一定两两互质,而中国剩余定理的条件是除数两两互质。 这是一般的模线性方程组...
分类:
其他好文 时间:
2016-01-26 15:19:12
阅读次数:
221