@51nod - 1196/1197/1198@ 字符串的数量
@description@
用 N 个不同的字符(编号1 - N),组成一个字符串,有如下要求:
(1) 对于编号为 i 的字符,如果2 * i > n,则该字符可以作为结尾字符。如果不作为结尾字符而是中间的字符,则该字符后面可以接任意字符。
(2) 对于编号为 i 的字符,如果2 * i <= n,则该字符不可以作为结尾字符。作为中间字符,那么后面接的字符编号一定要 >= 2 * i。
问有多少长度为M且符合条件的字符串,由于数据很大,只需要输出该数Mod 10^9 + 7的结果。
例如:N = 2,M = 3。则abb, bab, bbb是符合条件的字符串,剩下的均为不符合条件的字符串。
input
第1行:一个数T,表示后面用作输入测试的数的数量。(1 <= T <= 20)
第2 - T + 1行:每行2个数,N, M中间用空格分割,N为不同字符的数量,M为字符串的长度。(2 <= N, M <= 10^18)
output
共T行,对应符合条件的字符串的数量。由于数据很大,只需要输出该数Mod 10^9 + 7的结果。
sample input
1
6 3
sample output
73
@solution@
什么鬼题目。。。
@part - 1@
不难写出一个很 naive 的 dp:定义 dp[i][j] 长度为 i,末尾字符为 j 的方案数。
对 2*j > n 还是 <= n 进行分类讨论就可以得到转移式。
不难发现对于 2*j > n 的 dp 值是相同的,我们令 f[i] 表示这个值。
于是我们可以发现 dp[i][j] 总是可以表示为 f[i-1], f[i-2], ... 的线性组合。
考虑这个性质的实际意义:我们可以在一个字符后面添加若干的字符最终变为以 2*x > n 这一类字符结尾。
我们以 2*x > n 这一类字符作为隔断,将最终的字符串分成若干以这一类字符结尾的段。可以发现每一段的长度是 O(log) 的,因为后一个字符至少是前一个字符的两倍。
定义 g[i][x] 表示以字符 x 开头,产生长度 i 的段的方案数。于是根据定义就有 g[1][x] = 1,对于 2*x > n。
同时也有 f[i] = ∑g[j][x]*f[i-j]。因为段的长度不会超过 O(log),所以可以预处理 g 然后矩乘。
@part - 2@
先来考虑矩乘这一部分,顺便把我们上述的内容整理一下。
定义 f[i] 表示问题所想我们求解的答案,再定义 g[i][x] 表示以字符 x 开头,最后一位且只有最后一位是 2*x > n 这类字符,长度为 i 的字符串方案数。
于是有转移 \(f[i] =\sum_{j, k}g[j][k]*f[i-j]\)。
注意到 g[i][x] 这种字符串,每多一位,末尾字符就会变为前一个的字符两倍。故它的长度是 O(log) 的。
所以我们在预处理好 g 之后,可以写一个 O(log^3 n*log m) 的矩乘求出 f[n]。
但是这样有些卡,我们可能需要进一步的优化。
使用一种名为“常系数齐次线性递推”的黑科技即可。(可以点击我的这篇博客查看详情)
不过这道题可以直接暴力多项式乘法与多项式取模,因为 log 也不大,说不定 fft 的常数反而会跑得慢些。
最后时间可以优化到 O(log^2 n*log m)
@part - 3@
考虑怎么预处理好 g。不难发现 g 有如下的转移式:
\]
然后开始玄幻起来了。我们使用万能的归纳法,尝试证明 g 始终满足以下形式:
P1(x) & 1 \le x \le \lfloor{\frac{a}{2}}\rfloor \\
P2(x) & \lfloor{\frac{a}{2}}\rfloor + 1 \le x \le a
\end{cases}\]
即一个分段的多项式函数。
假设对于所有 <= i 的都成立,尝试证明对 i+1 成立。
如果 \(\lfloor{\frac{a}{2}}\rfloor + 1 \le 2*x \le a\),则:
\]
前一项是一个常数,而后一项呢?我们不妨令 \(P2(x) = a_0+a_1x+...\)。
于是:
\]
然后你发现那是个自然数幂和。然后根据你所了解到的知识,自然数幂和可以表示为多项式形式。然后得证。
如果 \(1 \le 2*x \le \lfloor{\frac{a}{2}}\rfloor\),证明是类似的,不再赘述。
关于自然数幂和,这里有两篇(我认为可供参考的)博客:blog1, blog2。
找时间可以专门研究这玩意儿(咕咕咕)。
然后就没有然后了。预处理一下自然数幂和的系数之类的即可。
具体可以看代码。
@accepted code@
#include<cstdio>
typedef long long ll;
const int MOD = int(1E9) + 7;
const int MSIZE = 64;
int pow_mod(int b, int p) {
int ret = 1;
while( p ) {
if( p & 1 ) ret = 1LL*ret*b%MOD;
b = 1LL*b*b%MOD;
p >>= 1;
}
return ret;
}
int len;
struct poly{
int k[2*MSIZE + 5];
poly() {for(int i=0;i<=2*MSIZE;i++) k[i] = 0;}
int get_val(int x) {
int ret = 0;
for(int i=len;i>=0;i--)
ret = (1LL*ret*x%MOD + k[i])%MOD;
return ret;
}
friend poly operator +(poly A, poly B) {
poly C;
for(int i=0;i<MSIZE;i++)
C.k[i] = (A.k[i] + B.k[i])%MOD;
return C;
}
friend poly operator -(poly A, poly B) {
poly C;
for(int i=0;i<MSIZE;i++)
C.k[i] = (A.k[i] + MOD - B.k[i])%MOD;
return C;
}
friend poly operator *(int k, poly B) {
poly C;
for(int i=0;i<MSIZE;i++)
C.k[i] = 1LL*k*B.k[i]%MOD;
return C;
}
friend poly mul_mod(poly A, poly B, poly C) {
poly D;
for(int i=0;i<=2*len;i++)
for(int j=0;j<=i;j++)
D.k[i] = (D.k[i] + 1LL*A.k[j]*B.k[i-j]%MOD)%MOD;
for(int i=2*len;i>=len;i--) {
for(int j=0;j<=len;j++)
D.k[i-len+j] = (D.k[i-len+j] + MOD - 1LL*D.k[i]*C.k[j]%MOD)%MOD;
}
return D;
}
void debug() {
for(int i=0;i<2*MSIZE;i++)
printf("%d ", k[i]);
puts("");
}
}pw[MSIZE + 5], a[MSIZE + 5];
int comb[MSIZE + 5][MSIZE + 5];
void init() {
for(int i=0;i<MSIZE;i++) {
comb[i][0] = 1;
for(int j=1;j<=i;j++)
comb[i][j] = (comb[i-1][j] + comb[i-1][j-1])%MOD;
}
for(int i=0;i<MSIZE;i++) {
for(int j=0;j<=i+1;j++)
pw[i].k[j] = comb[i+1][j];
for(int j=0;j<i;j++)
pw[i] = (pw[i] - comb[i+1][j]*pw[j]);
pw[i] = 1LL*pow_mod(i+1, MOD-2)*pw[i];
}
for(int i=0;i<MSIZE;i++)
for(int j=0;j<=i+1;j++)
for(int k=0;k<=j;k++)
a[i].k[k] = (a[i].k[k] + 1LL*pw[i].k[j]*(1LL*comb[j][k]*(1LL*pow_mod(2, k)*pow_mod(MOD-1, j-k)%MOD)%MOD)%MOD)%MOD;
}
int get_pw(ll x, int k) {return pw[k].get_val(x%MOD);}
int get_pw(ll l, ll r, int k) {return (get_pw(r, k) + MOD - get_pw(l - 1, k))%MOD;}
int get_val(poly A, ll l, ll r) {
int ret = 0;
for(int i=0;i<len;i++)
ret = (ret + 1LL*get_pw(l, r, i)*A.k[i]%MOD)%MOD;
return ret;
}
poly trans(poly A, ll r) {
poly B; B.k[0] = get_val(A, 0, r);
for(int i=0;i<len;i++)
B = B - A.k[i]*a[i];
return B;
}
poly get_M(ll n) {
poly A, B, M; B.k[0] = M.k[len] = 1;
for(int i=0;i<len;i++,n/=2){
int sA = get_val(A, 1, n/2), sB = get_val(B, n/2 + 1, n);
M.k[len-i-1] = (MOD - (sA + sB)%MOD)%MOD;
A = trans(A, n/2), A.k[0] = (A.k[0] + sB)%MOD;
B = trans(B, n);
}
return M;
}
poly p_pow(poly M, ll p) {
poly ret, b; b.k[1] = 1, ret.k[0] = 1;
while( p ) {
if( p & 1 ) ret = mul_mod(ret, b, M);
b = mul_mod(b, b, M);
p >>= 1;
}
return ret;
}
int solve(ll n, ll m) {
ll tmp = n; for(len = 0; tmp; len++, tmp /= 2);
return p_pow(get_M(n), m+len-1).k[len-1];
}
int main() {
init(); int T; scanf("%d", &T);
for(int i=1;i<=T;i++) {
ll n, m; scanf("%lld%lld", &n, &m);
printf("%d\n", solve(n, m));
}
}
@details@
我可能永远也想不到把 g 表示成分段多项式函数的形式这一步。。。
写起来还是比较愉快的吧。。。毕竟没有毒瘤的多项式操作。。。所有多项式操作都可以直接暴力搞。。。
矩乘过不了是真的令人生气。。。
@51nod - 1196/1197/1198@ 字符串的数量的更多相关文章
- 51nod1196 字符串的数量
用N个不同的字符(编号1 - N),组成一个字符串,有如下要求:(1) 对于编号为i的字符,如果2 * i > n,则该字符可以作为结尾字符.如果不作为结尾字符而是中间的字符,则该字符后面可以接 ...
- [51nod1197]字符串的数量 V2
用N个不同的字符(编号1 - N),组成一个字符串,有如下要求: (1) 对于编号为i的字符,如果2 * i > n,则该字符可以作为结尾字符.如果不作为结尾字符而是中间的字符,则该字符后面可以 ...
- 用map来统计数组中各个字符串的数量
1.背景 想要统计这一个字符串数组中每一个非重复字符串的数量,使用map来保存其key和value.这个需求在实际开发中经常使用到,我以前总是新建一个空数组来记录不重复字符串,并使用计数器计数,效率低 ...
- 51nod 1196 字符串的数量(DP+数论?)
这题好像是神题...V1 V2 V3分别涵盖了51nod 5级算法题 6级算法题 难题 讨论区的曹鹏神牛好强啊...一种做法切了V1 V2 V3,而且做法是一步一步优化的 还没去看优化的部分,未优化已 ...
- 51Nod 1196 字符串的数量
用N个不同的字符(编号1 - N),组成一个字符串,有如下要求: (1) 对于编号为i的字符,如果2 * i > n,则该字符可以作为结尾字符.如果不作为结尾字符而是中间的字符,则该字符后面可以 ...
- 51nod【1196】字符串的数量
超级神题! 有n种字符,若此种字符的编号( \(1\) ~ \(n\)),\(i*2>n\),则他后面可接任意字符.若不是,则他后面接的字符编号至少要是他的两倍. 问长度为m的字符串的个数. 这 ...
- 51nod 1197 字符串的数量 V2(矩阵快速幂+数论?)
接上一篇,那个递推式显然可以用矩阵快速幂优化...自己随便YY了下就出来了,学了一下怎么用LaTeX画公式,LaTeX真是个好东西!嘿嘿嘿 如上图.(刚画错了一发...已更新 然后就可以过V2了 or ...
- 51nod 1092 回文字符串 (dp)
http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1092 这个题是poj-3280的简化版,这里只可以增加字符,设 dp[i ...
- 51Nod 1003 阶乘后面0的数量(数学,思维题)
1003 阶乘后面0的数量 基准时间限制:1 秒 空间限制:131072 KB 分值: 5 难度:1级算法题 n的阶乘后面有多少个0? 6的阶乘 = 1*2*3*4*5*6 = 720 ...
随机推荐
- 在Spring应用中创建全局获取ApplicationContext对象
在Spring应用中创建全局获取ApplicationContext对象 1.需要创建一个类,实现接口ApplicationContextAware的setApplicationContext方法. ...
- python基础--基础数据类型
1.输入 python2中的输入: 关键字:input() --> 需要人为的告诉input你输入的是哪种类型的数据(声明数据类型) >>> name = input('pl ...
- Apache HttpComponents 工具类 [ HttpUtil ]
pom.xml <dependency> <groupId>org.apache.httpcomponents</groupId> <artifactId&g ...
- NoSQL最新现状和趋势:云NoSQL数据库将成重要增长引擎
NoSQL最早起源于1998年,但从2009年开始,NoSQL真正开始逐渐兴起和发展.回望历史应该说NoSQL数据库的兴起,完全是十年来伴随互联网技术,大数据数据的兴起和发展,NoSQL在面临大数据场 ...
- leetcode 198-234 easy
198. House Robber 相邻不能打劫,取利益最大化. 思想:当前值和前一个和的总数 与 前一个和 做大小比较,取最大值,重复该步骤. class Solution { publ ...
- HDU 1690 最短路
#include<stdio.h> #include<string.h> #include<queue> #include<algorithm> usi ...
- Sql基本知识回顾
一. SQL 基本语句 SQL 分类: DDL —数据定义语言 (Create , Alter , Drop , DECLARE) DML —数据操纵语言 (Select , Delete , Upd ...
- jQuery ajax请求struts action实现异步刷新
第一步:导入相关jar包,本样例需导入struts相关jar包,json-lib.jar,gson-2.1.jar可以任意选择,但是这里需要都导入,因为为了做测试,两种jar包的转换方式都用到了. 第 ...
- pycharm 2017 序列号失效问题解决(2016-2017版本都有效)
pycharm 序列号失效问题解决 this license BIG3CLIK6F has been cancelled 具体如下: 对,没错,这个激活码本来可以使用到2018年的,但是,忽然间 ...
- ubuntukylin-16.04安装
目录 安装 虚拟机安装 系统安装 vmware tools 安装 优化 1.切换软件源 2.软件源和软件更新 3.vim插件 4.安装Sublime Text 3 5.安装福昕阅读器 6.安装Anac ...