【Link】:

【Description】



输入正整数n和k(1≤n≤400,1≤k≤10),求长度为n的01串中有多少个不含长度至少 为k的回文连续子串。例如,n=k=3时只有4个串满足条件:001, 011, 100, 110。

【Solution】



状态压缩DP;

我们在做DP的时候要保证;

新添加的数,不会产生一个长度为k的回文子串;

但这还不够;

比如说

k=4

我们只保证不出现长度=4的回文子串是不是就能保证不会出现长度大于4的回文子串呢;

其实不然;

比如

10001

这是一个长度为5的回文子串;

但是它无法在去除长度为4的回文子串的情况下剔除掉;

因为

1000被认为是合法的;

下一次更新

去除首1;

只剩下000

然后再加上一个1

变成了0001

也会认为是合法的;

为什么呢?

因为我们没有剔除掉长度为3的回文子串

也即10001中间的3个0

因为长度为k的回文子串是由最左和最右两边的两个字符,加上中间一个回文子串组成的;

则只要我们能保证长度为k-2的回文子串没有出现;

就能保证新加的字符,不会组成长度为k的回文子串;

因此正确的做法应该是把长度为k和长度为k+1的回文子串都剔除掉;

这样;

长度为k+2,长度为k+3的回文子串也能剔除掉了;

从而所有长度大于等于k的回文子串也都能剔除掉了;

则设f[i][j]表示前i个数字,i-k..i这一段数字状态为j组成的符合要求的字符串的个数;

第i个数字有0和1两种情况;

分别枚举前一个状态是什么;

加上这个数字之后;

保证i-k+1..i不是回文子串且i-k..i也不是回文子串;

满足要求则这个状态就是可推的;

最后把f[n][0..2(k+1)−1]求和一下;



【NumberOf WA】



1



【Reviw】



在做DP的时候;

加上了一个字符;

要确定,它是不是一定能符合题意;

正确性要保证;

不要推到了错误的状态;



【Code】

#include <bits/stdc++.h>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define LL long long
#define rep1(i,a,b) for (int i = a;i <= b;i++)
#define rep2(i,a,b) for (int i = a;i >= b;i--)
#define mp make_pair
#define pb push_back
#define fi first
#define se second
#define ms(x,y) memset(x,y,sizeof x)
#define ri(x) scanf("%d",&x)
#define rl(x) scanf("%lld",&x)
#define rs(x) scanf("%s",x+1)
#define oi(x) printf("%d",x)
#define ol(x) printf("%lld",x)
#define oc putchar(' ')
#define all(x) x.begin(),x.end()
#define Open() freopen("F:\\rush.txt","r",stdin)
#define Close() ios::sync_with_stdio(0) typedef pair<int,int> pii;
typedef pair<LL,LL> pll; const int dx[9] = {0,1,-1,0,0,-1,-1,1,1};
const int dy[9] = {0,0,0,-1,1,-1,1,-1,1};
const double pi = acos(-1.0);
const int K = 11;
const int M = 2048;
const int N = 400;
const int MOD = 1e9+7; bool ban[K+5][M+10];
int two[K+10],dp[N+10][M+10];
bool bo[N+10][M+10]; bool ok(int k,int j){
vector <int> v,v1;
v.clear();
rep1(i,1,k){
v.pb(j%2);
j>>=1;
}
v1 = v;
reverse(all(v1));
if (v1!=v)
return false;
else
return true;
} int main(){
//Open();
//Close();
two[0] = 1;
rep1(i,1,11) two[i] = two[i-1]*2; rep1(k,1,11){
rep1(j,0,two[k]-1)
if (ok(k,j))
ban[k][j] = true;
} int T;
ri(T);
while (T--){
int n,k;
ri(n),ri(k); ms(dp,0),ms(bo,0);
dp[0][0] = 1,bo[0][0] = true; rep1(i,0,k-2)
rep1(j,0,two[k]-1)
if (bo[i][j]){
int x = (j<<1)&(two[k]-1);
//+0
bo[i+1][x] = 1;
dp[i+1][x] = (dp[i+1][x] + dp[i][j])%MOD;
//+1
x+=1;
bo[i+1][x] = 1;
dp[i+1][x] = (dp[i+1][x] + dp[i][j])%MOD;
} // i == k-1
rep1(j,0,two[k]-1)
if (bo[k-1][j]){
int x = (j<<1)&(two[k]-1);
//+0
if (!ban[k][x]){
bo[k][x] = 1;
dp[k][x] = (dp[k][x] + dp[k-1][j])%MOD;
}
//+1
x+=1;
if (!ban[k][x]){
bo[k][x] = 1;
dp[k][x] = (dp[k][x] + dp[k-1][j])%MOD;
}
} rep1(i,k,n-1)
rep1(j,0,two[k+1]-1)
if (bo[i][j]){
int x = (j<<1)&(two[k+1]-1);
int tx = x&(two[k]-1);
//+0
if (!ban[k][tx] && !ban[k+1][x]){
bo[i+1][x] = 1;
dp[i+1][x] = (dp[i+1][x] + dp[i][j])%MOD;
}
//+1
x+=1;
tx+=1;
if (!ban[k][tx] && !ban[k+1][x]){
bo[i+1][x] = 1;
dp[i+1][x] = (dp[i+1][x] + dp[i][j])%MOD;
}
} int ans = 0;
rep1(i,0,two[k+1]-1)
ans = (ans + dp[n][i])%MOD; oi(ans);puts("");
}
return 0;
}

【Uva 1633】Dyslexic Gollum的更多相关文章

  1. 【巧妙算法系列】【Uva 11464】 - Even Parity 偶数矩阵

    偶数矩阵(Even Parity, UVa 11464) 给你一个n×n的01矩阵(每个元素非0即1),你的任务是把尽量少的0变成1,使得每个元素的上.下.左.右的元素(如果存在的话)之和均为偶数.比 ...

  2. 【贪心+中位数】【UVa 11300】 分金币

    (解方程建模+中位数求最短累积位移) 分金币(Spreading the Wealth, UVa 11300) 圆桌旁坐着n个人,每人有一定数量的金币,金币总数能被n整除.每个人可以给他左右相邻的人一 ...

  3. 【UVa 10881】Piotr's Ants

    Piotr's Ants Porsition:Uva 10881 白书P9 中文改编题:[T^T][FJUT]第二届新生赛真S题地震了 "One thing is for certain: ...

  4. 【UVa 116】Unidirectional TSP

    [Link]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  5. 【UVa 1347】Tour

    [Link]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  6. 【UVA 437】The Tower of Babylon(记忆化搜索写法)

    [题目链接]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  7. 【uva 1025】A Spy in the Metro

    [题目链接]:https://uva.onlinejudge.org/index.php?option=com_onlinejudge&Itemid=8&page=show_probl ...

  8. 【Uva 11584】Partitioning by Palindromes

    [Link]:https://cn.vjudge.net/contest/170078#problem/G [Description] 给你若干个只由小写字母组成的字符串; 问你,这个字符串,最少能由 ...

  9. 【Uva 11400】Lighting System Design

    [Link]: [Description] 你要构建一个供电系统; 给你n种灯泡来构建这么一个系统; 每种灯泡有4个参数 1.灯泡的工作电压 2.灯泡的所需的电源的花费(只要买一个电源就能供这种灯泡的 ...

随机推荐

  1. webpack的理解

    webpack是一个模块打包工具,你可以使用webpack管理你的模块依赖,并编译输出模块们所需要的静态文件.它能够很好的管理.打包Web开发中所用到的HTML.Javascript.CSS以及各种静 ...

  2. ES6学习4 变量的解构赋值

    变量的解构赋值 一.数组结构赋值 1.数组结构赋值 let [a, b, c] = [1, 2, 3]; ES6 可以从数组中提取值,按照对应位置,对变量赋值. 1)  本质上,这种写法属于“模式匹配 ...

  3. 05002_Linux的安装

    1.虚拟机VMware workstation的安装 (1)下载链接:VMware workstation下载 密码:52wt: (2)双击VMware-workstation-full-10.0.2 ...

  4. ECNUOJ 2859 表达式的个数

    表达式的个数 Time Limit:5000MS Memory Limit:65536KBTotal Submit:47 Accepted:28 Description  世情薄,人情恶,雨送黄昏花易 ...

  5. Android怎样从外部跳进App

    解决问题有两个作用: 1.不用打开App直接进入某页面 2.实现App分享到外部,同一时候由外部进入App的闭环. 这个话题能够分双方面来讲.一方面是从微信进入App,还有一方面是从网页进入App. ...

  6. 【Oracle错误集锦】:ORA-00119 &amp; ORA-00132

    有时候老天就是爱和你开玩笑,昨天好不easy配置好Oracle.可以用PL/SQL正常登录使用,今天突然就不行了.而且错误十分诡异,没有提示什么错误代码.输入usernamepassword,点击登录 ...

  7. rac_grid自检提示缺少pdksh-5.2包

    原创作品,出自 "深蓝的blog" 博客,欢迎转载,转载时请务必注明下面出处,否则追究版权法律责任. 深蓝的blog:http://blog.csdn.net/huangyanlo ...

  8. Please ensure that adb is correctly located at &#39;D:\Android\android-sdk\platform-tools\adb.exe&#39; and

    1.启动任务管理器 2.找到百度安全组件杀掉进程. 3.一般都是组件给禁止了.

  9. vue9 计算属性 computed

    <!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...

  10. 智课雅思短语---一、be no exception

    智课雅思短语---一.be no exception 一.总结 一句话总结:…也不例外? …be no exception 1.经济的快速发展? the rapid development of ec ...