HDU 3336 Count the string (kmp+dp)

Count the string

Problem Description
It is well known that AekdyCoin is good at string problems as well as number theory problems. When given a string s, we can write down all the non-empty prefixes of this string. For example:
s: "abab"
The prefixes are: "a", "ab", "aba", "abab"
For each prefix, we can count the times it matches in s. So we can see that prefix "a" matches twice, "ab" matches twice too, "aba" matches once, and "abab" matches once. Now you are asked to calculate the sum of the match times for all the prefixes. For "abab", it is 2 + 2 + 1 + 1 = 6.
The answer may be very large, so output the answer mod 10007.
 
Input
The first line is a single integer T, indicating the number of test cases.
For each case, the first line is an integer n (1 <= n <= 200000), which is the length of string s. A line follows giving the string s. The characters in the strings are all lower-case letters.
 
Output
For each case, output only one number: the sum of the match times for all the prefixes of s mod 10007.
 
Sample Input
1 4 abab
 
Sample Output
6
 
Author
foreverlin@HNU
 
Source
 
题意: 求字符串的所有前缀在字符串中出现的次数和
 
dp+求kmp中的next数组 
状态转移方程:dp[i]+=dp[next[i]]+1,ans+=dp[i]
 
 1 #include <cctype>
 2 #include <cstdio>
 3 
 4 const int mod=10007;
 5 const int MAXN=10000010;
 6 
 7 int n,ans,T;
 8 
 9 int next[MAXN],dp[MAXN];
10 
11 char s[MAXN];
12 
13 inline void read(int&x) {
14     int f=1;register char c=getchar();
15     for(x=0;!isdigit(c);c=='-'&&(f=-1),c=getchar());
16     for(;isdigit(c);x=x*10+c-48,c=getchar());
17     x=x*f;
18 }
19 
20 inline void KMP_next() {
21     for(int i=1;i<=n;++i) next[i]=0;
22     next[1]=0;  
23     int j=0;  
24     for(int i=2;i<=n;i++){  
25         while(j>0&&s[j+1]!=s[i]) j=next[j];  
26         if(s[j+1]==s[i]) j+=1;  
27         next[i]=j;  
28     }  
29 }
30 
31 int hh() {
32     freopen("ha.in","r",stdin);
33     freopen("ha.out","w",stdout);
34     read(T);
35     while(T--) {
36         read(n);
37         scanf("%s",s+1);
38         KMP_next();
39         ans=0;
40         for(int i=1;i<=n;++i) {
41             dp[i]=(dp[next[i]]+1)%mod;
42             ans=(ans+dp[i])%mod;
43         }
44         printf("%d
",ans);
45     }
46     fclose(stdin);
47     fclose(stdout);
48     return 0;
49 }
50 
51 int sb=hh();
52 int main(int argc,char**argv) {;}
代码
原文地址:https://www.cnblogs.com/whistle13326/p/7605959.html