题目描述


NowCoder最近在研究一个数列:

  • F(0) = 7
  • F(1) = 11
  • F(n) = F(n-1) + F(n-2) (n≥2)

    他称之为NowCoder数列。请你帮忙确认一下数列中第n个数是否是3的倍数。

输入描述:

输入包含多组数据。

每组数据包含一个整数n,(0≤n≤1000000)。


输出描述

对应每一组输入有一行输出。

如果F(n)是3的倍数,则输出“Yes”;否则输出“No”。

输入例子:

0

1

2

3

4

5

输出例子:

No

No

Yes

No

No

No

题目分析

这是一个特殊初始条件的Fibonacci sequence.

1.初始条件n的值比较小,可以直接枚举每个F(n)的值,然后输入n查询:

#include <iostream>
#include <cstdio>
using namespace std;
int main (){
int n;
int *a=new int[N];
a[0]=7;
a[1]=11;
for(int i=2;i<N;i++)
a[i]=a[i-2]%3+a[i-1]%3;
while(~scanf("%d", &n)){
cout<<a[n]<<endl;
if(a[n]%3==0)
cout<<"Yes"<<endl;
else
cout<<"No"<<endl;
}
return 0;
}

2.可以寻找到规律:(牛客网友“我要过pat”提供)

这个真牛批,想不到还有这种规律。那Fibonacci数列是否有相似的规律呢?

#include <iostream>
using namespace std; int main (){
int n;
while(cin>>n){
if((n-2)%4==0)
cout<<"Yes"<<endl;
else
cout<<"No"<<endl;
}
return 0;
}

3.用矩阵快速幂来求解,这题n的范围比较小(0≤n≤1000000)所以用第1种方法可以不超时求解,但是,倘若n很大,达到(1<=n<=1000,000,000,000,000,000)这个范围,那么显而易见,必定超时。而且根本开不到那么大的数组。学校的算法课,有另外一题类似的,可以思考一下:

借此思路,用矩阵快速幂计算,稍微修改下就可以了。

 #include <iostream>
#include <cstdio>
#include <cstring>
#define N 3
#define maxn 2 using namespace std; struct Matrix{
long long a[maxn][maxn];
void init(){ //初始化为单位矩阵
memset(a, 0, sizeof(a));
for(int i=0;i<maxn;++i){
a[i][i] = 1;
}
}
}; //矩阵乘法
Matrix mul(Matrix a, Matrix b){
Matrix ans;
for(int i=0;i<maxn;++i){
for(int j=0;j<maxn;++j){
ans.a[i][j] = 0;
for(int k=0;k<maxn;++k){
ans.a[i][j] += a.a[i][k] * b.a[k][j];
ans.a[i][j] %= N;
}
}
}
return ans;
} //矩阵快速幂
Matrix qpow(Matrix a, long long n){
Matrix ans;
ans.init();
while(n){
if(n&1)
ans = mul(ans, a);
a = mul(a, a);
n /= 2;
}
return ans;
} int main (){
long long n; while(~scanf("%lld", &n)){
Matrix a;
a.a[0][0] = 1;
a.a[0][1] = 1;
a.a[1][0] = 1;
a.a[1][1] = 0;
long long s0=7,s1=11;
if (n>1)
{
Matrix ans= qpow(a, n-1);
long long res=ans.a[0][0]*s1+ans.a[0][1]*s0;
// printf("%lld",res);
if(res%3==0)
cout<<"Yes"<<endl;
else
cout<<"No"<<endl;
}
else
cout<<"No"<<endl; }
return 0;
}

测评网站真的恶心,有时候循环输入用"cin>>"就可以,有时候又超时。。。得改"成~scanf("%lld", &n)"。另外一开始没想到对中间结果对3取模,导致数值过大溢出。用到unsigned long long 都不够。。。

PS:检测一个数被3整除的算法

1.检测一个数能否被3整除----位运算

如果所有的偶数位出现1的次数为 even_count, 奇数位出现1的次数为 odd_count,两者只差如果是3的倍数,那么这个数就是3倍数。

2.不用除法和求模运算,判断一个数能否被3整除

现在给出一个数a,假设它能被3整除,结果是b,即a=3*b,那么从二进制乘法运算判断出,b的最低位与a的最低位一定是相同的,从而得到了b的最低位,将这个位左移1位变成次低位,那么a的次低位以上的比特减去这个位后在次低位上的结果一定是b的次低位。以此类推可以求出b的各个比特,如果最后能完成对b的各位的计算,那么a能够被3整除,否则不能被3整除。

那么算法原理是什么呢?

PAT乙级(Basic Level)练习题-NowCoder数列总结的更多相关文章

  1. 牛客网 PAT乙级(Basic Level)练习题 1023 考新郎

    题目描述 过年期间,老家举行了一场盛大的集体婚礼,为了使婚礼进行的丰富一些,司仪临时想出了有一个有意思的节目,叫做“考新郎”,具体的操作是这样的: 1. 首先,给每位新娘打扮得几乎一模一样,并盖上大大 ...

  2. C#版 - PAT乙级(Basic Level)真题 之 1021.个位数统计 - 题解

    版权声明: 本文为博主Bravo Yeung(知乎UserName同名)的原创文章,欲转载请先私信获博主允许,转载时请附上网址 http://blog.csdn.net/lzuacm. C#版 - P ...

  3. C#版 - PAT乙级(Basic Level)真题 之 1024.科学计数法转化为普通数字 - 题解

    版权声明: 本文为博主Bravo Yeung(知乎UserName同名)的原创文章,欲转载请先私信获博主允许,转载时请附上网址 http://blog.csdn.net/lzuacm. PAT Bas ...

  4. 牛客网编程练习之PAT乙级(Basic Level):1033 害死人不偿命的(3n+1)猜想

    3n+1水题.... AC代码: import java.util.Scanner; /** * @author CC11001100 */ public class Main { public st ...

  5. 牛客网编程练习之PAT乙级(Basic Level):1034 写出这个数

    AC代码: import java.util.*; /** * @author CC11001100 */ public class Main { public static void main(St ...

  6. 牛客网编程练习之PAT乙级(Basic Level):1032 选大王

    典型的约瑟夫环问题 AC代码: import java.util.ArrayList; import java.util.List; import java.util.Scanner; /** * @ ...

  7. 牛客网编程练习之PAT乙级(Basic Level):1041 说反话

    直接分隔取反即可 AC代码: import java.util.Scanner; /** * @author CC11001100 */ public class Main { public stat ...

  8. 牛客网PAT乙级(Basic Level)真题-组个最小数 (20)

    组个最小数 (20) 时间限制 1000 ms 内存限制 32768 KB 代码长度限制 100 KB 判断程序 Standard (来自 小小) 题目描述 给定数字0-9各若干个.你可以以任意顺序排 ...

  9. 牛客网PAT乙级(Basic Level)真题-数字分类 (20)

    题目描述 给定一系列正整数,请按要求对数字进行分类,并输出以下5个数字: A1 = 能被5整除的数字中所有偶数的和: A2 = 将被5除后余1的数字按给出顺序进行交错求和,即计算n1-n2+n3-n4 ...

随机推荐

  1. cpu高占用,线程堆栈,jstack,pstack,jmap, kill -3 pid,java(weblogic,tomcat)

    1 ps -mp pid -o THREAD,tid,time 2 printf "%x\n" tid 3 jstack pid |grep tid -A 30

  2. pip install Yellowfin失败的问题

    上面的方法简单有效 https://blog.csdn.net/quqiaoluo5620/article/details/80608474

  3. POJ 1144 Network(tarjan 求割点个数)

    Network Time Limit: 1000MS   Memory Limit: 10000K Total Submissions: 17016   Accepted: 7635 Descript ...

  4. bat输出重定向

    重定向符号主要有:>,>>,<,>&,<&和|,而本文只讨论前五个. 第一节 首先从一个经典问题开始,“1>nul 2>nul”的意思是 ...

  5. Pythoner使用的豆瓣pip源

    主要示例: sudo pip install -i http://pypi.douban.com/simple/  flask-script Flask的扩展: flask-script是一个可以在f ...

  6. Swift使用AlamoFire超时设置和事件处理

    一直在写swift项目,正好碰到服务器部署,请求超时或者请求失败的问题,页面就卡着不动了.顺手解决一下吧 差了些资料,说要设置超时时间 方法一: static let sharedSessionMan ...

  7. 20155237 2016-2017-2 《Java程序设计》第十周学习总结

    20155237 2016-2017-2 <Java程序设计>第十周学习总结 教材学习内容总结 计算机网络,是指分布在不同地理区域的计算机用通信线路互连起来的一个具有强大功能的网络系统.网 ...

  8. Educational Codeforces Round 41 (Rated for Div. 2)(A~D)

    由于之前打过了这场比赛的E题,而后面两道题太难,所以就手速半个多小时A了前4题. 就当练手速吧,不过今天除了C题数组开小了以外都是1A A Tetris 题意的抽象解释可以在Luogu里看一下(话说现 ...

  9. python 回溯法 子集树模板 系列 —— 16、爬楼梯

    问题 某楼梯有n层台阶,每步只能走1级台阶,或2级台阶.从下向上爬楼梯,有多少种爬法? 分析 这个问题之前用分治法解决过.但是,这里我要用回溯法子集树模板解决它. 祭出元素-状态空间分析大法:每一步是 ...

  10. Kubernetes学习之路(二十)之K8S组件运行原理详解总结

    目录 一.看图说K8S 二.K8S的概念和术语 三.K8S集群组件 1.Master组件 2.Node组件 3.核心附件 四.K8S的网络模型 五.Kubernetes的核心对象详解 1.Pod资源对 ...