CodeForces 107 D.Crime Management(状压DP+矩阵快速幂)

52 篇文章 0 订阅
40 篇文章 1 订阅

Description

要求用 26 个大写字母构造一个长度为 n 的字符串,使得其满足m个限制条件,每个限制条件要求某个字符出现次数必须是某个数的倍数,问满足条件的字符串个数

Input

第一行两个整数 n,m 表示字符串长度和限制数,之后 m 行每行输入一字符c和一整数 num 表示字符串中 c 字符出现次数必须是num的倍数 (0n1018,0m1000)

保证 num 乘积不超过 123

Output

输出满足条件的字符串个数,结果模 12345

Sample Input

5 2
A 1
B 2

Sample Output

16

Solution

c[i] num[i] 为第 i 个限制,val[i]=j=1inum[j],假设当前字符串中第 i 个字符出现次数为cnt[i],定义当前字符串状态为 S=i=1m(cnt[c[i]]%num[i])val[i1]

考虑往当前字符串后加一个字符 x ,状态S就会变成另一个状态 T <script type="math/tex" id="MathJax-Element-42">T</script>,以此为转移矩阵,用矩阵快速幂加速转移即可

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef int M[125][125];
const int p=12345;
void Mul(M &A,M B,int n)
{
    M C;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
        {
            C[i][j]=0;
            for(int k=0;k<n;k++)C[i][j]+=A[i][k]*B[k][j]%p; 
        }
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            A[i][j]=C[i][j]%p;
}
void Pow(M &A,ll k,int n)
{
    M B;
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            B[i][j]=(i==j);
    while(k)
    {
        if(k&1)Mul(B,A,n);
        Mul(A,A,n);
        k>>=1;
    }
    for(int i=0;i<n;i++)
        for(int j=0;j<n;j++)
            A[i][j]=B[i][j];
}
#define maxn 1005
int m,c[maxn],num[maxn],val[maxn],vis[33],flag[33];
ll n;
bool check(int S)
{
    memset(flag,0,sizeof(flag));
    for(int i=1;i<=m;i++)
        if(S%val[i]/val[i-1]==0)flag[c[i]]=1;
    for(int i=0;i<26;i++)
        if(vis[i]&&!flag[i])return 0;
    return 1;
}
int main()
{
    while(~scanf("%I64d%d",&n,&m))
    {
        memset(vis,0,sizeof(vis));
        val[0]=1;
        for(int i=1;i<=m;i++)
        {
            char s[3];
            scanf("%s%d",s,&num[i]);
            c[i]=s[0]-'A';
            val[i]=val[i-1]*num[i];
            vis[c[i]]=1;
        }
        M A;
        memset(A,0,sizeof(A));
        for(int S=0;S<val[m];S++)
        {
            for(int i=0;i<26;i++)
                if(vis[i])
                {
                    int T=S;
                    for(int j=1;j<=m;j++)
                        if(c[j]==i)
                        {
                            int temp=S%val[j]/val[j-1];
                            if(temp<num[j]-1)T+=val[j-1];
                            else T-=(num[j]-1)*val[j-1];
                        }
                    A[S][T]++;
                }
        }
        Pow(A,n,val[m]);
        int ans=0;
        for(int i=0;i<val[m];i++)
            if(check(i))ans=(ans+A[0][i])%p;
        printf("%d\n",ans);
    }   
    return 0;
}
  • 0
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值