【BZOJ5336】[TJOI2018]party(动态规划)

题面

BZOJ

洛谷

题解

这题好神仙啊。。。

考虑普通的\(LCS\)的\(dp\),\(f[i][j]=\max\{f[i-1][j],f[i][j-1],f[i-1][j-1]+(A_i==B_j)\}\)

发现对于固定的\(i\)而言,随着\(j\)的增长,相邻的两个数之间的差不超过\(1\),因此直接考虑一个\(2^k\)的状态记录差分的结果。

那么设\(f[i][S]\)表示当前考虑到了第\(i\)位,\(LCS\)的\(dp\)的差分数组为\(S\)的方案数。转移的时候枚举放那个字母就好了。因为还有连续字符的限制,所以再加一维就没有问题了。

这里转移因为会涉及差分数组的修改,所以可以提前把差分数组的转移给预处理出来,

然后大力\(dp\)就好了。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
#define MOD 1000000007
#define MAX 1010
char ch[20],QwQ[5]="NOI";
int f[2][1<<15][3],tr[1<<15][3],cnt[1<<15];
int n,k,ans[MAX];
void add(int &x,int y){x+=y;if(x>=MOD)x-=MOD;}
int dp[16],tmp[16];
int Get(int S,int c)
{
for(int i=1;i<=k;++i)dp[i]=dp[i-1]+((S>>(i-1))&1);
for(int i=1;i<=k;++i)
if(QwQ[c]==ch[i])tmp[i]=dp[i-1]+1;
else tmp[i]=max(dp[i],tmp[i-1]);
int ret=0;
for(int i=1;i<=k;++i)ret|=(tmp[i]-tmp[i-1])<<(i-1);
return ret;
}
int main()
{
scanf("%d%d",&n,&k);scanf("%s",ch+1);
for(int i=0;i<1<<k;++i)
for(int c=0;c<3;++c)
tr[i][c]=Get(i,c);
for(int i=0;i<(1<<k);++i)cnt[i]=cnt[i>>1]+(i&1);
f[0][0][0]=1;
for(int i=1,nw=1,pw=0;i<=n;++i,nw^=1,pw^=1)
{
memset(f[nw],0,sizeof(f[nw]));
for(int j=0;j<1<<k;++j)
for(int p=0;p<3;++p)
for(int c=0;c<3;++c)
{
int np;
if(c==0)np=1;else if(c==1)np=p==1?2:0;else if(c==2)np=p==2?3:0;
if(np==3)continue;
add(f[nw][tr[j][c]][np],f[pw][j][p]);
}
}
for(int i=0;i<1<<k;++i)
for(int c=0;c<3;++c)
add(ans[cnt[i]],f[n&1][i][c]);
for(int i=0;i<=k;++i)printf("%d\n",ans[i]);
return 0;
}

【BZOJ5336】[TJOI2018]party(动态规划)的更多相关文章

  1. BZOJ5336 TJOI2018 party 【状压DP】*

    BZOJ5336 TJOI2018 party Description 小豆参加了NOI的游园会,会场上每完成一个项目就会获得一个奖章,奖章 只会是N, O, I的字样.在会场上他收集到了K个奖章组成 ...

  2. BZOJ5336: [TJOI2018]party

    BZOJ5336: [TJOI2018]party https://lydsy.com/JudgeOnline/problem.php?id=5336 分析: 好题. 正常的思路是设\(f[i][j] ...

  3. [模板] dp套dp && bzoj5336: [TJOI2018]party

    Description Problem 5336. -- [TJOI2018]party Solution 神奇的dp套dp... 考虑lcs的转移方程: \[ lcs[i][j]=\begin{ca ...

  4. DP套DP

    DP套DP,就是将内层DP的结果作为外层DP的状态进行DP的方法. [BZOJ3864]Hero meet devil 对做LCS的DP数组差分后状压,预处理出转移数组,然后直接转移即可. tr[S] ...

  5. 【BZOJ5337】[TJOI2018]str(动态规划,哈希)

    [BZOJ5337][TJOI2018]str(动态规划,哈希) 题面 BZOJ 洛谷 题解 就很呆... 显然按层\(dp\),如果能够匹配上就进行转移,直接哈希判断是否能够匹配就好了... #in ...

  6. BZOJ5337 [TJOI2018] 碱基序列 【哈希】【动态规划】

    题目分析: 这道题的难点在于要取模,而题面没有写. 容易想到一个O(1E7)的dp.KMP或者哈希得到相关位置然后对于相关位置判断上一个位置有多少种情况. 代码: #include<bits/s ...

  7. BZOJ5336:[TJOI2018]游园会——题解

    https://www.lydsy.com/JudgeOnline/problem.php?id=5336 https://www.luogu.org/problemnew/show/P4590 小豆 ...

  8. 增强学习(三)----- MDP的动态规划解法

    上一篇我们已经说到了,增强学习的目的就是求解马尔可夫决策过程(MDP)的最优策略,使其在任意初始状态下,都能获得最大的Vπ值.(本文不考虑非马尔可夫环境和不完全可观测马尔可夫决策过程(POMDP)中的 ...

  9. 简单动态规划-LeetCode198

    题目:House Robber You are a professional robber planning to rob houses along a street. Each house has ...

随机推荐

  1. MySQL 优化集锦

    case 1: 如果筛选or条件有多个的时候,应该将最好判断的放在最前面,将最不好判断的放在最后面 比如,有一个学生表,想要找出其中年龄是20岁,住址中包含666这个数字的记录.可以下面这两个方案: ...

  2. 使用ajax请求后端程序时,关于目标程序路径问题

    这里涉及到和PHP中类似的问题,有待更新!!!

  3. #Leetcode# 985. Sum of Even Numbers After Queries

    https://leetcode.com/problems/sum-of-even-numbers-after-queries/ We have an array A of integers, and ...

  4. MySql concat与字符转义

    mysql函数之四:concat() mysql 多个字段拼接 - duanxz - 博客园https://www.cnblogs.com/duanxz/p/5098875.html mysql 多个 ...

  5. 关于PHP的那些坑

    因为PHP是弱类型语言,常常会发生许多意想不到的问题,所以,我们再次一一细数这些我们踏过的坑!!! 1) foreach中自动回将key为数值的转化成整型,造成无法匹配 function transl ...

  6. 解决多人开发时使用window.onload的覆盖问题

    通用型小函数:解决多人开发时,同时使用window.onload事件所出现的后面的window.onload函数覆盖前面一个window.onload函数的问题. function addLoadEv ...

  7. Hbase 架构体系

    有2个节点进程,一个是master,另一是regionserver.

  8. python之路--操作系统介绍,进程的创建

    一 .  操作系统的作用: 1:隐藏丑陋复杂的硬件接口,提供良好的抽象接口 2:管理.调度进程,并且将多个进程对硬件的竞争变得有序 二 多道技术: 所谓多道程序设计技术,就是指允许多个程序同时进入内存 ...

  9. app自动化测试Appium+python

    一.node.js安装 https://nodejs.org/en/download/ ##一直下一步 ###cmd查看 二.  .NET Framework安装 https://www.micros ...

  10. 莫烦scikit-learn学习自修第一天【scikit-learn安装】

    1. 机器学习的分类 (1)有监督学习(包括分类和回归) (2)无监督学习(包括聚类) (3)强化学习 2. 安装 (1)安装python (2)安装numpy >=1.6.1 (3)安装sci ...