[离散对数 FFT FastIO] ICPC 2016 Shenyang K & HDU 5958 New Signal Decomposition

最新推荐文章于 2019-04-07 00:18:53 发布
最新推荐文章于 2019-04-07 00:18:53 发布
阅读量867 收藏 2
点赞数
分类专栏： FFT&FWT 数论 高精度&FastIO模板
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/u014609452/article/details/68057978
版权

特判0 然后将所有下标取对数 那么就是
bk′=∑ah′∗ch′+k′
使用FFT解决

本来好好的一道题 卡常数+卡空间+卡精度？！

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cmath>
using namespace std;
typedef double ld;

namespace fastIO{
    #define BUF_SIZE 100000
    #define OUT_SIZE 100000
    #define ll long long
    //fread->read
    bool IOerror=0;
    inline char nc(){
        static char buf[BUF_SIZE],*p1=buf+BUF_SIZE,*pend=buf+BUF_SIZE;
        if (p1==pend){
            p1=buf; pend=buf+fread(buf,1,BUF_SIZE,stdin);
            if (pend==p1){IOerror=1;return -1;}
            //{printf("IO error!\n");system("pause");for (;;);exit(0);}
        }
        return *p1++;
    }
    inline bool blank(char ch){return ch==' '||ch=='\n'||ch=='\r'||ch=='\t';}
    inline void read(int &x){
        bool sign=0; char ch=nc(); x=0;
        for (;blank(ch);ch=nc());
        if (IOerror)return;
        if (ch=='-')sign=1,ch=nc();
        for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';
        if (sign)x=-x;
    }
    inline void read(ll &x){
        bool sign=0; char ch=nc(); x=0;
        for (;blank(ch);ch=nc());
        if (IOerror)return;
        if (ch=='-')sign=1,ch=nc();
        for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';
        if (sign)x=-x;
    }
    inline void read(double &x){
        bool sign=0; char ch=nc(); x=0;
        for (;blank(ch);ch=nc());
        if (IOerror)return;
        if (ch=='-')sign=1,ch=nc();
        for (;ch>='0'&&ch<='9';ch=nc())x=x*10+ch-'0';
        if (ch=='.'){
            double tmp=1; ch=nc();
            for (;ch>='0'&&ch<='9';ch=nc())tmp/=10.0,x+=tmp*(ch-'0');
        }
        if (sign)x=-x;
    }
    inline void read(char *s){
        char ch=nc();
        for (;blank(ch);ch=nc());
        if (IOerror)return;
        for (;!blank(ch)&&!IOerror;ch=nc())*s++=ch;
        *s=0;
    }
    inline void read(char &c){
        for (c=nc();blank(c);c=nc());
        if (IOerror){c=-1;return;}
    }
    //getchar->read
    inline void read1(int &x){
        char ch;int bo=0;x=0;
        for (ch=getchar();ch<'0'||ch>'9';ch=getchar())if (ch=='-')bo=1;
        for (;ch>='0'&&ch<='9';x=x*10+ch-'0',ch=getchar());
        if (bo)x=-x;
    }
    inline void read1(ll &x){
        char ch;int bo=0;x=0;
        for (ch=getchar();ch<'0'||ch>'9';ch=getchar())if (ch=='-')bo=1;
        for (;ch>='0'&&ch<='9';x=x*10+ch-'0',ch=getchar());
        if (bo)x=-x;
    }
    inline void read1(double &x){
        char ch;int bo=0;x=0;
        for (ch=getchar();ch<'0'||ch>'9';ch=getchar())if (ch=='-')bo=1;
        for (;ch>='0'&&ch<='9';x=x*10+ch-'0',ch=getchar());
        if (ch=='.'){
            double tmp=1;
            for (ch=getchar();ch>='0'&&ch<='9';tmp/=10.0,x+=tmp*(ch-'0'),ch=getchar());
        }
        if (bo)x=-x;
    }
    inline void read1(char *s){
        char ch=getchar();
        for (;blank(ch);ch=getchar());
        for (;!blank(ch);ch=getchar())*s++=ch;
        *s=0;
    }
    inline void read1(char &c){for (c=getchar();blank(c);c=getchar());}
    //scanf->read
    inline void read2(int &x){scanf("%d",&x);}
    inline void read2(ll &x){
        #ifdef _WIN32
            scanf("%I64d",&x);
        #else
        #ifdef __linux
            scanf("%lld",&x);
        #else
            puts("error:can't recognize the system!");
        #endif
        #endif
    }
    inline void read2(double &x){scanf("%lf",&x);}
    inline void read2(char *s){scanf("%s",s);}
    inline void read2(char &c){scanf(" %c",&c);}
    inline void readln2(char *s){gets(s);}
    //fwrite->write
    struct Ostream_fwrite{
        char *buf,*p1,*pend;
        Ostream_fwrite(){buf=new char[BUF_SIZE];p1=buf;pend=buf+BUF_SIZE;}
        void out(char ch){
            if (p1==pend){
                fwrite(buf,1,BUF_SIZE,stdout);p1=buf;
            }
            *p1++=ch;
        }
        void print(int x){
            static char s[15],*s1;s1=s;
            if (!x)*s1++='0';if (x<0)out('-'),x=-x;
            while(x)*s1++=x%10+'0',x/=10;
            while(s1--!=s)out(*s1);
        }
        void println(int x){
            static char s[15],*s1;s1=s;
            if (!x)*s1++='0';if (x<0)out('-'),x=-x;
            while(x)*s1++=x%10+'0',x/=10;
            while(s1--!=s)out(*s1); out('\n');
        }
        void print(ll x){
            static char s[25],*s1;s1=s;
            if (!x)*s1++='0';if (x<0)out('-'),x=-x;
            while(x)*s1++=x%10+'0',x/=10;
            while(s1--!=s)out(*s1);
        }
        void println(ll x){
            static char s[25],*s1;s1=s;
            if (!x)*s1++='0';if (x<0)out('-'),x=-x;
            while(x)*s1++=x%10+'0',x/=10;
            while(s1--!=s)out(*s1); out('\n');
        }
        void print(double x,int y){
            static ll mul[]={1,10,100,1000,10000,100000,1000000,10000000,100000000,
                1000000000,10000000000LL,100000000000LL,1000000000000LL,10000000000000LL,
                100000000000000LL,1000000000000000LL,10000000000000000LL,100000000000000000LL};
            if (x<-1e-12)out('-'),x=-x;x*=mul[y];
            ll x1=(ll)floor(x); if (x-floor(x)>=0.5)++x1;
            ll x2=x1/mul[y],x3=x1-x2*mul[y]; print(x2);
            if (y>0){out('.'); for (size_t i=1;i<y&&x3*mul[i]<mul[y];out('0'),++i); print(x3);}
        }
        void println(double x,int y){print(x,y);out('\n');}
        void print(char *s){while (*s)out(*s++);}
        void println(char *s){while (*s)out(*s++);out('\n');}
        void flush(){if (p1!=buf){fwrite(buf,1,p1-buf,stdout);p1=buf;}}
        ~Ostream_fwrite(){flush();}
    }Ostream;
    inline void print(int x){Ostream.print(x);}
    inline void println(int x){Ostream.println(x);}
    inline void print(char x){Ostream.out(x);}
    inline void println(char x){Ostream.out(x);Ostream.out('\n');}
    inline void print(ll x){Ostream.print(x);}
    inline void println(ll x){Ostream.println(x);}
    inline void print(double x,int y){Ostream.print(x,y);}
    inline void println(double x,int y){Ostream.println(x,y);}
    inline void print(char *s){Ostream.print(s);}
    inline void println(char *s){Ostream.println(s);}
    inline void println(){Ostream.out('\n');}
    inline void flush(){Ostream.flush();}
    //puts->write
    char Out[OUT_SIZE],*o=Out;
    inline void print1(int x){
        static char buf[15];
        char *p1=buf;if (!x)*p1++='0';if (x<0)*o++='-',x=-x;
        while(x)*p1++=x%10+'0',x/=10;
        while(p1--!=buf)*o++=*p1;
    }
    inline void println1(int x){print1(x);*o++='\n';}
    inline void print1(ll x){
        static char buf[25];
        char *p1=buf;if (!x)*p1++='0';if (x<0)*o++='-',x=-x;
        while(x)*p1++=x%10+'0',x/=10;
        while(p1--!=buf)*o++=*p1;
    }
    inline void println1(ll x){print1(x);*o++='\n';}
    inline void print1(char c){*o++=c;}
    inline void println1(char c){*o++=c;*o++='\n';}
    inline void print1(char *s){while (*s)*o++=*s++;}
    inline void println1(char *s){print1(s);*o++='\n';}
    inline void println1(){*o++='\n';}
    inline void flush1(){if (o!=Out){if (*(o-1)=='\n')*--o=0;puts(Out);}}
    struct puts_write{
        ~puts_write(){flush1();}
    }_puts;
    inline void print2(int x){printf("%d",x);}
    inline void println2(int x){printf("%d\n",x);}
    inline void print2(char x){printf("%c",x);}
    inline void println2(char x){printf("%c\n",x);}
    inline void print2(ll x){
        #ifdef _WIN32
            printf("%I64d",x);
        #else
        #ifdef __linux
            printf("%lld",x);
        #else
            puts("error:can't recognize the system!");
        #endif
        #endif
    }
    inline void println2(ll x){print2(x);printf("\n");}
    inline void println2(){printf("\n");}
    #undef ll
    #undef OUT_SIZE
    #undef BUF_SIZE
};
using namespace fastIO;

const int N=300005;

int P,G;
int ind[N],pw[N];

inline void Init(){
  G=P==103?5:2;
  int t=1;
  for (int i=0;i<P;i++)
    pw[i]=t,ind[t]=i,t=t*G%P;
}

double a[N];

struct E{
  ld real,imag;  
  E(ld real=0,ld imag=0):real(real) , imag(imag) { }
  void operator /=(ld a){ real/=a; imag/=a; }
  friend E operator+(E&a,E&b){ return E(a.real+b.real,a.imag+b.imag); }  
  friend E operator-(E&a,E&b){ return E(a.real-b.real,a.imag-b.imag); }  
  friend E operator*(E&a,E&b){ return E(a.real*b.real-a.imag*b.imag,a.imag*b.real+a.real*b.imag); }  
};
const ld PI=acos(-1.0);

int R[N];

inline void FFT(E *a,int n,int r){  
  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){  
    E wn(cos(PI/i),r*sin(PI/i));  
    for (int j=0;j<n;j+=(i<<1)){  
      E w(1,0);  
      for (int k=0;k<i;k++,w=w*wn){  
    E x=a[j+k],y=w*a[j+k+i];  
    a[j+k]=x+y; a[j+k+i]=x-y;  
      }  
    }  
  }  
  if (r==-1) for (int i=0;i<n;i++) a[i]/=n;  
}

int n1,n2,m,L;
double a1[N],a2[N];
E A[N],B[N];

double ans[N];

inline double C(int x){
  return pow(2,pow(sin(2*PI*x/P),3));
}

int main(){
  while (1){
    P=0; read2(P); if (!P) break;
    Init();
    ans[0]=0;
    for (int i=0;i<P;i++) read2(a[i]),ans[0]+=a[i];
    for (m=1,L=-1;m<=2*P;m<<=1,L++); 
    for (int i=0;i<P-1;i++) A[i]=E(a[pw[P-2-i]],0),B[i]=E(C(pw[i]),0);
    for (int i=P-1;i<m;i++) A[i]=B[i]=E(0,0);
    for (int i=1;i<m;i++) R[i]=(R[i>>1]>>1)|((i&1)<<L);
    FFT(A,m,1);
    FFT(B,m,1);
    for (int i=0;i<m;i++) A[i]=A[i]*B[i];
    FFT(A,m,-1);
    for (int i=P-1;i<m;i++) A[i-(P-1)].real+=A[i].real;
    for (int i=1;i<P;i++) ans[pw[i]]=a[0]+A[(i-1)%(P-1)].real;
    for (int i=0;i<P;i++)
      print(ans[i],3),print(' ');
    print('\n');
  }
  return 0;
}
里阿奴摩西
关注
  • 0
    点赞
  • 踩
  • 2
    收藏
    觉得还不错? 一键收藏
  • 2
    评论
专栏目录
fft.rar_功率曲线_功率谱_对数功率谱_离散型FFT_离散数据fft
09-23
标题提到的"fft.rar"通常是一个包含相关代码或示例的压缩包，用于演示如何通过离散傅立叶变换计算功率谱，并绘制双对数曲线。 傅立叶变换是一种数学工具，能够将一个在时间（或空间）域内的信号转换为频率域的表示...
ACM-数论完全总结（知识点+模板）
热门推荐
Ogmx的博客
08-30 3万+
一.整除的性质： 1.若a|b &amp;amp;amp;amp;amp;amp;amp;lt;-&amp;amp;amp;amp;amp;amp;amp;gt; -a|b &amp;amp;amp;amp;amp;amp;amp;lt;-&amp;amp;amp;amp;amp;amp;amp;gt; a|-b &amp;amp;amp;amp;amp;amp;amp;lt;-&amp;amp;amp;amp;amp;amp;amp;gt; |a| |
Bzoj 3122 离散对数
ACM, deep love
09-19 254
题目描述 HINT 这题的p是素数。。。尼玛想半天不会。。。会合数版本的请联系我 #include &lt;bits/stdc++.h&gt; using namespace std; typedef long long ll; map&lt;ll,int&gt;mp; void exgcd(ll a, ll b, ll &amp;x, ll &amp;y){ if(b==0){ x=1,...
离散对数求解算法
小周的专栏
01-21 1万+
///大步小步算法 struct baby///小步算法预存的结构体定义 { long long b,j; bool operator < (const baby &other)const{ return b<other.b; } }babyv[100005]; ///快速幂，x^y%mod long long q_pow(long long x,long
离散对数
01-15 2834
写在前面： 　　学习笔记，方便复习 　　非原创标明出处 我们都在努力奔跑，我们都是追梦人 byskecchiart 概念 对数 　　在数学中，对数是对求幂的逆运算，正如除法是乘法的倒数，反之亦然。 这意味着一个数字的对数是必须产生另一个固定数字（基数）的指数 　　如果a的x次方等于N（a>0，且a不等于1），那么数x叫做以a为底N的...
原根和离散对数BSGS求法（高次同余方程）
zjyang12345的博客
04-07 842
原根&离散对数 一.原根 1.定义: （a与m互质）使成立的最小的d(记住原根是a，不是d!) 2.原根的性质:一般给出p（有时叫m） 1.具有原根的数字仅有以下几种形式:,(p是奇质数) 2.一个数的最小原根的大小不超过 3.原根个数Φ(Φ(m))个，m为质数则原根个数Φ(m-1） 3.求解原根的基本步骤: 判断一个数是否有原根。(通过性质1,枚举质数即可) 求得最小...
test_离散信号fft_fft_matlab_离散FFT_earthcy8_
09-30
本文将详细介绍离散信号的FFT（快速傅立叶变换）以及如何在MATLAB中实现这一过程，并以"test_离散信号fft_fft_matlab_离散FFT_earthcy8_"为实例，解析相关知识点。 一、离散傅立叶变换（DFT） 离散傅立叶变换是一种...
sampe-fft_离散信号fft_离散_fft_
10-03
标题“sampe-fft_离散信号fft_离散_fft_”表明我们要讨论的是关于离散信号的FFT应用示例。MATLAB的`.mdl`文件扩展名通常用于Simulink模型文件，这暗示我们将在Simulink环境中探讨如何实现离散信号的FFT。 离散信号...
test1.rar_离散FFT_离散信号fft_离散数据fft
07-15
将离散的数据或信号，需要进行简单的傅里叶变换进行频谱分析。采用调用fft函数直接进行傅里叶变换。
fft.rar_离散信号 FFT
最新发布
09-14
离散傅立叶变换（Discrete Fourier Transform, DFT）是数字信号处理中的核心概念，它在音频处理、图像分析、通信系统等多个领域都有广泛应用。快速傅立叶变换（Fast Fourier Transform, FFT）是对DFT的一种高效计算...
2018ICPC沈阳网络赛题解
09-09
这是2018ICPC沈阳网络赛的题解，这是2018ICPC沈阳网络赛的题解
SPOJ MOD 离散对数
ACM, deep love
09-19 438
题目描述 求方程A^X=B%C的最小非负整数解，这里的C不一定是素数。 #include &lt;bits/stdc++.h&gt; using namespace std; typedef long long ll; //A^X=B%C unordered_map&lt;ll,int&gt;mp; ll exbsgs(ll a,ll b,ll c){ if(c==1)return 0; ll...
离散对数（关于方程x^A=B(mod C)的解)
ACdreamer
08-05 6755
求方程：的解个数 分析：设，那么上述方程解的个数就与同余方程组：的解等价。 设同于方程的解分别是：，那么原方程的解的个数就是 所以现在的关键问题是求方程：的解个数。 这个方程我们需要分3类讨论： 第一种情况： 对于这种情况，如果方程的某个解设为，那么一定有，可以得到，即 所以方程的解个数就是：，也就是 第二种情
[转载]FFT实践及频谱分析
da895的专栏
02-21 916
*************************************************************************% %                              FFT实践及频谱分析                          % %*****************************************************
离散对数:这个好难
duaduac.com
08-17 4153
定义 设g是m的一个原根，对于满足(k, m) = 1 的k, k关于g的离散对数(mod m)定义为一整数t,使gt≡k(modm)g^t \equiv k \pmod m 且t为一个最小剩余(modϕ(m))\pmod {\phi(m)}。 记作indgkind_gk。          //此处可类比整式里的loggklog_gk, 不过这里对ϕ(m)\phi(m)取模了。 举
【HDU】5958 New Signal Decomposition【离散对数下的FFT】
poursoul
11-01 1246
题目链接：【HDU】5958 New Signal Decomposition在此先感谢小q对我的指导，没有q老师的帮助，估计永远也做不出来了。首先我们考虑对这个式子做离散对数。令gg为pp的某个原根，则有： bi=∑p−1j=0aj⋅r(i,j)\quad b_i=\sum_{j=0}^{p-1}a_j\cdot r(i,j) bi=∑p−1j=0aj⋅2sin32πi⋅jp\quad b_i
HDU - 5958 New Signal Decomposition 原根下标变换+fft(未完待续)
我们做的队
07-14 463
题目链接点这里 过程看这里http://blog.csdn.net/u013368721/article/details/53001532 ，，然而，，这位大佬的程序也gg了，，可能改过数据了，，不过思路是这样的，， 然而我与标程对拍，，答案总是差0.00几。。反正，，我写的fft自带低精度，大内存，高常数。。，，以后再看吧 #include #include #include #incl
B&K PULSE系统全面指南与校准方法
B&K PULSE 是丹麦B&K公司在1996年推出的一款创新的噪声与振动多分析仪系统，这款系统的设计初衷是为了满足用户的多元化需求，包括多通道、实时分析、频谱分析（FFT）、连续波谱分析(CPB)以及总级值计算等。...
评论 2
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值