CSU1979 古怪的行列式

E(1979): 古怪的行列式

Submit Page    Summary    Time Limit: 1 Sec     Memory Limit: 128 Mb     Submitted: 41     Solved: 30    

---

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则不会改变运算规则

 

这道题两个知识点，首先逆序数是指在一行数字中，当数字的大小和排列的顺序相反的所有情况的总数叫逆序数（例：4213的逆序数为四  42 41 43 21）。第二，行列式值的计算。剩下的就是模拟了，最多矩阵只有8，所以不会超时。

 

#include<map>
#include<algorithm>
#include<cstdio>
#include <iostream>
#include<cmath>
#include<string.h>
#include<cstring>
using namespace std;
typedef long long ll;
const double PI = acos(-1.0);
const int INF = 1000000000;
const int maxn = 10+2;
int T=0,n,m;

int b[maxn];
ll a[maxn][maxn];

int Caculate(int x[]){
     int k=0;
     for(int i = 1 ; i<=n ; i++)
     {
         for(int j = i+1 ; j<=n ; j++)
         {
              if(x[j]<x[i])  k++;
         }
     }
     return k;
}

int main(){
     //freopen("in.txt","r",stdin);
     ll ans;
     scanf("%d",&T);
     while(T--){
          scanf("%d",&n);
          for(int i = 1 ; i<=n ; i++)
          {
              b[i] = i;
              for(int j = 1 ; j<=n ; j++)
              {
                    scanf("%lld",&a[i][j]);
              }
          }
          ans = 0;
          do{
               int r = Caculate(b);
               //cout<<"r="<<r<<"\n";
               ll temp = 1;
               for(int i = 1 ; i<=n ; i++)
               {
                    if(i+2<=n&&a[i][b[i]]==83&&a[i+1][b[i+1]]==83&&a[i+2][b[i+2]]==82)
                        i = i+2;
                    else
                      temp = temp * a[i][b[i]];
                      //cout<<"temp="<<temp<<"\n";
               }
               if(r&1)  ans -= temp;
               else  ans += temp;
               //cout<<"ans="<<ans<<"\n";
          }while(next_permutation(b+1,b+1+n));  //下一个排列组合
          printf("%lld\n",ans);
     }
     return 0;
}