zoukankan      html  css  js  c++  java
  • BZOJ 2194 FFT

    思路:
    题目中给的是差值一定的

    那么就把b数组倒一下 不就变成 i+j=k(k为定值的了嘛)
    嗯 然后搞个FFT

    //By SiriusRen
    #include <cstdio>
    #include <complex>
    using namespace std;
    typedef complex<double> cplxd;
    const int N=270000;const double pi=acos(-1);
    int n,m,R[N],L;
    cplxd a[N],b[N];
    char ch[N];
    void FFT(cplxd *a,int f){
        for(int i=0;i<n;i++)if(i<R[i])swap(a[i],a[R[i]]);
        for(int i=1;i<n;i<<=1){
            cplxd wn(cos(pi/i),f*sin(pi/i));
            for(int j=0;j<n;j+=(i<<1)){
                cplxd w(1,0);
                for(int k=0;k<i;k++,w*=wn){
                    cplxd x=a[j+k],y=w*a[j+k+i];
                    a[j+k]=x+y,a[j+k+i]=x-y; 
                }
            }
        }
        if(f==-1)for(int i=0;i<n;i++)a[i]/=n;
    }
    int main(){
        scanf("%d",&n),n--,m=2*n;
        for(int i=0;i<=n;i++)scanf("%lf%lf",&a[i].real(),&b[n-i].real());
        for(n=1;n<=m;n<<=1)L++;
        for(int i=0;i<n;i++)R[i]=(R[i>>1]>>1)|((i&1)<<(L-1));
        FFT(a,1),FFT(b,1);
        for(int i=0;i<n;i++)a[i]=a[i]*b[i];
        FFT(a,-1);
        for(int i=m/2;i<=m;i++)printf("%d
    ",(int)(a[i].real()+0.1));
    }

    这里写图片描述

  • 相关阅读:
    静态主席树
    uva 11107 Life Forms
    codeforce 605B. Lazy Student
    codeforce 606B Testing Robots
    codeforce 606C
    codeforce 606A
    uva 11019 Matrix Matcher
    uva 11468 Substring
    uvalive 4670 Dominating Patterns
    codeforce 603B
  • 原文地址:https://www.cnblogs.com/SiriusRen/p/6532063.html
Copyright © 2011-2022 走看看