[交流] Project Euler 48 欧拉工程 48 题已有5人参与

4楼: Originally posted by sudo at 2011-09-07 16:29:58:
因为
(a*b)%n == ((a%n)*(b%n))%n
故
[code]
#include <stdio.h>

int main()
{
    unsigned int i, j;
    unsigned long long sum=0, tmp;
    const unsigned long long k=10000000000ULL;
...

8楼: Originally posted by sudo at 2011-09-08 13:09:37:
+.+瞪了很久，有点牵强：
22  26  (20)  23  24  21  (23)  18  (27)

CODE:

22   26  (20)
23   24   21
(23) 18  (27)

正好每行每列的和都是68...

10楼: Originally posted by sudo at 2011-09-08 14:07:02:
不思考啦伤脑细胞...

不过这个“通项公式”是怎么定义的？是不是必须能写成f(n)的形式而且要求f(n)不能递归定义、不能为分段函数？或者还有什么别的限制条件（比如需不需要在n属于[0, +∞)上皆有定义）？
...

11楼: Originally posted by wangww2011 at 2011-09-08 16:17:27:
觉得出这种题的人最是无聊，这种题目有时候没有唯一解，只要你说的有道理，他就不能算错吧，譬如这个数列可以为

CODE:

22, 26, 23, 23, 24, 21, 25, 18,

理由是奇数项为等差数列，即22,23,24,25可以 ...

12楼: Originally posted by tieer at 2011-09-08 16:20:53:
22   26   （）   23   24    21  （）   
这几个数都减去18，
4 , 8，（）， 5, 6，3，（）

22   26   （）   23   24    21  （）  
这几个数的个位十位相加，
4 ，8 ，（），5 ，6 ，3， （）相同这样 ...

17楼: Originally posted by wangww2011 at 2011-09-08 19:04:51:
呵呵 你是没有仔细看啊
通项公式是

CODE:

a(2*n+1)=21+n
a(2*n)=28-5/2*n+(1/4)*(-1)^n+1/4

其中n取1,2,3,...

21楼: Originally posted by wangww2011 at 2011-09-08 21:55:35:
楼上的给出一个了，不用sin,cos也一样

a(2*n-1)=21+n
a(2*n)=28-5/2*n+(1/4)*(-1)^n+1/4
====>
[code]
a(n)=(1+(-1)^(n+1))/2*(21+(n+1)/2)+(1+(-1)^n)/2*(28-5/4*n-1/4*(-1)^(n/2)+1/4)
[/code ...

26楼: Originally posted by sudo at 2011-09-08 22:10:25:
其实多项式拟合就不复杂，就是求解痛苦而已，总而言之就是只要给出的数列项数有限，那么它就是万能的方法，而且对于空出来的位置还可以有无数多解

29楼: Originally posted by sudo at 2011-09-09 14:13:48:
这个好办

对于一个给出长度为n的序列，先待定最高项次数为n-1就好了嘛，这样必存在一个解。举个简答的例子

CODE:

n      1   2   3
a(n)   1   2   3

可以令a(n) = a2*n^2 + a1*n + a ...