题意
不想说了…
自己看下题目把
题解
首先,考虑状压枚举列的状态。
然后进行
d
p
dp
dp,我们还有维护递增递减的状态,否则没办法转移。
我就卡死在这里了,因为这两个状态的乘积显然是会
T
T
T的,但其实维护递增递减状态的时候,我们不需要枚举每一种状态,这样就能把对于每一种列的状态中的
d
p
dp
dp转移变为
n
2
n^2
n2。
对于每一种列的状态 c o l col col, d p [ i ] [ s t a ] dp[i][sta] dp[i][sta]表示前 i i i行,以 i i i为结尾,递增递减状态为 s t a sta sta的方案数, d p [ i ] [ c m p [ i ] [ j ] & c o l ] + = d p [ j ] [ c m p [ i ] [ j ] & c o l ] dp[i][cmp[i][j]\&col]+=dp[j][cmp[i][j]\&col] dp[i][cmp[i][j]&col]+=dp[j][cmp[i][j]&col], c m p [ i ] [ j ] cmp[i][j] cmp[i][j]表示 i i i和 j j j行的比较结果,与 c o l col col作 & \& &运算,表示两行只考虑抽出来的列的比较结果。
但是这种转移不好清空和计算,因为初值是 1 1 1,我们可以令 d p dp dp表示长度 > = 2 >=2 >=2的方案数,转移的时候(简略写) d p [ i ] + = d p [ j ] + 1 dp[i]+=dp[j]+1 dp[i]+=dp[j]+1,表示第 i i i位长度 > = 2 >=2 >=2的方案数由第 j j j位长度 > = 2 >=2 >=2的方案数和第 j j j位转移过来。
清空和统计答案的时候,直接模拟一遍转移过程,每次符合转移条件的时候把答案算上,然后清 0 0 0,防止再次计算。
#include<bits/stdc++.h>
#define FOR(i,a,b) for(int i=a;i<=b;i++)
using namespace std;
typedef long long ll;
const int maxn = 6e6+1000;
const ll mod = 1e9 + 7;
int n,m;
int cmp[105][105],A[105][105],h[105];
ll dp[25][1200050];
int main(){
cin>>n>>m;
for(int i=1;i<=n;i++)for(int j=1;j<=m;j++)scanf("%d",&A[i][j]);
for(int i=1;i<=n;i++){
for(int j=1;j<=n;j++){
for(int k=1;k<=m;k++){
if(A[i][k]<A[j][k])cmp[i][j]|=(1<<(k-1));
}
}
}
ll ans=0;
for(int sta=1;sta<(1<<m);sta++){
for(int i=1;i<=n;i++)h[i]=false;
for(int i=1;i<=n;i++){
int pre,ok;
ok=true,pre=-1;
for(int j=1;j<=m;j++)if(sta&(1<<(j-1))){
if(A[i][j]>pre)pre=A[i][j];
else ok=false;
}
if(ok)h[i]=true;
ok=true,pre=100000;
for(int j=1;j<=m;j++)if(sta&(1<<(j-1))){
if(A[i][j]<pre)pre=A[i][j];
else ok=false;
}
if(ok)h[i]=true;
if(h[i])ans++;
}
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(h[i]&&h[j]){
dp[i][sta&cmp[i][j]]=(dp[i][sta&cmp[i][j]]+dp[j][sta&cmp[i][j]]+1);
}
}
}
for(int i=1;i<=n;i++){
for(int j=1;j<i;j++){
if(h[i]&&h[j])ans=(ans+dp[i][sta&cmp[i][j]]),dp[i][sta&cmp[i][j]]=0;
}
}
}
cout<<ans<<endl;
}