BM模板

我之前可能学了假的BM

namespace BM{
	vector<int> h[N];
	int cnt,fail[N],d[N],mx;
	vector<int> work(int n,int *a) {
		h[cnt=mx=0].clear();
		fo(i,1,n) {
			int now=-a[i];
			for(unsigned j=0;j<h[cnt].size();j++) (now+=(ll)a[i-j-1]*h[cnt][j]%Mo)%=Mo;
			d[i]=now;if (!now) continue;
			fail[cnt]=i;
			if (!cnt) {
				h[++cnt].clear();
				h[cnt].resize(i);
				continue;
			}
			vector<int> r;r.resize(i-fail[mx]-1);
			int mul=-(ll)now*pwr(d[fail[mx]],Mo-2)%Mo;
			r.pb(-mul);for(unsigned j=0;j<h[mx].size();j++) r.pb((ll)h[mx][j]*mul%Mo);
			if (r.size()<h[cnt].size()) r.resize(h[cnt].size());
			for(unsigned j=0;j<h[cnt].size();j++) (r[j]+=h[cnt][j])%=Mo;
			if (i-fail[mx]+h[mx].size()>=h[cnt].size()) mx=cnt;
			h[++cnt]=r;
		}
		return h[cnt];
	}
}