2019暑假正睿集训8.10day7题解及总结

最新推荐文章于 2019-10-15 18:00:16 发布

Trilarflagz

最新推荐文章于 2019-10-15 18:00:16 发布

阅读量313

点赞数

分类专栏： 2019暑假正睿集训

本文链接：https://blog.csdn.net/flag_z/article/details/99124137

版权

2019暑假正睿集训专栏收录该内容

8 篇文章 0 订阅

订阅专栏

T1 松

50pts~80pts

这道题直接暴力搜索，每次在三个符号中任选一个，按题意合并数。当合并了还剩两个数的时候，就判断它们选哪种符号最终能为1，找到一种ans+1

100pts

优化1：位运算。把输⼊串按照FFT初始化⽅法⼀样bitrev⼀下就可以⽤位运算优化。

优化2：四⽑⼦。当串⻓之后16的时候预处理。 ~~~~ 可能要稍微松⼀下，也可能不⽤。~~

为此我今天学了三个多小时的FFT（理解思路＋代码实现）

其实学会后还~~挺简单的~~（看了无数篇题解后emmm）

就是把两个多项式相乘（即求卷积），直接乘的话复杂度是 $O(n^2)$ 所以把多项式的系数表示转化成点值表示。

举个例子 $f[x]=a_0x^0+a_1x^1+a_2x^2+......a_nx^n$ $a_i$ 是每一项的系数，所以这是多项式f[x]的系数表示法

$x_0,f[x_0])$ $ (x_1,f[x_1])$ … $x_n,f[x_n])$ 则为多项式的点值表示法相当于代入一些x的值，求出对应的y的值，而要代入的x的值的个数显然至少为n个才能确定一个多项式（可类比二元一次方程）。

而对于多项式 $f [x] * g [x]$ 若它们x代入的是相同的数，那么就把算出对于的y值相乘就可以了

这是个总体的思路，还有非常多的细节部分实现，例如引入复数和单位跟，DFT和IDFT，二进制转换，递归和迭代两种方法等等，一时半会真的说不清楚，之后再详细写一篇有关于FFT的题解QWQ。

我写的是更为高效的迭代算法：

附上代码如下

#include<bits/stdc++.h>
using namespace std;
typedef complex<double> cp;
#define N 2097153
const double pie=acos(-1);
int n;
cp a[N],b[N];
int rev[N],ans[N];
char s1[N],s2[N];
int read(){
	int sum=0,f=1;
	char ch=getchar();
	while(ch>'9'||ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0'&&ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
	return sum*f;
}
void init(int k)
{
    int len=1<<k;
	for(int i=0;i<len;i++)
	rev[i]=(rev[i>>1]>>1)|((i&1)<<(k-1));
}
void fft(cp *a,int n,int flag){
    for(int i=0;i<n;i++)
	{
	  if(i<rev[i])swap(a[i],a[rev[i]]);
	}
	for(int h=1;h<n;h*=2)
	{
	cp wn=exp(cp(0,flag*pie/h));
	 for(int j=0;j<n;j+=h*2)
	 {
	  cp w(1,0);
	   for(int k=j;k<j+h;k++)
	   {
	     cp x=a[k];
	     cp y=w*a[k+h];
         a[k]=x+y;
         a[k+h]=x-y;
         w*=wn;
	   }
	 }
	 }
	 if(flag==-1)
	 for(int i=0;i<n;i++)
     a[i]/=n;
}
int main(){
	n=read();
	scanf("%s%s",s1,s2);
    for(int i=0;i<n;i++)a[i]=(double)(s1[n-i-1]-'0');
	for(int i=0;i<n;i++)b[i]=(double)(s2[n-i-1]-'0');
	int k=1,s=2;
 	while((1<<k)<2*n-1)k++,s<<=1;
	init(k);
    fft(a,s,1);
    fft(b,s,1);
    for(int i=0;i<s;i++)
    a[i]*=b[i];
    fft(a,s,-1);
    for(int i=0;i<s;i++)
    {
	ans[i]+=(int)(a[i].real()+0.5);
	ans[i+1]+=ans[i]/10;
	ans[i]%=10;
	}
	while(!ans[s]&&s>-1)s--;
	if(s==-1)printf("0");
	else
	for(int i=s;i>=0;i--)
	printf("%d",ans[i]);
	return 0;
}