[COCI] Vještica （子集dp）

最新推荐文章于 2024-07-24 21:16:11 发布

Zimba_

最新推荐文章于 2024-07-24 21:16:11 发布

阅读量188

点赞数

CC 4.0 BY-SA版权

分类专栏： dp

本文链接：https://blog.csdn.net/weixin_43785386/article/details/108967450

dp 专栏收录该内容

4 篇文章

订阅专栏

本文讲解了一种高效的算法，通过状压DP和贪心策略，针对长度不超过10^6的n个字符串，最小化插入字典树后的结点数量。核心思路是利用集合状态和最长公共前缀计算。最终实现了一个O(3^n)复杂度的子集DP解决方案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意：

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

思路：

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

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

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

用 $P R E (S)$ 表示集合 $S$ 中所有串的最长公共前缀，考虑从 $(" a c ", " a b ")$ 和 $(" c ")$ 合并到 $(" a c ", " a b ", " c ")$ 。
$D P (" a a b c ", " a a b b ", " a b c ")$
$= D P (" a c ", " a b ", " c ") + P R E (“ a a b c ”, " a a b b " . " a b c ")$
$= D P (" a c ", " a b ") + D P (" c ") + P R E (" a a b c ", “ a a b b ”, " a b c ")$
$= 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 ")$
$= 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 ")$

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

代码：

#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;
}