【KMP找最小循环节长度】ZOJ 3785

题目：

What day is that day?

           

---

           

                 Time Limit: 2 Seconds                                     Memory Limit: 65536 KB                           

           

---

           

It's Saturday today, what day is it after 1¹ + 2² + 3³ + ... + N^N 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.

3、思路：

打表i的i次方%7后的值，就可以发现42一循环

然后用公式计算出结果

int ans=(n/42%7*(s[42]%7)%7+s[n%42]%7)%7;

最后输出

str[ans]

即可

思路：

这道题有两种做法。

第一种是暴力打表，然后去看看各个数模7的结果，最终找到最小循环节长度为 42.

#include <cstdio>

const char day[10][10] = {"Saturday", "Sunday", "Monday", "Tuesday", "Wednesday", "Thursday", "Friday"};
int s[300];

int work(int n)
{
    int sum = 1;
    for(int i = 1; i <= n; i++){
        sum = sum * n;
        sum %= 7;
    }
    return sum;
}

void init()
{
    s[0] = 0;
    for(int i = 1; i <= 294; i++){
        s[i] = s[i-1] + work(i);
        s[i] %= 7;
    }
}

int main()
{
    int T;
    int n;
    init();
    scanf("%d", &T);
    while(T--){
        scanf("%d", &n);
        n %= 294;
        printf("%s\n", day[ s[n]]);
    }
    return 0;
}

第二种是比较厉害的——通过KMP中的next数组来找最小循环节长度。

（via https://www.cnblogs.com/kevince/p/3887827.html）

不可否认的，对于周期较短的一组数字这样找周期并不难，可是如果周期大到数百数千甚至数万时，靠这种方法找周期恐怕是杯水车薪。

当时我就迷茫在了这一长串的数字中不知所措，猛然想起前不久看过的KMP中next数组的性质，当即想到了用KMP求最小重复子串长度的方法，于是脑洞大开……

该性质为：令j=leni-next[i],如果i%j==0且i/j>1，j就是Pi的最小循环节( Pi表示文本串的前i个字符，leni表示该字符串的长度，一般表示为leni = i + 1)

对此，我通过画图理解了这个性质：

代码：

#include<iostream>
#include<bits/stdc++.h>
using namespace std;

long long x[1000];
long long cal(int x)
{
	long long r=1;
	for(int i=0;i<x;i++) //x次 
	{
		r*=x%7;
		r%=7; 
	}
	return r;
}
int main()
{
	for(int i=0;i<1000;i++)
	{
		x[i]=cal(i+1);
	}
	long long next[1005];
	memset(next,0,sizeof(next));  //必须要memset！！！ 在main函数里定义的数组不会给你默认0的！！！！
	next[0] = 0;
	int len=1000;
	
	/*注：int next[]为next数组，int arr[]为要找规律的数组，len为数组长度*/
	//这里的kmp写法和我的模板不一样，试过我的模板发现结果出不来，反正Kmp我也没理解透彻过，就把这个也当模板记着吧。 
	for(int i = 1, q = 0; i < len; i++)
	{
    	while(q > 0 && x[i] != x[q])  
        	q = next[q-1];
    	if (x[i] == x[q])
        	q++;
    	next[i] = q;
    	if (q != 0 && (i + 1) % (i + 1 - q) == 0) //这里注意0~x[i]长度为i+1，且q的概念也是长度！ 
		{
    	    printf("%d\n", i+1-q); //输出最小循环节长度j 
        	break;
    	}
	}
	return 0;
}   //运行结果：:42