2007-08-23

四阶幻方

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

  4. using namespace std;  
  5.   
  6. struct uni  
  7. {  
  8.     int a,b,c,d;  
  9. };  
  10.   
  11. uni can[2065];  
  12. long sum = 0;  
  13. int sum_all = 0;  
  14.   
  15. void get_all_num()  
  16. {  
  17.     for(int i = 1;i <= 16;++i)  
  18.     {  
  19.         for(int j = 1;j <= 16;++j)  
  20.         {  
  21.             for(int k = 1; k <= 16; ++k)  
  22.             {  
  23.                 for(int l = 1;l <= 16;++l)  
  24.                 {                 
  25.                     if(i != j && i != k && i != l && j != k && j != l && k != l)  
  26.                         if(i + j + k + l == 34)  
  27.                         {  
  28.                             can[sum].a = i;  
  29.                             can[sum].b = j;  
  30.                             can[sum].c = k;  
  31.                             can[sum].d = l;  
  32.                             ++sum;  
  33.                         }  
  34.                 }  
  35.             }  
  36.         }  
  37.     }  
  38. }  
  39.   
  40. void p_rint(const int sq[])  
  41. {  
  42.     cout << "i wIll gIvE yOu sOme coLoR 2CC! " << sum_all <
  43.     cout << sq[0] << " " << sq[1] << " " << sq[2] << " " << sq[3] << endl;  
  44.     cout << sq[4] << " " << sq[5] << " " << sq[6] << " " << sq[7] << endl;  
  45.     cout << sq[8] << " " << sq[9] << " " << sq[10] << " " << sq[11] << endl;  
  46.     cout << sq[12] << " " << sq[13] << " " << sq[14] << " " << sq[15] << endl << endl;  
  47. }  
  48.   
  49. inline bool caca(const int sq[])  
  50. {  
  51.     if(sq[0] + sq[4] + sq[8] + sq[12] !=34)  
  52.         return false;  
  53.     if(sq[2] + sq[6] + sq[10] + sq[14] !=34)  
  54.         return false;  
  55.     return true;  
  56. }  
  57.   
  58. inline bool detec_2(const int sq[])  
  59. {  
  60.     int km2[17];  
  61.     int i,m;  
  62.     for(i = 0;i < 16;++i)  
  63.     {  
  64.         km2[i] = 0;  
  65.     }  
  66.   
  67.     for(m = 0;m < 16;++m)  
  68.     {  
  69.         km2[sq[m]]++;  
  70.         if(km2[sq[m]] > 1)  
  71.             return false;  
  72.     }  
  73.   
  74.     return true;  
  75. }  
  76. inline bool detec_1(const int c1[],const int c2[])  
  77. {  
  78.     int km2[17];  
  79.     for(int i = 0;i < 16;++i)  
  80.     {  
  81.         km2[i] = 0;  
  82.     }  
  83.   
  84.     km2[c1[0] - 1]++;  
  85.     km2[c1[1] - 1]++;  
  86.     km2[c1[2] - 1]++;  
  87.     km2[c1[3] - 1]++;  
  88.       
  89.     km2[c2[0] - 1]++;  
  90.     km2[c2[1] - 1]++;  
  91.     km2[c2[2] - 1]++;  
  92.     km2[c2[3] - 1]++;  
  93.   
  94.     for(i = 0;i < 16;i++)  
  95.     {  
  96.         if(km2[i] > 1)  
  97.             return false;  
  98.     }  
  99.   
  100.     return true;  
  101. }  
  102.   
  103. void fill_in_line(int c1[],int c2[])  
  104. {  
  105.     int sq[16];  
  106.   
  107.     sq[0] = c1[0];  
  108.     sq[5] = c1[1];  
  109.     sq[10] = c1[2];  
  110.     sq[15] = c1[3];  
  111.   
  112.     sq[3] = c2[0];  
  113.     sq[6] = c2[1];  
  114.     sq[9] = c2[2];  
  115.     sq[12] = c2[3];           
  116.   
  117.     for(int k = 1;k <= 16;k++)  
  118.     {  
  119.         sq[4] = k;  
  120.           
  121.         for(int l = 1;l <= 16;++l)  
  122.         {  
  123.             sq[1] = l;  
  124.                   
  125.             sq[7] = 34 - sq[4] - sq[5] - sq[6];  
  126.             sq[11] = 34 - sq[3] - sq[7] - sq[15];  
  127.             sq[8] = 34 - sq[9] - sq[10] - sq[11];  
  128.             sq[13] = 34 - sq[1] - sq[5] - sq[9];  
  129.             sq[14] = 34 - sq[12] - sq[13] - sq[15];  
  130.             sq[2] = 34 - sq[0] - sq[1] - sq[3];  
  131.           
  132.             if(sq[7] > 0 && sq[11] > 0 && sq[8] > 0 && sq[13] > 0 && sq[14] > 0 && sq[2] > 0)  
  133.                 if(sq[7] <= 16 && sq[11] <= 16 && sq[8] <= 16 && sq[13] <= 16 && sq[14] <= 16 && sq[2] <= 16)  
  134.                     if(detec_2(sq) && caca(sq))  
  135.                     {  
  136.                         sum_all++;  
  137.                         p_rint(sq);  
  138.                     }  
  139.         }  
  140.     }  
  141.   
  142. }  
  143. void fill_in_eight()  
  144. {  
  145.     int cse1[4],cse2[4];  
  146.   
  147.     for(int i = 0;i < sum;++i)  
  148.     {  
  149.         for(int j = 0;j < sum;++j)  
  150.         {  
  151.             if( i!=j )  
  152.             {  
  153.                 cse1[0] = can[i].a;  
  154.                 cse1[1] = can[i].b;  
  155.                 cse1[2] = can[i].c;  
  156.                 cse1[3] = can[i].d;  
  157.   
  158.                 cse2[0] = can[j].a;  
  159.                 cse2[1] = can[j].b;  
  160.                 cse2[2] = can[j].c;  
  161.                 cse2[3] = can[j].d;  
  162.                   
  163.                 if(detec_1(cse1,cse2))  
  164.                 {  
  165.                     fill_in_line(cse1,cse2);  
  166.                 }  
  167.             }  
  168.               
  169.         }  
  170.     }  
  171. }  
  172.   
  173. int main()  
  174. {  
  175.     time_t t;  
  176.     t = clock();  
  177.   
  178.     get_all_num();  
  179.     fill_in_eight();  
  180.   
  181.     cout << "time consumed: " << clock() - t << " ms"<< endl;  
  182.     cin.get();  
  183.     return 0;  
  184. }  

写得比较烂，咳咳咳


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

您还没有登录，请登录后发表评论

t0uch
搜索本博客
最近加入圈子
存档
最新评论