bzoj4259: 残缺的字符串

#链接
　　http://www.lydsy.com/JudgeOnline/problem.php?id=4259
#题解
　　这道题，我先是因为精度WA，然后又MLE，最后TLE。
　　不得不删掉原来的递归fft，学了非递归fft。
　　Claris太强了
　　考虑相等就会差等于0，全相等就是所有差的绝对值之和等于0，通配符可以考虑把它变成0然后在式子中乘上两个数的权值。你把第一个串翻过来，这个时候把两个串求下卷积，如果$c_i=0$就证明开头$i-m$是合法的。
#代码

//FFT
#include <cstdio>
#include <algorithm>
#include <cmath>
#define maxn (1048576+100)
using namespace std;
double pi=acos(-1);
int up, n, m, final[maxn], a[maxn], b[maxn], R[maxn];
char s[maxn];
struct cp
{
	double x, y;
	cp(double x, double y):x(x),y(y){}
	cp(){}
}ans[maxn], t1[maxn], t2[maxn];
inline cp operator+(cp a, cp b){return cp(a.x+b.x,a.y+b.y);}
inline cp operator-(cp a, cp b){return cp(a.x-b.x,a.y-b.y);}
inline cp operator*(cp a, cp b){return cp(a.x*b.x-a.y*b.y,a.x*b.y+a.y*b.x);}
void fft(cp *a, int n, int opt)
{
	int i, j, k, p; cp w, wn, x, y;
	for(i=0;i<n;i++)if(i>R[i])swap(a[i],a[R[i]]);
	for(i=1;i<n;i<<=1)
	{
		wn=cp(cos(pi/i),opt*sin(pi/i));
		for(p=i<<1,j=0;j<n;j+=p)
			for(k=0,w=cp(1,0);k<i;k++,w=w*wn)
			{
				x=a[j+k], y=a[j+k+i]*w;
				a[j+k]=x+y, a[j+k+i]=x-y;
			}
	}
	if(opt==-1)for(i=0;i<n;i++)a[i].x/=n;
}
void juan(cp *a, cp *b)
{
	fft(a,up,1), fft(b,up,1);
	for(int i=0;i<=up;i++)b[i]=a[i]*b[i];
	fft(b,up,-1);
}
void init()
{
	int i, L=0;
	scanf("%d%d",&m,&n);
	scanf("%s",s);
	for(i=0;i<m;i++)a[m-1-i]=s[i]=='*'?0:s[i]-'a'+1;
	scanf("%s",s);
	for(i=0;i<n;i++)b[i]=s[i]=='*'?0:s[i]-'a'+1;
	n--, m--;
	for(up=1;up<=n+m;up<<=1)L++;
	for(i=0;i<up;i++)R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
}
void work()
{
	int i, cnt=0;
	for(i=0;i<=up;i++)t1[i]=cp(a[i],0)*cp(a[i],0)*cp(a[i],0);
	for(i=0;i<=up;i++)t2[i]=cp(b[i],0);
	juan(t1,t2);for(i=0;i<=up;i++)ans[i]=t2[i];
	for(i=0;i<=up;i++)t1[i]=cp(a[i],0);
	for(i=0;i<=up;i++)t2[i]=cp(b[i],0)*cp(b[i],0)*cp(b[i],0);
	juan(t1,t2);for(i=0;i<=up;i++)ans[i]=t2[i]+ans[i];
	for(i=0;i<=up;i++)t1[i]=cp(a[i],0)*cp(a[i],0);
	for(i=0;i<=up;i++)t2[i]=cp(b[i],0)*cp(b[i],0);
	juan(t1,t2);for(i=0;i<=up;i++)ans[i]=ans[i]-cp(2,0)*t2[i];
	for(i=m;i<=n;i++)if((final[i]=(int)(ans[i].x+0.5))==0)cnt++;
	printf("%d\n",cnt);
	for(i=m;i<=n;i++)if(!final[i])printf("%d ",i-m+1);
}
int main()
{
	init();
	work();
	return 0;
}