【ZOJ】3785 What day is that day? ——KMP 暴力打表找规律
转自:http://www.cnblogs.com/kevince/p/3887827.html
首先声明一下,这里的规律指的是循环,即找到最小循环周期。
这么一说大家心里肯定有数了吧,“不就是next数组性质的应用嘛”,没错,正是如此。
在ACM的比赛中有些时候会遇到一些题目,可以或必须通过找出数据的规律来编写代码,这里我们专门来讨论下 如何运用KMP中next数组的性质 来寻找一个长数组中的最小循环周期。
先来看一道题
ZOJ 3785
What day is that day?
Time Limit: 2 Seconds Memory Limit: 65536 KB
It's Saturday today, what day is it after 11 + 22 + 33 + ... + NN days?
Input
There are multiple test cases. The first line of input contains an integer T indicating the number of test cases. For each test case:
There is only one line containing one integer N (1 <= N <= 1000000000).
Output
For each test case, output one string indicating the day of week.
Sample Input
2
1
2
Sample Output
Sunday
Thursday
Hint
A week consists of Sunday, Monday, Tuesday, Wednesday, Thursday, Friday and Saturday.
题目的大意是知道今天是周六,让你求 f = 11 + 22 + 33 + ... + NN 这么多天之后是星期几。
也就是求f % 7对于每个输入的N的值。这题在网上一搜题解,都说是打表找规律,当然这题有两种找法,一是对于每个ii % 7 的值都找规律。
这里我们打表可知 前100个值如下所示
1 4 6 4 3 1 0 1 1 4 2 1 6 0 1 2 5 1 5 1 0 1 4 1 4 4 6 0 1 1 3 2 6 1 0 1 2 2 1 2 6 0 1 4 6 4 3 1 0 1 1 4 2 1 6 0 1 2 5 1 5 1 0 1 4 1 4 4 6 0 1 1 3 2 6 1 0 1 2 2 1 2 6 0 1 4 6 4 3 1 0 1 1 4 2 1 6 0 1 2
有一种找规律的方法是当有数字等于第一个数的时候做个标记,再人工判断是否能够构成一个循环。
不可否认的,对于周期较短的一组数字这样找周期并不难,可是如果周期大到数百数千甚至数万时,靠这种方法找周期恐怕是杯水车薪。
当时我就迷茫在了这一长串的数字中不知所措,猛然想起前不久看过的KMP中next数组的性质,当即想到了用KMP求最小重复子串长度的方法,于是脑洞大开……
该性质为:令j=leni-next[i],如果i%j==0且i/j>1,j就是Pi的最小循环节( Pi表示文本串的前i个字符,leni表示该字符串的长度,一般表示为leni = i + 1)
还有一种找规律的方法是直接对 f % 7 的值进行打表找规律,按照上述方法找到的周期为294,下面要做的就很简单了~
#include <iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
#define N 605
#define M 200
#define ll long long
using namespace std;
int i,j;
int T;
int l;
int a[N],sum[N];
char s[N];
int next[N];
int fun(int n)
{
int re = ;
for(int i=;i<=n;i++){
re*=n;
re%=;
}
return re;
}
void ini()
{
a[]=sum[]=;
a[]=sum[]=;
s[]='';
s[]='';
for(i=;i<=N-;i++){
a[i]=fun(i);
sum[i]+=a[i]+sum[i-];
sum[i]%=;
s[i]=sum[i]+'';
}
s[i]='\n';
//printf("%s\n",s);
// for(i=1;i<=M;i++){
// printf(" %d %d %d\n",i,a[i],sum[i]);
// printf(" %d",sum[i]);
//}
return;
} void ini2()
{
int j = -, i = ;
next[] = -;
while(i < N-)
{
if(j == - || s[i] == s[j])
{ i++;
j++;
next[i] = j;
}
else
{
j = next[j];
}
}
} void ini3()
{
for(int i = ; i <= N-; ++i)
{
int length = i - next[i]; //循环节的长度
if(i != length && i % length == ) //如果有多个循环
{
printf("%d %d\n", i, i / length);
break;
} }
} int main()
{
int x;
int ans;
ini();
//ini2();
// ini3();
//freopen("data.txt","r",stdin);
scanf("%d",&T);
//while(scanf("%d",&n)!=EOF)
while(T--)
{
scanf("%d",&x);
ans=sum[x%];
if(ans==) printf("Saturday\n");
else if(ans==) printf("Sunday\n");
else if(ans==) printf("Monday\n");
else if(ans==) printf("Tuesday\n");
else if(ans==) printf("Wednesday\n");
else if(ans==) printf("Thursday\n");
else if(ans==) printf("Friday\n");
}
return ;
}
后来又练了一次,用了map
4172110 | 2016-03-15 15:26:39 | Accepted | 3785 | C++ | 820 | 280 | njczy2010 |
#include <iostream>
#include <cstdio>
#include <map>
#include <cstring>
#include <algorithm>
#include <cmath>
#include <string> using namespace std; #define N 1005
#define ll long long
#define mod 7 int TT;
int n;
int s[N];
int next[N];
map<int,string> mp; void get_next()
{
int i,j;
i = ;
j = -;
next[] = -;
while(i<N)
{
if(j == - || s[i] == s[j]){
i++;j++;next[i] = j;
}
else{
j = next[j];
}
}
} int quickpow(int x,int n)
{
int re = ;
while(n)
{
if(n&){
re = (re * x) % mod;
}
n /= ;
x = (x * x) % mod;
}
return re;
} void fi()
{
int i;
for(i = ;i < N;i++){
int length = i - next[i]; //循环节的长度
if(i != length && i % length == ) //如果有多个循环
{
printf("%d %d %d\n", i, length, i / length);
break;
}
}
} void ini()
{
int i;
s[] = ;
for(i = ;i < N;i++){
s[i] = (s[i - ] + quickpow(i,i) )%mod;
}
mp[] = "Saturday";
mp[] = "Sunday";
mp[] = "Monday";
mp[] = "Tuesday";
mp[] = "Wednesday";
mp[] = "Thursday";
mp[] = "Friday";
//get_next();
//fi();
//
} int main()
{
//freopen("in.txt","r",stdin);
//freopen("out.txt","w",stdout);
ini();
scanf("%d",&TT);
while(TT--){
//while(scanf("%d",&n)!=EOF){
scanf("%d",&n);
cout << mp[ s[ n% ] ] << endl;
}
return ;
}
【ZOJ】3785 What day is that day? ——KMP 暴力打表找规律的更多相关文章
- 【ZOJ】3785 What day is that day? ——浅谈KMP在ACM竞赛中的暴力打表找规律中的应用
转载请声明出处:http://www.cnblogs.com/kevince/p/3887827.html ——By Kevince 首先声明一下,这里的规律指的是循环,即找到最小循环周期. 这 ...
- ZOJ 3622 Magic Number 打表找规律
A - Magic Number Time Limit:2000MS Memory Limit:32768KB 64bit IO Format:%lld & %llu Subm ...
- zoj 3785 What day is that day? (打表找规律)
题目 思路:比赛的时候有想过找循环节,但是,打表打错了. 后来,看着过了挺多人,就急了, 看了一下别人的时间 耗时都挺长的,就以为不是找规律, 没想到真是找规律,不过,这个题的数据可能挺大的. AC代 ...
- zoj 3629 Treasure Hunt IV 打表找规律
H - Treasure Hunt IV Time Limit:2000MS Memory Limit:65536KB 64bit IO Format:%lld & %llu ...
- 【ARC077F】SS kmp+打表找规律
Description 如果某个串可以由两个一样的串前后连接得到,我们就称之为"偶串".比如说"xyzxyz"和"aaaaaa"是偶串, ...
- ZOJ 3785 What day is that day?(今天是星期几?)
Description 题目描述 It's Saturday today, what day is it after 11 + 22 + 33 + ... + NN days? 今天是星期六,11 + ...
- ZOJ 3587 Marlon's String 扩展KMP
链接:http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemCode=3587 题意:给出两个字符串S和T.S,T<=100000.拿出 ...
- zoj 3785 What day is that day?
http://acm.zju.edu.cn/onlinejudge/showProblem.do?problemId=5272 打表找规律. #include <cstdio> #incl ...
- ZOJ 3785 What day is that day?(数论:费马小定理)
What day is that day? Time Limit: 2 Seconds Memory Limit: 65536 KB It's Saturday today, what da ...
随机推荐
- BZOJ 3876:支线剧情(有下界最小费用最大流)
3876: [Ahoi2014]支线剧情 Description [故事背景]宅男JYY非常喜欢玩RPG游戏,比如仙剑,轩辕剑等等.不过JYY喜欢的并不是战斗场景,而是类似电视剧一般的充满恩怨情仇的剧 ...
- [BZOJ5417] [NOI2018]你的名字
Description 小A 被选为了ION2018 的出题人,他精心准备了一道质量十分高的题目,且已经把除了题目命名以外的工作都做好了. 由于ION 已经举办了很多届,所以在题目命名上也是有规定的, ...
- 第五届华中区程序设计邀请赛暨武汉大学第十四届校赛 网络预选赛 A
Problem 1603 - Minimum Sum Time Limit: 2000MS Memory Limit: 65536KB Total Submit: 564 Accepted: ...
- Codeforces Round #350 (Div. 2) A
A. Holidays time limit per test 1 second memory limit per test 256 megabytes input standard input ou ...
- 洛谷P1546 最短网络 Agri-Net
P1546 最短网络 Agri-Net 526通过 959提交 题目提供者JOHNKRAM 标签图论贪心USACO 难度普及/提高- 提交该题 讨论 题解 记录 最新讨论 50分C++代码,求解 请指 ...
- 记录一发wm_concat()函数排序的问题
需求:需要将列转行之后的工序按照待执行工序号排序,如果一样按工序号排 解决方法如下: select part_no, max(ywggx) ywggx from(select mt.part_no , ...
- 获取html元素内容
html: <!DOCTYPE ><html> <head> <meta http-equiv="Content-Type" conten ...
- sperman系数
https://baike.baidu.com/item/spearman%E7%9B%B8%E5%85%B3%E7%B3%BB%E6%95%B0/7977847?fr=aladdin https:/ ...
- javaScript获取文档中所有元素节点的个数
HTML+JS 代码: <!DOCTYPE html> <html lang="en"> <head> <meta charset=&qu ...
- php模板引擎smarty
一. smarty的特点 速度:相对于其他模板引擎,速度较快 编译型:在下次访问模板时直接访问编译文件,不再进行模板重新编译 缓存技术:可以将用户最终看到的HTML文件缓存成一个静态HTML 插件技术 ...