题目描述

小豆参加了\(NOI\)的游园会,会场上每完成一个项目就会获得一个奖章,奖章 只会是\(N\), \(O\), \(I\)的字样。在会场上他收集到了\(K\)个奖章组成的串。

兑奖规则是奖章串和兑奖串的最长公共子序列长度为小豆最后奖励的等级。

现在已知兑奖串长度为\(N\),并且在兑奖串上不会出现连续三个奖章为\(NOI\),即奖章中不会出现子串\(NOI\)。

现在小豆想知道各个奖励等级会对应多少个不同的合法兑奖串。

题解

我们可以先考忽略那些奇奇怪怪的限制条件,直接考虑如何统计序列数。

我们先考虑\(LCS\)的dp方法。

\[dp[i][j]=max(dp[i-1][j],dp[i][j-1],dp[i-1][j-1]+(s[i]==t[j]))
\]

我们如果把第二维单独拿出来考虑,我们可以发现第二维是一个单调不降的数列而且前后两项的差是小于等于1的。

因为第二维的长度非常小,所以我们可以把第二维差分之后状压起来作为我们的状态。

而且有一个好处,我们知道了这个状态,就能发当前\(dp\)值还原出来,假设我们知道了下一个字符时什么,我们就可以知道转移后的状态是什么了。

所以我们预处理\(g[i][j]\)表示当前状态为\(i\),新加入字符为\(j\),能够转移的状态。

然后我们再设\(dp[i][j][k]\)表示做到第\(i\)为,当前状态为\(j\),匹配到k个字符(这个是判断那个特殊限制用的),随便转移一下就好了。

代码

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
const int mod=1e9+7;
int tran[200],g[3][1<<16],dp[2][1<<16][3],n,k,nw[20],a[20],ma,ans[20];
char s[20];
inline int rd(){
int x=0;char c=getchar();bool f=0;
while(!isdigit(c)){if(c=='-')f=1;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
return f?-x:x;
}
inline void MOD(int &x){x=x>=mod?x-mod:x;}
int main(){
n=rd();k=rd();
tran['N']=0;tran['O']=1;tran['I']=2;
scanf("%s",s+1);
for(int i=1;i<=k;++i)a[i]=tran[(int)s[i]];
ma=(1<<k+1);
for(int o=0;o<3;++o){
for(int i=0;i<ma;++i){
for(int j=0;j<=k;++j){
nw[j]=(i&(1<<j))!=0;
if(j)nw[j]+=nw[j-1];
}
for(int j=k;j>=1;--j){
nw[j]=max(nw[j],nw[j-1]);
if(a[j]==o)nw[j]=max(nw[j],nw[j-1]+1);
}
for(int j=1;j<=k;++j)nw[j]=max(nw[j],nw[j-1]);
int s=0;
for(int j=k;j>=1;--j)nw[j]=nw[j]-nw[j-1];
for(int j=0;j<=k;++j)s|=(1<<j)*nw[j];
g[o][i]=s;
}
}
int now=1,pre=0;
dp[now][0][0]=1;
for(int i=0;i<n;++i){
swap(now,pre);
memset(dp[now],0,sizeof(dp[now]));
for(int j=0;j<ma;++j)
for(int l=0;l<3;++l)if(dp[pre][j][l])
for(int p=0;p<3;++p){
int no=(p==l)?l+1:(p==0);
if(no==3)continue;
MOD(dp[now][g[p][j]][no]+=dp[pre][j][l]);
}
}
for(int i=0;i<ma;++i){
for(int j=0;j<=k;++j){
nw[j]=(i&(1<<j))!=0;
if(j)nw[j]+=nw[j-1];
}
MOD(ans[nw[k]]+=(1ll*dp[now][i][0]+dp[now][i][1]+dp[now][i][2])%mod);
}
for(int i=0;i<=k;++i)printf("%d\n",ans[i]);
return 0;
}

TJOI2018Party的更多相关文章

随机推荐

  1. Mysql使用event,类似oracle job

    MySQL从5.1开始支持event功能,类似oracle的job功能.有了这个功能之后我们就可以让MySQL自动的执行数据汇总等功能,不用像以前需要操作的支持了.如linux crontab功能. ...

  2. 如何将数据库中的值经过servlet传入到jsp页面,并且用EL表达式显示出值

    方法一:通过id查询某一数据库表中具体的行,将值封装在相应的对象中,如下面的对象Notice servlet中 String noticeId=request.getParameter("n ...

  3. jsp内置对象-config对象

    1.概念:config对象中存储了一些Servlet初始化的数据结构,当Servlet初始化时,JSP容器通过config对象将这些信息传递给这个Servlet.一般在web.xml文件中配置Serv ...

  4. Spring boot入门(三):SpringBoot集成结合AdminLTE(Freemarker),利用generate自动生成代码,利用DataTable和PageHelper进行分页显示

    关于SpringBoot和PageHelper,前篇博客已经介绍过Spring boot入门(二):Spring boot集成MySql,Mybatis和PageHelper插件,前篇博客大致讲述了S ...

  5. Netty 核心容器之ByteBuf 结构详解

    原文链接 Netty 核心容器之ByteBuf 结构详解 代码仓库地址 Java的NIO模块提供了ByteBuffer作为其字节存储容器,但是这个类的使用过于复杂,因此Netty实现了ByteBuf来 ...

  6. 深入理解group by 语句的执行顺序 from→where→group by→select(含聚合函数)

    由于之前没有对group by 语句的执行顺序(执行原理)做深入的了解,所以导致在实际应用过程中出现了一些问题.举个简单的粟子,比如一个表testA中的所有数据如下图: 我现在想从testA中查询us ...

  7. <%@ Register TagPrefix="uc1" TagName="user" Src="../Control/user.ascx" %>什么意思?

    TagPrefix定义控件位置的命名空间.有了命名空间制约,就可以在同一个网页里使用不同功能的同名控件. TagName指向所使用控件的名字.在同一个命名空间里的控件名是唯一的.控件名一般都表明控件的 ...

  8. ext前后台数据传输的标准化

    一.标准化的数据传输是什么 这里所说的标准化主要是指,使用代理提交数据时,格式必须统一化.标准化,而服务器返回的数据格式也必须是标准化的数据. 简言之,使用代理提交数据时,前台--->后台,后台 ...

  9. anaconda 环境新建/删除/拷贝 jupyter notebook上使用python虚拟环境 TensorFlow

    naconda修改国内镜像源 国外网络有时太慢,可以通过配置把下载源改为国内的通过 conda config 命令生成配置文件,这里使用清华的镜像: https://mirrors.tuna.tsin ...

  10. 如何快速搭建一个基于ServiceStack框架的web服务

    ServiceStack是一个高性能的.NET Web Service 平台,能够简化开发高性能的REST (支持JSON,XML,JSV,HTML,MsgPack,ProtoBuf,CSV等消息格式 ...