C语言一天一练--22.august

gloomyakon

计算起来还蛮费劲儿的，所以只等它算出个部分结果。不知道有没有更快的计算方法，求教。

1. #include<stdio.h>
   
2. #include<stdlib.h>
   
3. #include<math.h>
   
4. /*
   
5. 1.12位以下的雷劈数，则超出int的表示范围了。使用long long才行
   
6. 2.如果是雷劈数则(a+b)^2=a*[log10(b)]+b，其中[]是向上取整。
   
7. 3.但是C语言里的log10等函数没有针对long long类型的所以要强制转换。*/
   
8. int main()
   
9. {
   
10. 
    
11.         long long sum;
    
12.         long long product;
    
13.         long long a,b;
    
14.         for (a=1;a<999999999999;a++)
    
15.         {
    
16.                 for(b=0;b<999999999999;b++)
    
17.                 {
    
18.                         product = (a+b)*(a+b);
    
19.                         if(log10((long double)product)>12)
    
20.                                 break;
    
21.                         sum=a*pow(10,floor(log10((long double)b)+1))+b;
    
22.                         if(product==sum)
    
23.                                 printf("a=%lld,b=%lld,\t%lld\n",a,b,product);
    
24.                 }
    
25.         }
    
26.         system("pause");
    
27.         return 0;
    
28. }

复制代码