四阶幻方

我写的四阶幻方出来了，运行结果正确，7040个的版本。虽然好像还有些问题，不过算了，能跑就可以了，呵呵。
不要输出的话我的p4 3.2G上耗时5s，算还过得去吧 ^_^
以下是程序：

cpp 代码

 

#include <iostream></iostream> #include <time.h></time.h>
using namespace std;      struct uni   {       int a,b,c,d;   };      uni can[2065];   long sum = 0;   int sum_all = 0;      void get_all_num()   {       for(int i = 1;i <= 16;++i)       {           for(int j = 1;j <= 16;++j)           {               for(int k = 1; k <= 16; ++k)               {                   for(int l = 1;l <= 16;++l)                   {                                      if(i != j && i != k && i != l && j != k && j != l && k != l)                           if(i + j + k + l == 34)                           {                               can[sum].a = i;                               can[sum].b = j;                               can[sum].c = k;                               can[sum].d = l;                               ++sum;                           }                   }               }           }       }   }      void p_rint(const int sq[])   {       cout << "i wIll gIvE yOu sOme coLoR 2CC! " << sum_all <     cout << sq[0] << " " << sq[1] << " " << sq[2] << " " << sq[3] << endl;       cout << sq[4] << " " << sq[5] << " " << sq[6] << " " << sq[7] << endl;       cout << sq[8] << " " << sq[9] << " " << sq[10] << " " << sq[11] << endl;       cout << sq[12] << " " << sq[13] << " " << sq[14] << " " << sq[15] << endl << endl;   }      inline bool caca(const int sq[])   {       if(sq[0] + sq[4] + sq[8] + sq[12] !=34)           return false;       if(sq[2] + sq[6] + sq[10] + sq[14] !=34)           return false;       return true;   }      inline bool detec_2(const int sq[])   {       int km2[17];       int i,m;       for(i = 0;i < 16;++i)       {           km2[i] = 0;       }          for(m = 0;m < 16;++m)       {           km2[sq[m]]++;           if(km2[sq[m]] > 1)               return false;       }          return true;   }   inline bool detec_1(const int c1[],const int c2[])   {       int km2[17];       for(int i = 0;i < 16;++i)       {           km2[i] = 0;       }          km2[c1[0] - 1]++;       km2[c1[1] - 1]++;       km2[c1[2] - 1]++;       km2[c1[3] - 1]++;              km2[c2[0] - 1]++;       km2[c2[1] - 1]++;       km2[c2[2] - 1]++;       km2[c2[3] - 1]++;          for(i = 0;i < 16;i++)       {           if(km2[i] > 1)               return false;       }          return true;   }      void fill_in_line(int c1[],int c2[])   {       int sq[16];          sq[0] = c1[0];       sq[5] = c1[1];       sq[10] = c1[2];       sq[15] = c1[3];          sq[3] = c2[0];       sq[6] = c2[1];       sq[9] = c2[2];       sq[12] = c2[3];                   for(int k = 1;k <= 16;k++)       {           sq[4] = k;                      for(int l = 1;l <= 16;++l)           {               sq[1] = l;                                  sq[7] = 34 - sq[4] - sq[5] - sq[6];               sq[11] = 34 - sq[3] - sq[7] - sq[15];               sq[8] = 34 - sq[9] - sq[10] - sq[11];               sq[13] = 34 - sq[1] - sq[5] - sq[9];               sq[14] = 34 - sq[12] - sq[13] - sq[15];               sq[2] = 34 - sq[0] - sq[1] - sq[3];                          if(sq[7] > 0 && sq[11] > 0 && sq[8] > 0 && sq[13] > 0 && sq[14] > 0 && sq[2] > 0)                   if(sq[7] <= 16 && sq[11] <= 16 && sq[8] <= 16 && sq[13] <= 16 && sq[14] <= 16 && sq[2] <= 16)                       if(detec_2(sq) && caca(sq))                       {                           sum_all++;                           p_rint(sq);                       }           }       }      }   void fill_in_eight()   {       int cse1[4],cse2[4];          for(int i = 0;i < sum;++i)       {           for(int j = 0;j < sum;++j)           {               if( i!=j )               {                   cse1[0] = can[i].a;                   cse1[1] = can[i].b;                   cse1[2] = can[i].c;                   cse1[3] = can[i].d;                      cse2[0] = can[j].a;                   cse2[1] = can[j].b;                   cse2[2] = can[j].c;                   cse2[3] = can[j].d;                                      if(detec_1(cse1,cse2))                   {                       fill_in_line(cse1,cse2);                   }               }                          }       }   }      int main()   {       time_t t;       t = clock();          get_all_num();       fill_in_eight();          cout << "time consumed: " << clock() - t << " ms"<< endl;       cin.get();       return 0;   }  

写得比较烂，咳咳咳

找个时间，来个ruby版的。
ps：最近发现牛人们都很臭屁，臭不可闻啊，唉，风气啊风气。做技术的人都很臭屁，果然啊果然，但愿我以后不会这样。