Project Euler 21

这个问题虽然放了挺久……但是似乎是挺简单的~

我曾经想过有没有什么高效算法，但最后没有想出来……呵呵……

看了一下forum，似乎做法都差不多……硬算~

上一下我的code

#include <stdio.h>

int
main()
{
    int i = 0, now, j;
    int sum = 0;
    int f[10000];
    for (i=1; i<10000; i++)
    {
	now = i;
	f[i] = 1;
	for (j=now/2; j>1; j--)
	    if (now%j == 0)
		f[i] +=j;
    }
    for (i=1; i<10000; i++)
	if ((f[i] < 10000) && (f[f[i]] == i) && (f[i]!=i))
	{
	    sum += i;
	    printf("%d %d\n",i, f[i]);
	}
    printf("%d", sum);
    return 0;
}

别嫌我的代码丑……我在emacs里面显示的很好的……擦……

一出来就不行了……

再来就是jason的~

不过有个问题，我觉得我的代码使用了空间换时间的方法

用 f 储存，这样减少了检查这部分。

谁想到……jason大牛竟然对着自己这个不怎么快速的代码优化了半天……还用来比较python，cython，c……

上一下，人家写的还是好看……

#include <stdio.h>
#include <stdlib.h>
 
unsigned int sum_divisors(unsigned int n) {
    unsigned int s = 0, i = 1;
    while(i < n) {
        if(n % i == 0) s = s + i;
        i++;
    }
    return s;
}
 
unsigned int sum_amicable_pairs(unsigned int low, unsigned int high) {
    unsigned int a = low, b, sum = 0;
    while(a <= high) {
        b = sum_divisors(a);
        if(b > a && sum_divisors(b) == a) sum = sum + a + b;
        a++;
    }
    return sum;
}
 
int main(int argc, char **argv) {
    unsigned int ans = sum_amicable_pairs(1,10000);
    printf("result: %d\n",ans);
    return 0;
}

看这个main……就觉得像python上来一句……

if '__name__' == '__main__':
	main

显得就高端了……

后面说个闲话，今天CE1004，great ideas in computing 考final quiz……我竟然迟到了……啊……一天了……都没释怀……

还有就是math lab的 scilab 实在是不怎么好用……不喜欢，不科学！

直接用python多好~跟1003统一一下，还能练习编程……当然我只是吐个槽~呵呵……