找回密码
 欢迎注册
查看: 11202|回复: 9

[擂台] n阶幻方全解

[复制链接]
发表于 2009-1-5 10:12:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
求指定的n阶幻方全部不同解。(n<100)
这个问题上次在火车上想了一下,笔算了一下3阶应该是8种情况,大概我想到了两种适合编程的思路(不知道对不对),
一、解线性方程组+搜索。
二、将1-n*n的数,分成每组和都为幻和的n组,找到全部的分法,然后将每组内每行的数进行位置调整,得到全部解
不过个人觉得这两种方法运算量太大,想了解下各位的想法。
这方面GXQCN应该是“专家”,不知道有没有好的方法?
==补充=========
不一定要求答案,提供一点理论上可行的思路就行。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-1-5 10:15:53 | 显示全部楼层
5以上很难
4的我做过
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-1-5 10:19:19 | 显示全部楼层
http://qzc.zgz.cn/X-huanfang7.htm
================================
幻方的数量研究简介

3阶幻方只有一个。由对称与旋转,有8种形式。

4阶幻方因为幻和是34,可以用穷举法获得所有的幻方。早在1693年,弗兰尼克尔(Bernard Frenicle de Bessy)就已经知道4阶幻方共有880个,经对称可得7040个不同的形式。 本站给出了这本人编程得到的这全部4阶幻方880个, 同时还有一个Flash软件产生这些幻方。

4 阶泛对角幻方,共有384个,有趣的是,由任一个4阶泛对角幻方出发,经一定的变换,可以得出所有这384个泛对角幻方。近年来,一些学者讨论这个问题。参见徐桂芳等《纯幻方的构造原理与方法》(西安交大出版社,1993),潘凤雏《构造所有4阶泛对角线幻方的统一公式》(合肥工大,大学数学,2005.3)。

5阶幻方的数量,人们原先以为可以类似3,4阶得出,但经过数百年仍无结果。产生5阶幻方有各种方法,除前面介绍的两个方法,还有一些,如拉伊尔法,尔仅用此法就可得出57600种不同的幻方。有人估计总数在130000000(一亿三千万)以上,但很多人不相信。1973年,有人用计算机得出5阶幻方的总数是275 305 224个。

至于6阶和更高阶数的幻方数量,现在仍是研究的课题,并已经有一些了不起的结果
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-1-5 10:35:18 | 显示全部楼层
幻方更多的是研究其美学性质,比如和谐美、奇异美等。
而对于其数目的有限性一般不做研究,因为它是个天文数字,对于当前的计算能力等同于无限。
所以人们更多的研究如何快速构造、判定存在与否等问题,很少去计数统计,因为尚没有好的理论支持。

我现在可以比较随心所欲地构造幻方,
如果不要求数字的连续性,还可以定制非常美妙的幻方,
比如将一个电话号码或身份证数字分成几段后隐藏在幻方的某行的连续格中。

唉,不碰幻方已多年。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2009-1-5 10:50:37 | 显示全部楼层

啊,比我想像的多很多啊!

我也觉得比较多,没想到有这么多
怪不得网上找不到答案呢!
不好意思出了个这样的题目,那换成大家想点可能可以解决的想法,不一定要解出来。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2009-1-5 11:12:52 | 显示全部楼层

回复 3# 无心人 的帖子

你的那个链接我看了,
跟我以前看过的一些构造方法还不一样呢。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-9-23 11:56:59 | 显示全部楼层

RE: n阶幻方全解

SnakeD.rar (25.45 KB, 下载次数: 13) 贴一个显示5×5幻方的旋转型程序!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-9-24 20:27:17 | 显示全部楼层
求指定的n阶幻方全部不同解。(n
winxos 发表于 2009-1-5 10:12

据说当年计算全部275305224个5×5幻方花了100个小时,几十年过去了,计算机的发展突飞猛进,而算法的研究却进展不大。按现在计算机的速度,我的算法也要用60多个小时,不过因为能同时运行3个程序,所以也能在24小时以内搞定(配置为酷睿2 四核)。
我是按4角排列组合的,角上4个数共有24种排列,其中8个1组是同构的,实际上只有3种,再加上逆幻方,搜索的范围是总数的16分之1。算法的关键在于剪枝,根据幻和分组的意义不大,因为求出全部满足幻和的组合也不需要1秒钟,所以预先搞定能节约的时间实在有限。
3种幻方分为:交叉型X,12与34交叉;平行型Z,12与34平行;旋转型D,12与34反向。
X:89661808,Z:96617928,D:89025488。合计:275305224。
楼上给出的是旋转型D的显示程序,统计时是关闭显示的!
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2011-9-24 20:31:36 | 显示全部楼层
1234指的是由小到大的四个数,
|        X      |         Z        |           D           |
|           1   3       |      1    2     |        1      2        |
|           4   2       |     3    4     |       4      3        |
简称XZD法。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2012-3-12 21:02:26 | 显示全部楼层
以前看了个生成幻方的程序,不知道谁写的,在给大家看一下(c++的)

#include<iostream.h>
#include <stdio.h>
#include<iomanip.h>
void jishu(int n)   //任意奇数阶幻方制作的函数(阶梯法)
{
        int p=n-1,q=0,s=1,i,j,t=n-1;
        int b[50][50];
        int a[100][100];
        for(i=0;i<2*n-1;i++)       //给数组a赋初值
                for(j=0;j<2*n-1;j++)
                        a[i][j]=0;
        while(p<2*n-1)            //从a的中间起向下运动一直到最后
        {
                {
                        for(i=p;i>=p+1-n;i--)   //根据行的变换,列数每次右移五格
                                for(j=q;j<q+n;j++)
                                        if(i+j==t)
                                        {
                                                a[i][j]=s;  //每次赋值都加1
                                                s++;
                                        }
                }
                        t+=2;
                        q++;
                        p++;
        }

        for(i=0;i<2*n-1;i++)   //把中间框以外的数字按规律填到中间的框内
                for(j=0;j<2*n-1;j++)
                {
                        if(i<(n-1)/2&&a[i][j]!=0)   //格子上方的数填入数组中
                        {
                                a[i+n][j]=a[i][j];
                        }
                                else if(i>=(3*n-1)/2&&a[i][j]!=0)   //格子下方的数填入数组中
                                {
                                        a[i-n][j]=a[i][j];
                                }
                                        else if(j<(n-1)/2&&a[i][j]!=0)   //格子左方的数填入数组中
                                        {
                                                a[i][j+n]=a[i][j];
                                        }
                                                else if(j>=(3*n-1)/2&&a[i][j]!=0) //格子右方的数填入数组中
                                                {
                                                        a[i][j-n]=a[i][j];
                                                }
                }
        for(i=0;i<n;i++)     //用一个新数组把a的中间提取出来
                for(j=0;j<n;j++)
                        {
                                b[i][j]=a[i+(n-1)/2][j+(n-1)/2];
                        }
        for(i=0;i<n;i++)    //换行输出,每行5个
        {
                for(j=0;j<n;j++)
                        cout<<setw(3)<<b[i][j]<<" ";
                cout<<endl;
        }
}
void f1(int n)    //设计的一个通用双偶函数
{
        int a[100][100];
   int i,j;
        int x,k=1,m=3,t;
t=(n-1)/4;
        if(t==0)   //即n=4
        {
                for(i=0;i<=m;i++)         //给数组依次赋初值1,2,3…
                        for(j=0;j<=m;j++)     
                        {
                                a[i][j]=k;k++;
                        }
                        for(i=1;i<=2;i++)   //用走日字的方法交换数组中的元素
                        {
                         x=a[0][i];a[0][i]=a[m][m-i];a[m][m-i]=x;
             x=a[i][0];a[i][0]=a[m-i][m];a[m-i][m]=x;
                        }
                        for(i=0;i<=m;i++)    //输出数组
                        {
                                for(j=0;j<=m;j++)
                     cout<<a[i][j]<<" ";
                                cout<<endl;
                        }goto loop;      //跳到loop
        }

for(i=0;i<8;i++)        //给数组赋初值0
                for(j=0;j<8;j++)
                        a[i][j]=0;
        if(n%4==0&&n/4>1)
        {
        for(i=0;i<2;i++)     //取左上角四个格子把主对角线上的元素赋值为1
                for(j=0;j<2;j++)
                {
if(i==j)
                                a[i][j]=1;
                }
        }
                while(t>0)   
        {
        for(i=0;i<2;i++)    //水平复制刚才的四格   
        for(j=2;j<4;j++)
                        a[i][j]=a[i][j-2];
        for(i=2;i<4;i++)       //将上述八格水平翻转
for(j=0;j<4;j++)
                        a[i][j]=a[3-i][j];
             


for(i=0;i<n;i++)           //水平对称翻转  
for(j=n/2;j<n;j++)
                        a[i][j]=a[i][n-1-j];


for(i=n/2;i<n;i++)   //垂直对称翻转            for(j=0;j<n;j++)
                        a[i][j]=a[n-1-i][j];
                t--;
        }
        for(i=0;i<n;i++)     //从左上角起,给值为0的元素重新赋值,值为n*i+j+1
                for(j=0;j<n;j++)
                        if(a[i][j]==0)
                                a[i][j]=n*i+j+1;
        for(i=n-1;i>=0;i--)   //从右下角起,给值为1的元素重新赋值,值为n*n-n*i-j
                for(j=n-1;j>=0;j--)
                        if(a[i][j]==1)
                                a[i][j]=n*n-n*i-j;
                for(i=0;i<n;i++)   //输出数组
{
                        for(j=0;j<n;j++)
                                cout<<setw(3)<<a[i][j]<<" ";
                                cout<<endl;
                        }
loop:;
}
void danou(int n)    //设计的一个通用单偶函数
{
        int p=n-1,q=0,s=1,i,j,t=n-1,v=n/2;
        int a[50][50],b[50][50],c[50][50],d[50][50],e[100][100];
        int m=(n-2)/4;
        for(i=0;i<n/2;i++)   //给数组赋初值
                for(j=0;j<n/2;j++)
                        a[i][j]=0;
                int x,y;
        int max;
        int mid;
        max=v*v;
        mid=v/2;
        x=v-1;
        y=mid+1;
        for(i=0;i<n/2;i++)   //将a做一次奇数阶幻方变换
                for(j=0;j<n/2;j++)
                        a[i][j]=0;     //设置初始值是0
                a[0][mid]=1;
                for(i=2;i<=max;i++)
                {
                        if(x<0&&y>n/2-1)
                        {
                                x=0;y=n/2-1;x=x+1;     //在最右上角超出
                        }
                        if(x<0)
                        {
                                x=n/2-1;      //在上方超出
                        }
                        if(y>n/2-1)     //在右方超出
                        {
                                y=0;
                        }
                        if(a[x][y]!=0)
                        {
                                x=x+2;y=y-1;     //在右上角有数字
                        }
                        a[x][y]=i;    //正常附值
                        x=x-1;   //正常附值
                        y=y+1;   //正常附值
                }

        for(i=0;i<n/2;i++)    //给b数组赋值,在a的对应数上加2*v*v
                for(j=0;j<n/2;j++)
                        b[i][j]=a[i][j]+2*v*v;
        for(i=0;i<n/2;i++)    //给c数组赋值,在a的对应数上加v*v
                for(j=0;j<n/2;j++)
                        c[i][j]=a[i][j]+v*v;
        for(i=0;i<n/2;i++)      //给d数组赋值,在a的对应数上加3*v*v
                for(j=0;j<n/2;j++)
                        d[i][j]=a[i][j]+3*v*v;
                        int temp;   //定义交换变量
                for(i=0;i<=n/2-1;i++)     //将a和d的前面进行交换
                        for(j=0;j<m;j++)
                        {
                                temp=a[i][j];
                                a[i][j]=d[i][j]        ;
                                d[i][j]=temp;
                        }
                            temp=a[(n-2)/4][0];     //a和d中间一行第一列的数再交换过来,即不交换
                            a[(n-2)/4][0]=d[(n-2)/4][0];       
d[(n-2)/4][0]=temp;

                        temp=a[(n-2)/4][m];       //a和d中间一行第二列的数交换过来
                        a[(n-2)/4][m]=d[(n-2)/4][m];
                        d[(n-2)/4][m]=temp;

                        for(i=0;i<=n/2-1;i++)    //将c和b的最后列进行交换
                                for(j=n/2-1;j>(n-2)/4+1;j--)
                                {
                                        temp=b[i][j];
                                        b[i][j]=c[i][j];
                                        c[i][j]=temp;
                                }
        for(i=0;i<n/2;i++)     //用大数组e【100】【100】把四个小数组拼接起来
                for(j=0;j<n/2;j++)
                        e[i][j]=a[i][j];
        for(i=0;i<n;i++)
                for(j=n/2;j<n;j++)
                        e[i][j]=b[i][j-n/2];
        for(i=n/2;i<n;i++)
                for(j=0;j<n/2;j++)
                        e[i][j]=d[i-n/2][j];
        for(i=n/2;i<n;i++)
                for(j=n/2;j<n;j++)
                        e[i][j]=c[i-n/2][j-n/2];
                for(i=0;i<n;i++)   //输出数组
{
                        for(j=0;j<n;j++)
                                cout<<setw(3)<<e[i][j]<<" ";
                        cout<<endl;
                }
}


void main()  //主函数
{    bool f=true;
        while(f)
{
        int n;                      //界面设计
        cout<<endl;
        cout<<setw(48)<<"       幻方制作"<<endl;
        cout<<endl;
        cout<<setw(50)<<"============="<<endl;
        cout<<endl;
        cout<<setw(50)<<"      1.奇数阶幻方"<<endl;
        cout<<setw(50)<<"      2.偶数阶幻方"<<endl;
        cout<<endl;
        cout<<setw(56)<<"请选择(1或2,0:退出)"<<endl;
        int k;
        cin>>k;
        //int k;
//        cin>>k;
        switch(k)   
        {
        case 1:
                cout<<"请输入奇数n(2<n<=10)"<<endl;
                cin>>n;
                if(n<3||n>=10||n%2==0)
                {
                        do{
                            cout<<"对不起,您的输入有误,请重新来过!"<<endl;
                           cin>>n;
                   }while(n>=3&&n<10&&n%2!=0);
                }
                else
                        jishu(n);
                break;
        case 2:
                        cout<<"请输入需要的偶数n(2<n<=10)"<<endl;
                cin>>n;
                if(n<3||n>10||n%2!=0)
                {
                        do{
                cout<<"对不起,您的输入有误,请重新来过!"<<endl;
                           cin>>n;
                        }while(n>3&&n<=10&&n%2==0);
                }
                else{
                       if(n%4==0)
                                 f1(n);
                        else if(n%2==0&&n%4!=0)
                                  danou(n);
                }
                break;
        case 0:f=false;break;
        default:
                cout<<"查无此号"<<endl;
}
}
}

幻方制作.txt

6.04 KB, 下载次数: 2, 下载积分: 金币 -1 枚, 经验 1 点, 下载 1 次

毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-3-28 17:36 , Processed in 0.049357 second(s), 19 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表