如何求出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. }
    
19. 
    
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