[COCI] Vještica (子集dp)

题意:

给定 n ( n ≤ 16 ) n(n\leq 16) n(n16)个串。你可以把串内的字母任意排列,使得将这些串插入字典树后的结点最少。
串总长不超过 1 0 6 10^{6} 106.

思路:

只有 16 16 16位,考虑状压表示集合。
D P ( ( 10111 ) 2 ) DP((10111)_{2}) DP((10111)2)表示第 1 , 3 , 4 , 5 1,3,4,5 1,3,4,5串被插入后的最小结点数。

先考虑两个串插入字典树,贡献的结点数是两个串的结点数之和减去两个串的前缀。
所以,我们考虑两个串的前缀要最长,就是每个字母的个数取 m i n min min的和。

那么考虑几个串 " a a b c " , " a a b b " . " a b c " "aabc","aabb"."abc" "aabc","aabb"."abc"。可以贪心地拿出 " a b " "ab" "ab"放前面,然后剩下 " a c " , " a b " . " c " "ac","ab"."c" "ac","ab"."c"再进行考虑。

P R E ( S ) PRE(S) PRE(S)表示集合 S S S中所有串的最长公共前缀,考虑从 ( " a c " , " a b " ) ("ac","ab") ("ac","ab") ( " c " ) ("c") ("c")合并到 ( " a c " , " a b " , " c " ) ("ac","ab","c") ("ac","ab","c")
D P ( " a a b c " , " a a b b " , " a b c " ) DP("aabc","aabb","abc") DP("aabc","aabb","abc")
= D P ( " a c " , " a b " , " c " ) + P R E ( “ a a b c ” , " a a b b " . " a b c " ) =DP("ac","ab","c")+PRE(“aabc”,"aabb"."abc") =DP("ac","ab","c")+PRE(aabc,"aabb"."abc")
= D P ( " a c " , " a b " ) + D P ( " c " ) + P R E ( " a a b c " , “ a a b b ” , " a b c " ) =DP("ac","ab")+DP("c")+PRE("aabc",“aabb”,"abc") =DP("ac","ab")+DP("c")+PRE("aabc",aabb,"abc")
= D P ( " a a b c " , " a a b b " ) − P R E ( " a b " , " a b " ) + D P ( " a b c " ) − P R E ( " a b " ) + P R E ( " a a b c " , “ a a b b ” , " a b c " ) =DP("aabc","aabb")-PRE("ab","ab")+DP("abc")-PRE("ab")+PRE("aabc",“aabb”,"abc") =DP("aabc","aabb")PRE("ab","ab")+DP("abc")PRE("ab")+PRE("aabc",aabb,"abc")
= D P ( " a a b c " , " a a b b " ) + D P ( " a b c " ) − P R E ( " a a b c " , “ a a b b ” , " a b c " ) =DP("aabc","aabb")+DP("abc")-PRE("aabc",“aabb”,"abc") =DP("aabc","aabb")+DP("abc")PRE("aabc",aabb,"abc")

那么我们考虑 S S S从它的子集 A A A转移过来,就有 D P ( S ) = D P ( A ) + D P ( S − A ) + P R E ( S ) DP(S)=DP(A)+DP(S-A)+PRE(S) DP(S)=DP(A)+DP(SA)+PRE(S)
然后就可以做子集 d p dp dp了,复杂度 o ( 3 n ) o(3^{n}) o(3n)

代码:

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;

int num1[1000005];
int getNum(int x)
{
    if(num1[x])return num1[x];
    int ans=0;
    while(x!=0)
    {
        x-=x&-x;
        ans++;
    }
    return num1[x]=ans;
}

int len[20];
int num[20][26];
char ss[1000005];
int n;
int pre(int s)
{
    int nu[26];
    for(int i=0;i<26;i++)nu[i]=1e9;
    for(int i=0;i<n;i++)if((s>>i)&1)
    {
        for(int j=0;j<26;j++)nu[j]=min(nu[j],num[i][j]);
    }
    int sum=0;
    for(int i=0;i<26;i++)sum+=nu[i];
    return sum;
}
int dp[1000005];
int main()
{
    scanf("%d",&n);
    for(int i=0;i<n;i++)
    {
        scanf("%s",ss);
        len[i]=strlen(ss);
        for(int j=0;j<len[i];j++)num[i][ss[j]-'a']++;
    }
    for(int s=1;s<1<<n;s++)
    {
        dp[s]=1e9;
        if(getNum(s)==1)
        {
            for(int i=0;i<n;i++)if((s>>i)&1)dp[s]=len[i];
            continue;
        }
        int preS=pre(s);
        for (int i=(s-1)&s;i;i=(i-1)&s)
        {
            dp[s]=min(dp[s],dp[i]+dp[s^i]-preS);
        }
    }
    printf("%d\n",dp[(1<<n)-1]+1);
    return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值