首页 > 学院 > 开发设计 > 正文

ZOJ 3785 What day is that today? (费马小定理)

2019-11-08 19:38:09
字体:
来源:转载
供稿:网友

题目大意

今天是周六,问1^1 + 2^2 + 3^3 + … + N^N天后是周几?(1<=N<=1000000000)

解题思路

第一眼看到题目天真的以为是快速幂取模的简单应用,在飞速敲完代码看到火红的 TLE 之后意识到自己想多了,由于N的范围达到了10^9,所以快速幂取模的方法行不通。

接下来就是使用万能打表A题法,但是由于本题中的循环数太大(294)导致在眼睛看瞎之前很有可能找不到循环点(反正我是没有找到),既然打表也不靠谱,倒不如用数学的方式推导出公式。既靠谱又不伤眼睛。

首先我们关注公式中的单项 N^N % 7 的值, 当 N 的值小于7时,毫无疑问底数 N 有6种情况(1-6),当 N 的值大于7时,可以得到 N=7 * k + b, 则

(7 * k + b)^N % 7 = (7 * K+B) % 7 * (7 * K + B)^ (N-1) % 7=b * (7 * K + B) ^ (N-1) % 7 ,经过递推最终可以得到

(7 * K + B)^N % 7 = B ^ N % 7** (0<=B<=6),根据费马小定理 在 gcd (A,P)=1 的条件下存在 A ^ (P-1)=1(mod P)。

在本题中 ,N % 7 !=0 时,恒有 B^6 % 7 = 1 ,则 B^(6*K2+B2)=B ^ B2(0<=B2<=5)。当 N % 7==0 时 , 毫无疑问 N^N%7==0 恒成立, 则 我们可以得出一个结论 N^N % 7的数值 共有 B * B2 =6 * 7=42种情况。但是此时求出的只是单项的循环结,所以将结果再乘前缀和的循环结7即可得到最终的循环结 294 。


代码

#include <cstdio>#include <cstring>#include <algorithm>using namespace std;typedef long long int ll;int ans[1001];ll quick(ll a,ll b,ll c){ ll ans=1; while(b) { a=a%c; ans=ans%c; if(b%2==1) ans=(ans*a)%c; b/=2; a=(a*a)%c; } return ans%c;}void init(){ ans[1]=1; for(int i=2;i<=294;i++){ ans[i]=(ans[i-1]+quick(i,i,7))%7; }}int main(){ ll t,sum; int n;init(); scanf("%lld",&t); while(t--) { sum=0; scanf("%d",&n); sum=ans[n%294]; if(sum==0) PRintf("Saturday/n"); else if(sum==1) printf("Sunday/n"); else if(sum==2) printf("Monday/n"); else if(sum==3) printf("Tuesday/n"); else if(sum==4) printf("Wednesday/n"); else if(sum==5) printf("Thursday/n"); else printf("Friday/n"); } return 0;}
发表评论 共有条评论
用户名: 密码:
验证码: 匿名发表