第四平方和定理，又稱為拉格朗日定理：每個正整數都可以表示為至多4個正整數的平方和。如果把 0包括進去，就正好可以表示為4個數的平方和。比如：5=0^2+0^2+1^2+2^2  7=1^2+1^2+1^2+2^2（^符號表示乘方的意思） 對於一個給定的正整數，可能存在多種平方和的表示法。要求你對4個數排序：0 <= a<=b<=c<=d 並對所有的可能表示法按a,b,c,d為聯合主鍵升序排列，最後輸出第一個表示法。  

【輸入】一個正整數N(N<5000000)  

【輸出】4個非負整數，按從小到大排序，中間用空格分開  

【輸入範例 1】5  

【輸出範例 1】0 0 1 2  

【輸入範例 2】773535

【輸出範例 2】1 1 267 838

程式程式碼：

#include<stdio.h>

#include<math.h>

#include<stdlib.h>                                                                  

int main()

{

    int N,a,b,c,d;       //N表示輸入的數，a,b,c,d分別為對應所找的數

    scanf("%d",&N);      //輸入N

    for(a=0;a<=pow(N,1.0/8);a++)    //縮小a的範圍，減少迴圈的次數,a從0開始，以逐次加一來迴圈

    {

       for(b=a;b<=pow(N,1.0/4);b++)    //縮小b的範圍，減少迴圈的次數，b>=a,b從a開始，以逐次加一來迴圈

       {             //因為a,b較小，所以a,b從最小的數開始進行篩選

           for(c=(int)(2*sqrt(N)-1);c>=b;c--) //縮小c的範圍，減少迴圈的次數，c>=b,因為c比較大，所以從最大的值開始逐步遞減來篩選

           {

              for(d=(int)sqrt(N);d>=c;d--)       //縮小d的範圍，減少迴圈的次數，d>=c,因為d比較大，所以從最大的值開始逐步遞減來篩選 

              {

                  if(a*a+b*b+c*c+d*d-N==0) //判斷是否與輸入的值相等

                  {

                     printf("%d %d %d%d\n",a,b,c,d);   //輸出

                     exit(0);             //輸出第一個後退出程式

                  }

              }

           }

        }

    }

    return 0;

}