GYM-101667-H Rock Paper Scissors FFT

GYM-101667-H
规则就是剪子包袱锤，机器按顺序出，你也按顺序出，问你最多能拿多少分！你可以选择从那个位置开始出。

乍一看，字符串问题，然后你就做不出来了。

实际上是FFT，FFT卷积这个优美的性质。

每次把我的出手序列反转，然后做卷积的时候就可以一一对应了！

一共做三次FFT，每次把我能赢对面的赋值为1，比如我是剪刀，那我就把机器的包袱赋值1，我的剪刀赋值1，这样一遍FFT一下就得到了答案。

举个例子，机器序列是
$$a1,a2,a3,a4,a5,a6$$
我的序列是
$$b3,b2,b1（已经逆序）$$
然后a和b做卷积，c3=a1b1+a2b2+a3b3，这就是起点为1的答案

c4=a2b1+a3b2+a4b3 这就是起点为2的，剪刀的答案！

然后每次FFT结果存下来，然后每次我取当前ans和同样起点的三次FFT结果的和取最大值。

ans = max(ans,fft1[pos]+fft2[pos]+fft3[pos]);

因为，每次FFT的结果都只是一种出法下的结果，一共有三种出法。

模板是kuangbin的FFT模板

#include <stdio.h>
#include <string.h>
#include <iostream>
#include <algorithm>
#include <math.h>
using namespace std;
typedef long long ll;
const double PI = acos(-1.0);
//复数结构体
struct complex
{
    double r,i;
    complex(double _r = 0.0,double _i = 0.0)
    {
        r = _r; i = _i;
    }
    complex operator +(const complex &b)
    {
        return complex(r+b.r,i+b.i);
    }
    complex operator -(const complex &b)
    {
        return complex(r-b.r,i-b.i);
    }
    complex operator *(const complex &b)
    {
        return complex(r*b.r-i*b.i,r*b.i+i*b.r);
    }
};
/*
 * 进行FFT和IFFT前的反转变换。
 * 位置i和 （i二进制反转后位置）互换
 * len必须去2的幂
 */
void change(complex y[],int len)
{
    int i,j,k;
    for(i = 1, j = len/2;i < len-1; i++)
    {
        if(i < j)swap(y[i],y[j]);
        //交换互为小标反转的元素，i<j保证交换一次
        //i做正常的+1，j左反转类型的+1,始终保持i和j是反转的
        k = len/2;
        while( j >= k)
        {
            j -= k;
            k /= 2;
        }
        if(j < k) j += k;
    }
}
/*
 * 做FFT
 * len必须为2^k形式，
 * on==1时是DFT，on==-1时是IDFT
 */
void fft(complex y[],int len,int on)
{
    change(y,len);
    for(int h = 2; h <= len; h <<= 1)
    {
        complex wn(cos(-on*2*PI/h),sin(-on*2*PI/h));
        for(int j = 0;j < len;j+=h)
        {
            complex w(1,0);
            for(int k = j;k < j+h/2;k++)
            {
                complex u = y[k];
                complex t = w*y[k+h/2];
                y[k] = u+t;
                y[k+h/2] = u-t;
                w = w*wn;
            }
        }
    }
    if(on == -1)
        for(int i = 0;i < len;i++)
            y[i].r /= len;
}
const int MAXN = 300010;
char s1[MAXN],s2[MAXN];
complex x1[MAXN],x2[MAXN];
int len1,len2;
ll num1[MAXN/2],num2[MAXN/2];
ll sum[MAXN];
ll ans1[MAXN],ans2[MAXN],ans3[MAXN];
int main()
{
    scanf("%d%d",&len1,&len2);
    getchar();
    scanf("%s%s",s1,s2);
    int len = 1;
    while(len < (len1+len2) ) len<<=1;//补长
    for(int i = 0;i < (len2+1)/2;i ++)	swap(s2[i],s2[len2-i-1]);
 
    for(int i = 0;i < len1;i ++)
	{
		if(s1[i] == 'S')	x1[i] = complex(1,0);
		else	x1[i] = complex(0,0);
	}
	for(int i = 0;i < len2;i ++)
	{
		if(s2[i] == 'R')	x2[i] = complex(1,0);
		else	x2[i] = complex(0,0);
	}
	for(int i = len1;i < len;i++)
        x1[i] = complex(0,0);//补0
    for(int i = len2;i < len;i++)
        x2[i] = complex(0,0);//补0
    
    //求DFT
    fft(x1,len,1);
    fft(x2,len,1);
 
    for(int i = 0;i < len;i++)//计算F(x)
        x1[i] = x1[i]*x2[i];
    //IDFT
    fft(x1,len,-1);
    for(int i = 0;i < len;i++)
    {
        sum[i] = (ll)(x1[i].r+0.5);//计算答案
//    	printf("i and sum: %d %lld\n",i,sum[i]);
    }
    for(int i = len2-1;i < len1+len2;i ++)
    {
    	ans1[i-len2+1] = sum[i];
    //	printf("i and sum and sub : %d %lld %d\n",i,sum[i],i-len2+1);
    }
 
    for(int i = 0;i < len1;i ++)
	{
		if(s1[i] == 'P')	x1[i] = complex(1,0);
		else	x1[i] = complex(0,0);
	}
	for(int i = 0;i < len2;i ++)
	{
		if(s2[i] == 'S')	x2[i] = complex(1,0);
		else	x2[i] = complex(0,0);
	}
	for(int i = len1;i < len;i++)
        x1[i] = complex(0,0);//补0
    for(int i = len2;i < len;i++)
        x2[i] = complex(0,0);//补0
    
    //求DFT
    fft(x1,len,1);
    fft(x2,len,1);
 
    for(int i = 0;i < len;i++)//计算F(x)
        x1[i] = x1[i]*x2[i];
    //IDFT
    fft(x1,len,-1);
    for(int i = 0;i < len;i++)
        sum[i] = (ll)(x1[i].r+0.5);//计算答案
    for(int i = len2-1;i < len1+len2;i ++)
    {
    	ans2[i-len2+1] = sum[i];
    }
 
    for(int i = 0;i < len1;i ++)
	{
		if(s1[i] == 'R')	x1[i] = complex(1,0);
		else	x1[i] = complex(0,0);
	}
	for(int i = 0;i < len2;i ++)
	{
		if(s2[i] == 'P')	x2[i] = complex(1,0);
		else	x2[i] = complex(0,0);
	}
	for(int i = len1;i < len;i++)
        x1[i] = complex(0,0);//补0
    for(int i = len2;i < len;i++)
        x2[i] = complex(0,0);//补0
    
    //求DFT
    fft(x1,len,1);
    fft(x2,len,1);
 
    for(int i = 0;i < len;i++)//计算F(x)
        x1[i] = x1[i]*x2[i];
    //IDFT
    fft(x1,len,-1);
    for(int i = 0;i < len;i++)
        sum[i] = (ll)(x1[i].r+0.5);//计算答案
    for(int i = len2-1;i < len1+len2;i ++)
    {
    	ans3[i-len2+1] = sum[i];
    }
 
    ll ans = 0;
    for(int i = 0;i < len1;i ++)
    {
    	ans = max(ans,ans1[i] + ans2[i] + ans3[i]);
    }
 
    printf("%lld\n",ans);
    return 0;
}