Kattis - monotony状压dp

题意

不想说了…
自己看下题目把

题解

首先，考虑状压枚举列的状态。
然后进行$dp$，我们还有维护递增递减的状态，否则没办法转移。
我就卡死在这里了，因为这两个状态的乘积显然是会$T$的，但其实维护递增递减状态的时候，我们不需要枚举每一种状态，这样就能把对于每一种列的状态中的$dp$转移变为$n^2$。

对于每一种列的状态$col$，$dp[i][sta]$表示前$i$行，以$i$为结尾，递增递减状态为$sta$的方案数，$dp[i][cmp[i][j]\&col]+=dp[j][cmp[i][j]\&col]$，$cmp[i][j]$表示$i$和$j$行的比较结果，与$col$作$\&$运算，表示两行只考虑抽出来的列的比较结果。

但是这种转移不好清空和计算，因为初值是$1$，我们可以令$dp$表示长度$>=2$的方案数，转移的时候(简略写)$dp[i]+=dp[j]+1$，表示第$i$位长度$>=2$的方案数由第$j$位长度$>=2$的方案数和第$j$位转移过来。

清空和统计答案的时候，直接模拟一遍转移过程，每次符合转移条件的时候把答案算上，然后清$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;
}