hdu5745 La Vie en rose dp+bitset优化

多校联合第12题

贴一个官方题解 http://bestcoder.hdu.edu.cn/blog/2016-multi-university-training-contest-2-solutions-by-zimpha/

题解的DP很容易理解，然后用bitset把i优化掉，也就是现在dp[j][k]存的是一个100000位的二进制值，第i位为1就代表官方题解中的dp[i][j][k]=true，0则为false。同理s[i]=p[i]这部分也可以直接用bitset表示。

以后遇到dp值为bool，而且状态转移只需要进行与或非运算的时候，估摸着都可以用bitset优化的感觉。。。

然而此题正解用时5s多，纯暴力跑只需2.5s。。。。（瞎）

#include <bits/stdc++.h>
#define fi first
#define se second
using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
const int inf = 0x3f3f3f3f;


const int maxn=100000+10;
int n,m;
char s[maxn], p[maxn];
bitset<maxn> letter[26],l[26],r[26];
bitset<maxn> dp[2][3];
bitset<maxn> ans,aux;


void get_letter(){
    memset(letter,0,sizeof letter);
    for (int i = 0; i < n; ++i)
    {
        letter[s[i]-'a'][i]=1;
    }
    for (int i = 0; i < 26; ++i)
    {
        r[i]=letter[i]<<1;
        l[i]=letter[i]>>1;
    }
}


int main() {
#ifndef ONLINE_JUDGE
    freopen("test.in", "r", stdin);
#endif


    int T;
    scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        scanf("%s%s",s,p);


        get_letter();


        int cur=0, pre=1;
        memset(dp[pre],0,sizeof dp[pre]);


        dp[pre][1]=letter[p[0]-'a'];
        dp[pre][2]=letter[p[0]-'a']>>1;


        for (int j = 1; j < m; ++j)
        {
            memset(dp[cur],0,sizeof dp[cur]);
            int ch=p[j]-'a';


            aux=(dp[pre][0]|dp[pre][1])<<1;
            dp[cur][0]=(dp[pre][2]<<1) & r[ch];
            dp[cur][1]=aux & letter[ch];
            dp[cur][2]=aux & l[ch];


            swap(cur,pre);
        }


        ans=(dp[pre][0]|dp[pre][1])>>(m-1);


        for(int i=0; i<n; i++){
            printf("%d",ans.test(i));
        }
        puts("");
    }
}

好吧数据改了，现在暴力已废，而且代码需要少些位运算才能过