如何求出1-n以内勾股数之和

wayne

30# zeroieme 这个改得好， 因为a-b是奇数，所以a取 b+1,b+3,b+5,........ 仔细想了一下，如果a取步长2的话，那么if语句里面的其实就可以用gcd(a,b)==1 了

1. #include <stdio.h>
2. typedef unsigned long int DW;
3. inline DW gcd(DW m,DW n){DW t; while(m!=0) t=m,m=n%m,n=t;return n;}
4. DW Test(DW n){
5. DW t,s,ii,cnt=1,sum=0,c=1;
6. DW a,b;
7. for(b=1;b*b<=n/2;b++)
8. for(a=b+1;t=a*a+b*b,t<=n;a+=2){
9. if(gcd(a,b)==1){
10. s=n/t;
11. // for(ii=1;ii<=s;ii++)
12. // printf("%5u.%-5u\t%u,%u,%u\n",c++,cnt,ii*(a*a-b*b),2*ii*a*b,ii*t);
13. // cnt++;
14. sum+=s*(s+1)*a*(a+b);
15. }
16. }
17. return sum;
18. }
20. int main(int c ,char **v)
21. { if (c==2) {
22. printf("b: total of triples less than %s :\t%lu\n",v[1],Test(atoi(v[1])));
23. }
24. return 0;
25. }

复制代码

〇〇

看来用好算法很重要

wayne

32# 〇〇 在装有vxworks 的板子上运行了上面的程序，晒晒截图：

plp626

刚才发现〇〇 和wayne 互为马甲。。。

wayne

34# plp626

gxqcn

刚才发现〇〇 和wayne 互为马甲。。。 plp626 发表于 2011-5-12 18:28

这个似乎没有依据。 倒是感觉两人的行文差异很大，不会是同一人。

wayne

36# gxqcn 这个话题很有趣，大家可以深度挖掘一下

G-Spider

按本论坛的潜规则。头像不太对劲，应该不是马甲。

liangbch

两个人的IP不同

plp626

我是看它的百度空间猜测的。。 那两个圈圈。。。还有后面哪个date -d ... 很面熟。。。 http://hi.baidu.com/%B0%D7%C0%CB/home