hdu 1005 用mod的循环性质

最新推荐文章于 2022-07-09 15:36:58 发布

linleiqin

最新推荐文章于 2022-07-09 15:36:58 发布

阅读量476

点赞数

CC 4.0 BY-SA版权

分类专栏：解题报告文章标签：扩展

本文链接：https://blog.youkuaiyun.com/linleiqin/article/details/5616213

解题报告专栏收录该内容

12 篇文章

订阅专栏

本文介绍了一种高效求解特定递推数列f(n)模7的方法，利用状态压缩和循环特性避免重复计算，适用于给定的大规模n值。通过构建状态转移矩阵，快速定位数列中的任意项。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

http://acm.hdu.edu.cn/showproblem.php?pid=1005
A number sequence is defined as follows:
f(1) = 1, f(2) = 1, f(n) = (A * f(n - 1) + B * f(n - 2)) mod 7.
Given A, B, and n, you are to calculate the value of f(n).
1 <= A, B <= 1000, 1 <= n <= 100,000,000

解:
f(n) = (A * f(n - 1) + B * f(n - 2)) %7
= (A * f(n - 1) %7 + B * f(n - 2) %7) %7
所以对于给定的A和B,可以先打表,找出数列的循环部分. 鸽巢原理知,状态总数不会超过7*7
注意循环节不一定从f(3)开始...
#include <iostream> using namespace std; int a,b,n,x,y,l,h,m[7][7],q[300]; int main() { while(scanf("%d%d%d",&a,&b,&n)!=EOF && (a||b||n)) { memset(m, 0, sizeof(m)); q[1] = q[2] = 1; x = y = 1; l=3; while(m[x][y]==0) { //该状态还未经历过,则扩展 q[l] = (a*x+b*y)%7; m[x][y] = l; y = x; x = q[l++]; } //此时,q[1...h-1]为前面的非循环部分 //q[h...l-1]为循环节 h = m[x][y]; //循环节的起始位置 if(n<h) { printf("%d/n",q[n]); } else { printf("%d/n",q[((n-h)%(l-h))+h]); } } return 0; }