hdu 3658 How many words
来源:互联网 发布:手机跳墙软件 编辑:程序博客网 时间:2024/06/05 06:00
点击打开hdu 3658
思路: 递推+矩阵快速幂
分析:
1 题目的意思是在52个英文字母里面选择m个字母组成一个字符串,满足以下两个条件。第一是相邻的两个字符的ASCLL码的绝对值小于等于32,第二至少要有一对的字符的绝对值为32
2 那么不考虑第二个条件的时候,我们可以求出所有的符合的个数。假设f(n)(j)表示的是前n个字符最后一个字符为j,那么我们可以求出所有满足第一个条件的所有个数。因为至少需要有一对相邻的字符的绝对值为32,那么我们只要把第一次求出的所有的个数减去“相邻的两个字符的ASCLL码的绝对值小于等于31”的即可
3 那么我们考虑“相邻的两个字符的ASCLL码的绝对值小于等于32”这种情况,f(n)(j) = Σ(f(n-1)(k)) , abs(j-k) <= 32
那么我们可以构造出如下的矩阵
4 那么相邻的两个字符的ASCLL码的绝对值小于等于31就和上面的类似
代码:
/************************************************ * By: chenguolin * * Date: 2013-08-31 * * Address: http://blog.csdn.net/chenguolinblog * ************************************************/#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;typedef long long int64;const int MOD = 1e9+7;const int N = 52; struct Matrix{ int64 mat[N][N]; Matrix operator*(const Matrix &m)const{ Matrix tmp; for(int i = 0 ; i < N ; i++){ for(int j = 0 ; j < N ; j++){ tmp.mat[i][j] = 0; for(int k = 0 ; k < N ; k++) tmp.mat[i][j] += mat[i][k]*m.mat[k][j]%MOD; tmp.mat[i][j] %= MOD; } } return tmp; }};void init(Matrix &m1 , Matrix &m2){ // m1 memset(m1.mat , 0 , sizeof(m1.mat)); int x = 25; for(int i = 0 ; i < 26 ; i++){ x++; for(int j = 0 ; j <= x ; j++) m1.mat[i][j] = 1; } x = -1; for(int i = 26 ; i < N ; i++){ x++; for(int j = x ; j < N ; j++) m1.mat[i][j] = 1; } // m2 memset(m2.mat , 0 , sizeof(m2.mat)); x = 24; for(int i = 0 ; i < 26 ; i++){ x++; for(int j = 0 ; j <= x ; j++ ) m2.mat[i][j] = 1; } x = 0; for(int i = 26 ; i < N ; i++){ x++; for(int j = x ; j < N ; j++) m2.mat[i][j] = 1; }}int64 Pow(Matrix m , int n){ Matrix ans; n--; memset(ans.mat , 0 , sizeof(ans.mat)); for(int i = 0 ; i < N ; i++) ans.mat[i][i] = 1; while(n){ if(n&1) ans = ans*m; n >>= 1; m = m*m; } int64 sum = 0; for(int i = 0 ; i < N ; i++){ for(int j = 0 ; j < N ; j++){ sum += ans.mat[i][j]; sum %= MOD; } } return sum;}void solve(Matrix &m1 , Matrix &m2 , int n){ int64 x = Pow(m1 , n); int64 y = Pow(m2 , n); printf("%lld\n" , (x-y+MOD)%MOD);}int main(){ int cas , n; Matrix m1 , m2; init(m1 , m2); scanf("%d" , &cas); while(cas--){ scanf("%d" , &n); solve(m1 , m2 , n); } return 0;}
- HDU 3658 How many words
- hdu 3658 How many words
- hdu 3658 How many words dp+矩阵连乘
- HDU - 3658 How many words 矩阵快速幂
- HDU 3658 How many words (矩阵快速幂&递推)
- hdu 3658 HDU 3658 How many words(矩阵快速幂)
- 杭电acm 3658How many words
- hdu3658 How many words
- hdu 2609 How Many
- hdu How Many Tables
- hdu 2609 How many
- HDU How Many Tables
- hdu 2609 How many
- HDU How Many Tables
- hdu how many ways
- HDU How many days?
- hdu how many tables
- HDU 2609 How Many
- 在Linux内核中单独编译KO文件
- 浮点一般都是不可重入的,printf()经常有重入和性能上的问题
- LoadRunner 用户行为模拟(一)
- how to install dictionary in ubuntu.
- 多校前四版刷题
- hdu 3658 How many words
- ffmpeg 新老接口问题及对照集锦
- 函数模板与类模板
- js分页插件:首页,上页,快退,页码,快进,下页末页,跳转,使用模板,灵活自由
- how to set weblogic password length
- 100个常用vim命令(转)
- mkdir() 和 mkdirs() 这两个方法有什么区别
- js OO写的一个键盘字母游戏
- 【转】获取网页表单信息