POJ 1166 The Clocks 解题报告(高斯消元法 & 逆矩阵)

    题目大意:9种操作可以让不同的种转动90度,求最小上升的操作方式。

    解题报告:这题解法众多,也很有趣。可以BFS,DFS,9重循环暴搜也没问题= =。当然,为了学习还是用高斯消元法做的。

    Discuss也有人讨论了,4不是质数,求解过程中不能模4,不一定有解的问题。按照我的理解,题目既然说了有唯一解,就不用考虑这个问题了。

    另外,寻找当前列的对应行时不能选绝对值最大的,会WA。具体原因不详……代码如下:

#include 
#include 
#include 
#include 
using namespace std;

int a[10][10];
int x[10];
bool free_x[10];
int sum;

int Gauss(int equ,int var)
{
    memset(x,0,sizeof(x));
    memset(free_x,1,sizeof(free_x));

    int k=0,col=0;
    while(k=0;i--)
    {
        int temp=a[i][var];
        for(int j=i+1;j0?" ":"\n");
        }
}

    当然,也可以事先求出转动矩阵的逆矩阵,直接求解= =代码如下:

// 逆矩阵
#include 
#include 
using namespace std;

int a[9][9]={
{3,2,3,2,2,1,3,1,0,},
{3,3,3,3,3,3,2,0,2,},
{3,2,3,1,2,2,0,1,3,},
{3,3,2,3,3,0,3,3,2,},
{3,2,3,2,1,2,3,2,3,},
{2,3,3,0,3,3,2,3,3,},
{3,1,0,2,2,1,3,2,3,},
{2,0,2,3,3,3,3,3,3,},
{0,1,3,1,2,2,3,2,3,}};

int x[9];
int res[9];

int main()
{
    for(int i=0;i<9;i++)
    {
        scanf("%d",x+i);
        x[i]=(4-x[i])%4;
    }

    for(int i=0;i<9;i++)
        for(int j=0;j<9;j++)
            res[i]+=a[i][j]*x[j];

    for(int i=0;i<9;i++) while(res[i]%4 && res[i]--)
        printf("%d ",i+1);
    puts("");
}


你可能感兴趣的:(数学,ACM)