suoi62 网友跳 (暴搜+dp)

传送门

sbw太神啦orz

首先N<=20可以直接暴搜

然后玄学剪枝可以过18个点

那么N<=40的时候，就把它拆成两半分别暴搜，再用dp拼起来

对于前半段，设f[i][j]是开始高度为i，获得金币为j的方案数；对于后半段，设g[i][j]是结束高度为i，获得金币为j的方案数（离散化一下高度）

然而V<=4e7，并不能直接记

但其实每一段最多只有$2^{20}$种金币数，状压一下每一位选不选，再预处理出来这样的金币数是多少

然后统计答案，$ans=sum{g[i][j]*f[k][l]}，i>k,v[j]+v[l]>=M$

给f做一个前缀和，给每个状态按v排序，可以$O(n*2^{n/2})$统计答案

然后会MLE，注意到对于每种状态其实每个开始高度/结束高度 都只有选或者不选，做个前缀和最多也只有20

把int改成unsigned char 就行了(逃

#include<bits/stdc++.h>
#define pa pair<int,int>
#define CLR(a,x) memset(a,x,sizeof(a))
using namespace std;
typedef long long ll;
const int maxn=45,maxs=(1<<20)+2;
inline ll rd(){
    ll x=0;char c=getchar();int neg=1;
    while(c<'0'||c>'9'){if(c=='-') neg=-1;c=getchar();}
    while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar();
    return x*neg;
}
int H[maxn],N,M,v1[maxs],v2[maxs],rk1[maxs],rk2[maxs];
unsigned char f[42][maxs],g[42][maxs];
void dfs1(int x,int tar,int h,int y){
    if(x>tar) return;
    dfs1(x+1,tar,h,y);
    if(H[x]>h) f[H[x]][y|(1<<(x-1))]++,dfs1(x+1,tar,H[x],y|(1<<(x-1)));
}
void dfs2(int x,int tar,int h,int y){
    if(x<tar) return;
    dfs2(x-1,tar,h,y);
    if(H[x]<h) g[H[x]][y|(1<<(x-tar))]++,dfs2(x-1,tar,H[x],y|(1<<(x-tar)));
}

inline bool cmp1(int a,int b){return v1[a]<v1[b];}
inline bool cmp2(int a,int b){return v2[a]<v2[b];}
int main(){
    // freopen("62.in","r",stdin);
    // freopen("62.out","w",stdout);
    int i,j,k;
    N=rd(),M=rd();
    for(i=1;i<=N;i++){
        rk1[i]=H[i]=rd(),rk2[i]=rd();
    }
    sort(rk1+1,rk1+N+1);
    int t=unique(rk1+1,rk1+N+1)-rk1-1;
    for(i=1;i<=N;i++)
        H[i]=lower_bound(rk1+1,rk1+t+1,H[i])-rk1;
    int m=N>>1;
    int n1=(1<<m)-1,n2=(1<<(N-m))-1;
    dfs1(1,m,0,0);
    dfs2(N,m+1,100,0);
    
    f[0][0]=1;g[t+1][0]=1;
    for(i=1;i<=t;i++){
        for(j=0;j<=n1;j++)
            f[i][j]+=f[i-1][j];
    }
    for(i=1;i<=n1;i++){
        for(j=1;j<=m;j++)
            if((i>>(j-1))&1) v1[i]+=rk2[j];
    }for(i=1;i<=n2;i++){
        for(j=1;j<=N-m;j++)
            if((i>>(j-1))&1) v2[i]+=rk2[j+m];
    }
    for(i=1;i<=n1;i++) rk1[i]=i;
    for(i=1;i<=n2;i++) rk2[i]=i;
    sort(rk1+1,rk1+n1+1,cmp1);
    sort(rk2+1,rk2+n2+1,cmp2);
    ll ans=0;
    for(i=1;i<=t+1;i++){
        ll s=0;
        for(j=0,k=n1;j<=n2;j++){
            for(;k>=0&&v1[rk1[k]]+v2[rk2[j]]>=M;k--)
                s+=f[i-1][rk1[k]];
            ans+=s*g[i][rk2[j]];
        }
    }
    printf("%lld
",ans);
    return 0;
}