Description
这几天,子浩君潜心研究线性代数。 行列式的值定义如下:
其中,τ(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;
}