百科问答小站 logo
百科问答小站 font logo



x^4+y^4+z^4+w^4=a^4有正整数解吗? 第1页

  

user avatar   yu-guang-ting-92 网友的相关建议: 
      

以下是一些 的正整数解

利用多线程库在28核56线程CPU上跑大概几分钟就能出第一组解,后面的几组是在64核128线程CPU上跑了8小时找出来的。

       #include<math.h> #include<stdio.h> #include<stdlib.h> #include<threads.h>  const unsigned long NUM_THREADS = 1024; void solution(void*);  int main(){         thrd_t threads[NUM_THREADS];         for(unsigned long t=4; t<NUM_THREADS; t++){                 thrd_create(&threads[t], (thrd_start_t)solution, (void*)t);         }         thrd_exit(EXIT_SUCCESS);         return EXIT_SUCCESS; }  void solution(void* p){         long s = (long)p;         printf("%ld: pthread ID - %lu
", s, thrd_current());         for(long x=1; x<s-3; ++x){                 for(long y=x; x+y<s-2; ++y){                         for(long z=y; x+y+z<=s-z; ++z){                                 long w = s-x-y-z;                                 double d = pow(x,4)+pow(y,4)+pow(z,4)+pow(w,4);                                 if(abs(pow(round(sqrt(sqrt(d))), 4)-d)<1e-3){                                         s = (long)round(sqrt(sqrt(d)));                                         printf("%ld:	%ld	%ld	%ld	%ld
", s,x,y,z,w);                                 }                         }                 }         } }     

编译命令(文件名为main.c)

       cc main.c -lm -pthread -o main     




  

相关话题

  数学系的鄙视链是什么? 
  由AB=BA=O可以得出什么结论? 
  这两个级数该怎么解答? 
  有哪些看似简单其实非常精妙的代码? 
  减一个负数,为什么是加这个数的相反数呢? 
  一个一般的二次型等于0,这个方程应该如何求通解? 
  有哪些经典的反直觉数学结论? 
  高斯的博士论文是不是太简单了? 
  数学大牛是怎么看待悖论和无穷的? 
  如何评价第12届全国大学生数学竞赛初赛(数学类A组)试题? 

前一个讨论
量子力学老师提出了一个问题:为什么 Schrödinger 方程里有虚数 i ?
下一个讨论
据报道称“浏览器内核有上千万行代码”,浏览器内核真的很复杂吗?





© 2024-05-05 - tinynew.org. All Rights Reserved.
© 2024-05-05 - tinynew.org. 保留所有权利