博客
关于我
2016vijos 1-3 兔子的晚会(生成函数+倍增FWT)
阅读量:798 次
发布时间:2023-04-16

本文共 1683 字,大约阅读时间需要 5 分钟。

求出序列的生成函数后,倍增FWT

在解决某些序列问题时,我们需要求出序列的生成函数,并对其进行倍增快速傅里叶变换(FWT)。以下是实现这一过程的详细方法。

首先,我们在代码中定义了常量和模数:

#define N 2048const int mod=1e9+7;int inv;int f[N+1];

接下来,定义了一个快速幂的函数来计算大数的幂取模:

int Pow(int a,int b){    int res=1;    for(;b;a=1LL*a*a%mod,b >>=1)        if(b & 1) res = 1LL*res*a % mod;    return res;}

然后,实现了FWT变换和逆变换的函数。这些函数用于对序列进行快速傅里叶变换和逆变换:

void FWT(int *a,int n){    int x,y;    for(int d=1;d < n;d <<=1){        for(int m=d; m < n; m <<=1){            for(int i=0; i < n; i += m){                for(int j=0; j < d; j++){                    x = a[i+j];                    y = a[i+j+d];                    a[i+j] = (x + y) % mod;                    a[i+j+d] = (x - y + mod) % mod;                    a[i+j] -= a[i+j];                    a[i+j] %= mod;                    a[i+j+d] += a[i+j+d] < 0 ? mod : 0;                }            }        }    }}void IFWT(int *a,int n){    int x,y;    for(int d=1;d < n;d <<=1){        for(int m=d; m < n; m <<=1){            for(int i=0; i < n; i += m){                for(int j=0; j < d; j++){                    x = a[i+j];                    y = a[i+j+d];                    a[i+j] = (1LL*(x + y) * inv % mod);                    a[i+j+d] = (1LL*(x - y + mod) % mod * inv % mod);                }            }        }    }}

最后,在主函数中,我们读取输入并初始化变量:

int main(){    freopen("xor.in","r",stdin);    freopen("xor.out","w",stdout);    int n,m,l,r;    scanf("%d%d%d%d",&n,&m,&l,&r);    n = 2*n +1;    inv = Pow(2,mod-2);    len = 1;    ans = 0;    while(len <= x + m) len <<=1;    for(i=0; i < len; i++)        f[i] = Pow(f[i],n);    fwt(f,len);    ifwt(f,len);    ans += f[0];    printf("%d",ans);}

通过以上步骤,我们可以高效地求出序列的生成函数并对其进行倍增FWT变换,解决实际问题。

转载地址:http://fjgfk.baihongyu.com/

你可能感兴趣的文章
mysql中出现Incorrect DECIMAL value: '0' for column '' at row -1错误解决方案
查看>>
mysql中出现Unit mysql.service could not be found 的解决方法
查看>>
mysql中出现update-alternatives: 错误: 候选项路径 /etc/mysql/mysql.cnf 不存在 dpkg: 处理软件包 mysql-server-8.0的解决方法(全)
查看>>
Mysql中各类锁的机制图文详细解析(全)
查看>>
MySQL中地理位置数据扩展geometry的使用心得
查看>>
Mysql中存储引擎简介、修改、查询、选择
查看>>
Mysql中存储过程、存储函数、自定义函数、变量、流程控制语句、光标/游标、定义条件和处理程序的使用示例
查看>>
mysql中实现rownum,对结果进行排序
查看>>
mysql中对于数据库的基本操作
查看>>
Mysql中常用函数的使用示例
查看>>
MySql中怎样使用case-when实现判断查询结果返回
查看>>
Mysql中怎样使用update更新某列的数据减去指定值
查看>>
Mysql中怎样设置指定ip远程访问连接
查看>>
mysql中数据表的基本操作很难嘛,由这个实验来带你从头走一遍
查看>>
Mysql中文乱码问题完美解决方案
查看>>
mysql中的 +号 和 CONCAT(str1,str2,...)
查看>>
Mysql中的 IFNULL 函数的详解
查看>>
mysql中的collate关键字是什么意思?
查看>>
MySql中的concat()相关函数
查看>>
mysql中的concat函数,concat_ws函数,concat_group函数之间的区别
查看>>