【codeforces 590div3】F.Yet Another Substring Reverse

题意:

现在有一个字符串 s s s 1 ≤ ∣ s ∣ ≤ 1 0 6 1 \leq |s| \leq 10^6 1s106),现在选择一个区间 [ l , r ] [l,r] [l,r],反转一次或者不反转。若某区间的字母各不相同,则该区间为完美区间。你要做的就是执行完操作后使完美区间的长度最大。(整个字符串的字母种类数不超过 20 20 20

题解:

这个问题可以转换为寻找两个不相交的完美区间使得他们的长度和最大。

首先枚举所有的完美区间,复杂度为 O ( n ⋅ s i z ) O(n \cdot siz) O(nsiz), s i z siz siz为字符种类数。

对于每个完美区间,用二进制形式来表示,若第 i i i个字母存在,则第 i i i位为 1 1 1,如 a b c f abcf abcf即可表示为:

11100100000000000000 11100100000000000000 11100100000000000000

接下来令 d p m a s k dp_{mask} dpmask代表 m a s k mask mask所有子集的长度的最大值。这个过程可以采用二进制枚举子集的方式进行 d p dp dp。复杂度 O ( 2 s i z ⋅ s i z ) O(2^{siz} \cdot {siz}) O(2sizsiz)

接下来其实就是找两个不相交的状态 m a s k 1 mask1 mask1 m a s k 2 mask2 mask2使得 d p m a s k 1 + d p m a s k 2 dp_{mask1}+dp_{mask2} dpmask1+dpmask2最大化,枚举 m a s k 1 mask1 mask1, m a s k 2 mask2 mask2其实就是 m a s k 1 mask1 mask1按位取反。枚举更新答案即可。

代码:

#include

using namespace std;
int dp[1<<21];//dp[mask]表示mask所有子状态的最大值
int cal(int x){
    int ans=0;
    while(x){
        ans+=x%2;
        x/=2;
    }
    return ans;
}
int main(){
    string s;
    cin>>s;
    int l=(int)s.length();
    memset(dp,0,sizeof(dp));
    for(int i=0;i<l;i++){//枚举s[i:j]满足s[i:j]各个字母不相同
        bool vis[20];
        memset(vis,0,sizeof(vis));
        int mask=0;
        for(int j=0;i+j<l;j++){
            if(vis[s[i+j]-'a']) break;
            vis[s[i+j]-'a']=1;
            mask|=(1<<(s[i+j]-'a'));
            dp[mask]=cal(mask);
            //cout<
        }
    }
    for(int mask=0;mask<(1<<20);mask++){
        for(int j=0;j<20;j++){
            if((mask>>j)&1){
                dp[mask]=max(dp[mask],dp[mask^(1<<j)]);//枚举子集更新最大值
                //cout<
            }
        }
    }
    int ans=0;
    for(int mask=0;mask<(1<<20);mask++){
        if(dp[mask]==cal(mask)){
            int now=~mask&((1<<20)-1);
            //cout<
            ans=max(ans,dp[mask]+dp[now]);
        }
    }
    cout<<ans<<endl;
    return 0;
}

你可能感兴趣的:(OJ--codeforces,动态规划--状压DP)