csuoj 古怪的行列式

Description

这几天,子浩君潜心研究线性代数。 行列式的值定义如下:

DET 其中,τ(j1j2...jn)为排列j1j2...jn的逆序数。

子浩君很厉害的,但是头脑经常短路,所以他会按照行列式值的定义去计算,这个行列式子浩君也还是能算对的。但是,在计算的过程中,如果出现连续三行选取的元素为83(S),83(S),82(R)的话,子浩君会忍不住拿走它们:-D,然后这三个数的乘积将被视为1,而其它数值计算不变。那么在子浩君的计算下,最后得到的行列式的值会为多少呢?

Input

数据第一行为一个整数T(T<=50)。 接下来有T组测试数据,每组数据开始有一个整数n(2<=n<=8)。 接下来有n行数字,每行有n个数字,第ith行第jth个数字代表矩阵的第ith行第jth列的数字,保证每个数字在int范围内的非负整数。

Output

输出一个整数,保证在[-(2^63-1), 2^63-1]范围内,即使在子浩君计算过程中也是。

Sample Input

4

2
1 1
0 1

3
83 1 1
0 83 1
0 0 82

3
83 1 1
0 82 1
0 0 83

3
83 1 1
0 83 1
0 1 82

Sample Output

1
1
564898
-82

Hint

例如,当子浩君遇到a11 * a22 * a33 * a44 = 83 * 83 * 82 * 1,会计算成1 * 1 = 1,而83 * 82 * 83 * 1或者83 * 83 * 1 * 82则不会改变运算规则


思路:考虑数据量,暴力即可。比较繁琐的是按行列式的规则计算

#include<cstdio>

const int MAX = 10;

long A[MAX][MAX];   //行列式
int n;

bool mark[MAX];     //每列是否被选择,false为未选
long number[MAX];   //每行选择的列号
long sum = 0;

void choose(int, int);  //穷举行列式的所有选择情况
int inversion_number();    //计算逆序数

void choose(int row)
{
    if (row == 1)
        for (int i = 1; i <= n; i++)
            mark[i] = false;

    if (row == n)
    {
        for (int i = 1; i <= n; i++)
        {
            if (mark[i] == false)
            {
                number[row] = i;
                break;
            }
        }
        long res = 1;
        long temp[MAX];
        for (int i = 1; i <= n; i++)
            temp[i] = A[i][number[i]];

        for (int i = 1; i <= n - 2; i++)
        {
            if (temp[i] == 83 && temp[i + 1] == 83 && temp[i + 2] == 82)
            {
                temp[i] = 1;
                temp[i + 1] = 1;
                temp[i + 2] = 1;
            }
        }

        for (int i = 1; i <= n; i++)
            res *= temp[i];

        int j = inversion_number();
        if (j % 2 == 1)
            res = -res;
        sum += res;
        return;
    }

    bool use[MAX];
    for (int i = 1; i <= n; i++)
        use[i] = mark[i];

    for (int i = 1; i <= n - row + 1;)  //穷举的次数
    {
        for (int j = 1; j <= n; j++)    //每次都从1开始找
        {
            if (mark[j] == false && use[j] == false)
            {
                //选择该列
                mark[j] = true;
                number[row] = j;
                choose(row + 1);
                //恢复,考虑下一种情况
                mark[j] = false;
                use[j] = true;
                i++;
            }
        }
    }
}

int inversion_number()
{
    int cnt = 0;
    for (int i = 1; i <= n - 1; i++)
    {
        for (int j = i + 1; j <= n; j++)
            if (number[i] > number[j])
                cnt++;
    }
    return cnt;
}

int main()
{
    int T;
    scanf("%d", &T);

    while (T--)
    {
        scanf("%d", &n);

        sum = 0;

        for (int i = 1; i <= n; i++)
            for (int j = 1; j <= n; j++)
                scanf("%ld", &A[i][j]);

        choose(1);
        printf("%ld\n", sum);
    }

    return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值