HDU - 6153 A Secret 擴展kmp

題目鏈接:https://vjudge.net/problem/HDU-6153

題意:對於s2的每一個後綴,假設長度爲l,在s1出現的次數爲k,求l*k的和

題解:我們把兩個串都倒過來,變爲s1,s2,那麼問題就變爲,對於s2的前綴匹配s1的每個位置的後綴的總匹配長度和,舉個栗子,s1: aaa s2: aa  那麼在s1的位置1,能匹配的長度有1 2,因此我們擴展KMP求出最長長度l,把(1+l)*l/2累加起來即可

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 1000100;
const ll mod = 1e9 + 7;
int nex[N], ex[N];
void get_nex(char *s) {
    int i = 0, j, pos, len = strlen(s);
    nex[0] = len;
    while(s[i] == s[i + 1] && i + 1 < len) i++;
    nex[1] = i;
    pos = 1;
    for(i = 2; i < len; i++) {
        if(nex[i - pos] + i < nex[pos] + pos) 
            nex[i] = nex[i - pos];
        else {
            j = nex[pos] + pos - i;
            if(j < 0) j = 0;
            while(i + j < len && s[j] == s[j + i]) j++;
            nex[i] = j;
            pos = i;
        }
    }
}
void ex_kmp(char *s, char *t) {
    int i = 0, j, pos, l1 = strlen(s), l2 = strlen(t);
    get_nex(s);
    while(s[i] == t[i] && i < l1 && i < l2) i++;
    ex[0] = i;
    pos =  0;
    for(int i = 1; i < l2; i++) {
        if(nex[i - pos] + i < ex[pos] + pos)
            ex[i] = nex[i - pos];
        else {
            j = ex[pos] + pos - i;
            if(j < 0) j = 0;
            while(i + j < l2 && j < l1 && s[j] == t[i + j]) j++;
            ex[i] = j;
            pos = i;
        }
    }
}
char cnt[N];
int main() {
	char s[N], t[N];
	int ls, lt;
	int pos, len;
	int T;
	ll ans;
	scanf("%d", &T);
	while(T--) {
		scanf("%s", cnt);
		len = strlen(cnt);
		lt = len;
		for(int i = 0; i < len; i++)
			t[len - i - 1] = cnt[i];
		scanf("%s", cnt);
		len = strlen(cnt);
		for(int i = 0; i < len; i++) 
			s[len - i - 1] = cnt[i];
		get_nex(s);
		ex_kmp(s, t);
		ans = 0;
		for(int i = 0; i < lt; i++) {
		//	cout << ex[i] << endl;
			ans = (ans + 1LL*ex[i] * (ex[i] + 1) / 2) % mod;
		}
			
		printf("%lld\n", ans);
	}
	return 0;
}

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章