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则不会改变运算规则
一开始想着用dfs做,结果一直wa(逆序数定义记不清了我也很绝望)
后来结果大佬指点,了解到一种新的操作 叫next_permutation 详细请看 这个大佬的博客 http://www.cnblogs.com/kuangbin/archive/2012/03/30/2424482.html
总之就是可以持续生成新的排列组合,而我们的行列式本来就是一层层取的,所以对应起来十分方便
代码:
#include<iostream>
#include<algorithm>
using namespace std;
long long nmap[10][10]; //nmap[行][列]
int posi[10]; //代表哪一列
int N;
void init()
{
for(int i=1;i<=N;i++)
for(int j=1;j<=N;j++)
cin>>nmap[i][j];
for(int i=1;i<=N;i++)
posi[i]=i; //因为要得到所有的排列组合,所以一开始得是排序好了的
}
int main()
{
int T;
long long sum;
cin>>T;
while(T--)
{
cin>>N;
init();
sum=0;
do
{
long long ans=1;
int ni=0;int po=1;
for(int i=1;i<=N;i++)
for(int j=i+1;j<=N;j++)
if(posi[i]>posi[j]) //获得逆序数
ni++;
if(ni%2) //逆序数不为偶,那整个行列式就是负的
po=-1;
for(int i=1;i<=N;i++)
{
if(nmap[i][posi[i]]==82&&nmap[i-1][posi[i-1]]==83&&nmap[i-2][posi[i-2]]==83)
ans/=6889; //这里检查一下前两个数是否为条件限制
else
ans*=nmap[i][posi[i]];
}
ans*=po;
sum+=ans;
}while(next_permutation(posi+1,posi+N+1)); //愉快的获得下一组排列
cout<<sum<<endl;
}
return 0;
}