数据结构之魔方程序

魔方构造问题


一个魔方(magic square)就是由1到n平方的整数构成的n*n矩阵,其中每行每列以及主对角线上的数字之和都相等。
当n为奇数时,Coxeter法则可以产生魔方:
把1放入第一行最中间的方格中。向左上方移动,并按照数字的递增顺序,把数字填入空方格。如果移出了魔方(即越过了魔方边界),则进入魔方对边的对应放歌。继续填写方格。如果一个方格已被填入数字,则向下继续填写。
魔方程序涉及到的Coxeter法则属于群论。
魔方程序迭代写法:

/*魔方程序-迭代*/
#include <stdio.h>
#include <stdlib.h>
#define MAX_SIZE 15 /* maximum size of square */
void main(){
    /* construct a magic square, iteratively */
    static int square[MAX_SIZE][MAX_SIZE];
    int i, j, row, column; /* indices */
    int count; /* counter */
    int size;

    printf("Enter the size of square: ");
    scanf("%d", &size);
    /* check for input errors */
    if(size < 1 || size > MAX_SIZE + 1){
        fprintf(stderr, "Error! Size is out of range\n");
        exit(1);
    }
    if(!(size % 2)){
        fprintf(stderr, "Error! Size is even\n");
        exit(1);
    }
    for(i = 0; i < size; i++)
        for(j = 0; j < size; j++)
            square[i][j] = 0;
    square[0][(size-1) / 2] = 1; /* middle of first row */
    /* i and j are current position */
    i = 0;
    j = (size-1) / 2;
    for(count = 2; count <= size*size; count++){
        row
    }
}

魔方程序递归写法:

/*魔方程序-递归*/
#include <stdio.h>
#include <stdlib.h>
#define MAX_SIZE 15 /* maximum size of square */

void checkInput(int size){
    /* check for input errors */
    if(size < 1 || size > MAX_SIZE + 1){
        fprintf(stderr, "Error! Size is out of range\n");
        exit(1);
    }
    if(!(size % 2)){
        fprintf(stderr, "Error! Size is even\n");
        exit(1);
    }
}
void initSquare(int size, int sq[][MAX_SIZE]){
    int i,j;
    for(i = 0; i < size; i++)
        for(j = 0; j < size; j++)
            sq[i][j] = 0;
}
void magicSquare(int i, int j, int n, int size, int sq[][MAX_SIZE]){
    if(n == size*size){
        sq[i][j] = n;
    }
    else{   
        sq[i][j] = n;
        i = (i-1 < 0) ? (size - 1):(i - 1);
        j = (j-1 < 0) ? (size - 1):(j - 1);
        if(sq[i][j]){
            i = (i+2) % size;
            j = (j+1) % size;
        }
        magicSquare(i,j,n+1,size,sq);
    }
}
void output(int size, int sq[][MAX_SIZE]){
    int i,j;
    printf(" Magic Square of size %d: \n\n", size);
    for(i = 0; i < size; i++){
        for(j = 0; j < size; j++)
            printf("%5d", sq[i][j]);
        printf("\n");
    }
    printf("\n\n");
}
void main(){
    static int square[MAX_SIZE][MAX_SIZE];
    int i, j;
    int size;
    printf("Enter the size of square: ");
    scanf("%d", &size);
    checkInput(size);
    /* init the magic square */
    initSquare(size, square);
    i = 0;
    j = (size-1) / 2;
    magicSquare(i,j,1,size,square);
    /* output the magic square */
    output(size, square);
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值