Happy Necklace

Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 131072/131072 K (Java/Others)
Total Submission(s): 1146 Accepted Submission(s): 491

Problem Description

Little Q wants to buy a necklace for his girlfriend. Necklaces are single strings composed of multiple red and blue beads.
Little Q desperately wants to impress his girlfriend, he knows that she will like the necklace only if for every prime length continuous subsequence in the necklace, the number of red beads is not less than the number of blue beads.
Now Little Q wants to buy a necklace with exactly n beads. He wants to know the number of different necklaces that can make his girlfriend happy. Please write a program to help Little Q. Since the answer may be very large, please print the answer modulo 109+7.
Note: The necklace is a single string, {not a circle}.

Input

The first line of the input contains an integer T(1≤T≤10000), denoting the number of test cases.
For each test case, there is a single line containing an integer n(2≤n≤1018), denoting the number of beads on the necklace.

Output

For each test case, print a single line containing a single integer, denoting the answer modulo 109+7.

Sample Input

223

Sample Output

Source

2017中国大学生程序设计竞赛 - 女生专场

题目描述：你有一个有n个颜色为红色或蓝色的珠子的项链，你可以将项链截断成长度为素数的珠子，如果截取出的一段红色的珠子的个数大于蓝色的珠子，则成为good，问一共有多少种good的可能性。

比赛的时候，当拿到这道题的时候，完全一脸懵逼。打了好一会才意识到可以打表，故打只，可以发现

n=2，res=3；n=3，res=4；n=4，res=6；n=5，res=9；

故有递推式 An=An-1+An-3；

赛后问了下师兄，发现正确的推断方法：

如果用a表示红色，用b表示蓝色。题意明显可以看出只需要管长度2和3的连续序列是否符合！

如果以b结尾，那么下一个必须是a，或者加个aab就可以了！

所以同理就可以推出递推式An=An-1+An-3;

看一下题目的数据范围，n最大1e18，因此常规O(n)的递推显然不可行，因此直接上O(logn)的矩阵快速幂。

发现递推式是四阶的递推式，故所得的常数矩阵应该是四维的。之后只需带入矩阵快速幂模板即可。

#include <bits/stdc++.h>
using namespace std;
const int mod=1e9+7;
typedef long long ll;
struct martix{ll mo[4][4];martix(){memset(mo,0,sizeof(mo));}
};
martix mul(martix a,martix b){martix c;for(int i=0;i<4;i++){for(int j=0;j<4;j++){for(int k=0;k<4;k++){c.mo[i][j]=(c.mo[i][j]+a.mo[i][k]*b.mo[k][j])%mod;}}}return c;
}
ll powmod(martix a,ll n){martix T;for(int i=0;i<4;i++){T.mo[i][i]=1;}while(n){if(n&1) T=mul(a,T);n>>=1;a=mul(a,a);}return (T.mo[0][0]*6+T.mo[0][1]*4+T.mo[0][2]*3)%mod;
}
int main()
{ll t;cin>>t;while(t--){ll n;cin>>n;int ans[5]={0,0,3,4,6};if(n<=4){cout<<ans[n]%mod<<endl;}else{martix q;q.mo[0][0]=q.mo[0][2]=q.mo[3][2]=q.mo[1][0]=q.mo[2][1]=q.mo[3][2]=1;cout<<powmod(q,n-4)%mod<<endl;}}return 0;
}

转载于:https://www.cnblogs.com/Chen-Jr/p/11007315.html

HDU 6030（矩阵快速幂+规律）相关推荐

HDU 5411(矩阵快速幂)
本题目意思: 给定n个点,给定一个n*n的转移矩阵,要求求出所有长度不大于M的序列个数,(M<=1E5 , n <= 50) 分析: 记 d[i][j] 生成至多长度为i的且以j开头的所有 ...
hdu 1757(矩阵快速幂)
题意:公式如下 If x < 10 f(x) = x. If x >= 10 f(x) = a0 * f(x-1) + a1 * f(x-2) + a2 * f(x-3) + -- + a ...
【HDU - 4990】 Reading comprehension （构造+矩阵快速幂）
题干: Read the program below carefully then answer the question. #pragma comment(linker, "/STACK ...
HDU 6185 Covering 矩阵快速幂递推
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=6185 题目描述: 一个4*n的矩形, 你用1*2的矩形覆盖有多少种方案, n <= 1e18 ...
hdu 5451 Best Solver 矩阵循环群+矩阵快速幂
http://acm.hdu.edu.cn/showproblem.php?pid=5451 题意:给定x 求解思路: 由斐波那契数列的两种表示方法, 之后可以转化为线性表示 F[n] = ...
HDU - 4990 Reading comprehension(矩阵快速幂，水题)
题目链接:点击查看题目大意:给出一段程序,进行优化后提交题目分析:其实就是找规律,大水题一个,偶尔也是需要做做水题找找自信(逃) 先将题目中的程序拿下来,跑上100项,然后拿到oeis里找一下规律 ...
HDU - 6185 Covering(暴搜+递推+矩阵快速幂/杜教BM)
题目链接:点击查看题目大意:规定宽度为4,给定长度为n,求用1*2和2*1的瓷砖,将其完全铺满能有多少种方法. 分析:自从学会了矩阵快速幂之后,看到1e18的数据量都会下意识的往递推上面想,但是以前 ...
HDU 1757 A Simple Math Problem（矩阵快速幂）
题目链接题意 :给你m和k, 让你求f(k)%m.如果k<10,f(k) = k,否则 f(k) = a0 * f(k-1) + a1 * f(k-2) + a2 * f(k-3) + -- ...
HDU 2276 Kiki Little Kiki 2 （位运算+矩阵快速幂）
HDU 2276 Kiki & Little Kiki 2 (位运算+矩阵快速幂) ACM 题目地址:HDU 2276 Kiki & Little Kiki 2 题意: 一排灯,开关 ...
hdu 5411 CRB and Puzzle（矩阵快速幂）
题目链接: http://acm.hdu.edu.cn/showproblem.php?pid=5411 解题思路: 题目大意: 给定n个点常数m 下面n行第i行第一个数字表示i点的出边数,后面给出 ...

HDU 6030（矩阵快速幂+规律）

Happy Necklace

HDU 6030（矩阵快速幂+规律）相关推荐

最新文章

热门文章