【HDU5730】Shell Necklace（多项式运算，分治FFT）

最新推荐文章于 2024-11-08 12:34:34 发布
最新推荐文章于 2024-11-08 12:34:34 发布
阅读量77 收藏
点赞数
文章标签： shell
原文链接：http://www.cnblogs.com/cjyyb/p/8799206.html
版权

【HDU5730】Shell Necklace（多项式运算，分治FFT）

题面

Vjudge
翻译：
有一个长度为\(n\)的序列
已知给连续的长度为\(i\)的序列装饰的方案数为\(a[i]\)
求将\(n\)个位置全部装饰的总方案数。
答案\(mod\ 313\)

题解

很明显，是要求：
\(f[n]=\sum_{i=0}^na[i]\times f[n-i],f[0]=0\)
卷积的形式啊。。

然后就可以开始搞了

忍不住的方法一

好明显啊，把生成函数\(F,A\)给搞出来
然后就有\(F*A+1=F\)
然后\((1-A)F=1\)
然后\(F=\frac{1}{1-A}\)
多项式求逆，没了。。。
但是这个\(mod\ 313\)很蛋疼啊。。
直接蒯模板了

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define ull unsigned long long
#define RG register
#define MAX 288888
#define MOD (313)
const double Pi=acos(-1);
const int m=sqrt(MOD);
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int fpow(int a,int b){int s=1;while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}return s;}
struct Complex{double a,b;}W[MAX],A[MAX],B[MAX],C[MAX],D[MAX];
Complex operator+(Complex a,Complex b){return (Complex){a.a+b.a,a.b+b.b};}
Complex operator-(Complex a,Complex b){return (Complex){a.a-b.a,a.b-b.b};}
Complex operator*(Complex a,Complex b){return (Complex){a.a*b.a-a.b*b.b,a.a*b.b+a.b*b.a};}
int r[MAX];
void FFT(Complex *P,int N,int opt)
{
    for(int i=0;i<N;++i)if(i<r[i])swap(P[i],P[r[i]]);
    for(int i=1;i<N;i<<=1)
        for(int p=i<<1,j=0;j<N;j+=p)
            for(int k=0;k<i;++k)
            {
                Complex w=(Complex){W[N/i*k].a,W[N/i*k].b*opt};
                Complex X=P[j+k],Y=P[i+j+k]*w;
                P[j+k]=X+Y;P[i+j+k]=X-Y;
            }
    if(opt==-1)for(int i=0;i<N;++i)P[i].a/=1.0*N;
}
void Multi(int *a,int *b,int len,int *ret)
{
    for(int i=0;i<(len<<1);++i)A[i]=B[i]=C[i]=D[i]=(Complex){0,0};
    for(int i=0;i<len;++i)
    {
        a[i]%=MOD;b[i]%=MOD;
        A[i]=(Complex){(a[i]/m)*1.0,0};
        B[i]=(Complex){(a[i]%m)*1.0,0};
        C[i]=(Complex){(b[i]/m)*1.0,0};
        D[i]=(Complex){(b[i]%m)*1.0,0};
    }
    int N,l=0;
    for(N=1;N<=len;N<<=1)++l;
    for(int i=0;i<N;++i)r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
    for(int i=1;i<N;i<<=1)
        for(int k=0;k<i;++k)W[N/i*k]=(Complex){cos(k*Pi/i),sin(k*Pi/i)};
    FFT(A,N,1);FFT(B,N,1);FFT(C,N,1);FFT(D,N,1);
    for(int i=0;i<N;++i)
    {
        Complex tmp=A[i]*C[i];
        C[i]=B[i]*C[i],B[i]=B[i]*D[i],D[i]=D[i]*A[i];
        A[i]=tmp;C[i]=C[i]+D[i];
    }
    FFT(A,N,-1);FFT(B,N,-1);FFT(C,N,-1);
    for(int i=0;i<len;++i)
    {
        ret[i]=0;
        ret[i]=(ret[i]+1ll*(ll)(A[i].a+0.5)%MOD*m%MOD*m%MOD)%MOD;
        ret[i]=(ret[i]+1ll*(ll)(C[i].a+0.5)%MOD*m%MOD)%MOD;
        ret[i]=(ret[i]+1ll*(ll)(B[i].a+0.5)%MOD)%MOD;
        ret[i]=(ret[i]+MOD)%MOD;
    }
}
int c[MAX],d[MAX];
void Inv(int *a,int *b,int len)
{
    if(len==1){b[0]=fpow(a[0],MOD-2);return;}
    Inv(a,b,len>>1);
    Multi(a,b,len,c);
    Multi(c,b,len,d);
    for(int i=0;i<len;++i)b[i]=(b[i]+b[i])%MOD;
    for(int i=0;i<len;++i)b[i]=(b[i]+MOD-d[i])%MOD;
}
int n,a[MAX],b[MAX];
int main()
{
    while(n=read())
    {
        for(int i=1;i<=n;++i)a[i]=read()%MOD;
        for(int i=1;i<=n;++i)a[i]=(MOD-a[i])%MOD;
        a[0]++;
        int N;for(N=1;N<=n;N<<=1);
        Inv(a,b,N);
        printf("%d\n",b[n]);
        memset(a,0,sizeof(a));memset(b,0,sizeof(b));
    }
    return 0;
}

很生疏的方法二

这种东西显然也可以\(CDQ\)分治+\(FFT\)来做
简称分治\(FFT\)
怎么考虑？
式子长成这样：
\(f[n]=\sum_{i=0}^na[i]\times f[n-i],f[0]=0\)
一定要求出前面的才能计算后面的。
每次考虑左半段对于右半段的贡献
直接拿左半段和右半段，（要求的东西）做卷积
这样，对于每一项的系数，都是右半段每个点的一部分贡献。累加即可。
复杂度\(O(nlog^2n)\)
模数\(313\)依旧蛋疼。

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define ll long long
#define ull unsigned long long
#define RG register
#define MAX 288888
#define MOD (313)
const double Pi=acos(-1);
const int m=sqrt(MOD);
inline int read()
{
    RG int x=0,t=1;RG char ch=getchar();
    while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
    if(ch=='-')t=-1,ch=getchar();
    while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
    return x*t;
}
int fpow(int a,int b){int s=1;while(b){if(b&1)s=1ll*s*a%MOD;a=1ll*a*a%MOD;b>>=1;}return s;}
struct Complex{double a,b;}W[MAX],A[MAX],B[MAX],C[MAX],D[MAX];
Complex operator+(Complex a,Complex b){return (Complex){a.a+b.a,a.b+b.b};}
Complex operator-(Complex a,Complex b){return (Complex){a.a-b.a,a.b-b.b};}
Complex operator*(Complex a,Complex b){return (Complex){a.a*b.a-a.b*b.b,a.a*b.b+a.b*b.a};}
int r[MAX];
void FFT(Complex *P,int N,int opt)
{
    for(int i=0;i<N;++i)if(i<r[i])swap(P[i],P[r[i]]);
    for(int i=1;i<N;i<<=1)
        for(int p=i<<1,j=0;j<N;j+=p)
            for(int k=0;k<i;++k)
            {
                Complex w=(Complex){W[N/i*k].a,W[N/i*k].b*opt};
                Complex X=P[j+k],Y=P[i+j+k]*w;
                P[j+k]=X+Y;P[i+j+k]=X-Y;
            }
    if(opt==-1)for(int i=0;i<N;++i)P[i].a/=1.0*N;
}
void Multi(int *a,int *b,int len,int *ret)
{
    for(int i=0;i<(len<<1);++i)A[i]=B[i]=C[i]=D[i]=(Complex){0,0};
    for(int i=0;i<len;++i)
    {
        a[i]%=MOD;b[i]%=MOD;
        A[i]=(Complex){(a[i]/m)*1.0,0};
        B[i]=(Complex){(a[i]%m)*1.0,0};
        C[i]=(Complex){(b[i]/m)*1.0,0};
        D[i]=(Complex){(b[i]%m)*1.0,0};
    }
    int N,l=0;
    for(N=1;N<=len;N<<=1)++l;
    for(int i=0;i<N;++i)r[i]=(r[i>>1]>>1)|((i&1)<<(l-1));
    for(int i=1;i<N;i<<=1)
        for(int k=0;k<i;++k)W[N/i*k]=(Complex){cos(k*Pi/i),sin(k*Pi/i)};
    FFT(A,N,1);FFT(B,N,1);FFT(C,N,1);FFT(D,N,1);
    for(int i=0;i<N;++i)
    {
        Complex tmp=A[i]*C[i];
        C[i]=B[i]*C[i],B[i]=B[i]*D[i],D[i]=D[i]*A[i];
        A[i]=tmp;C[i]=C[i]+D[i];
    }
    FFT(A,N,-1);FFT(B,N,-1);FFT(C,N,-1);
    for(int i=0;i<len;++i)
    {
        ret[i]=0;
        ret[i]=(ret[i]+1ll*(ll)(A[i].a+0.5)%MOD*m%MOD*m%MOD)%MOD;
        ret[i]=(ret[i]+1ll*(ll)(C[i].a+0.5)%MOD*m%MOD)%MOD;
        ret[i]=(ret[i]+1ll*(ll)(B[i].a+0.5)%MOD)%MOD;
        ret[i]=(ret[i]+MOD)%MOD;
    }
}
int n,a[MAX],x[MAX],y[MAX],f[MAX];
void CDQ(int l,int r)
{
    if(l==r){f[l]=(f[l]+a[l])%MOD;return;}
    int mid=(l+r)>>1,N;
    CDQ(l,mid);for(N=1;N<=(r-l+1);N<<=1);
    for(int i=0;i<=N;++i)x[i]=y[i]=0;
    for(int i=l;i<=mid;++i)x[i-l]=f[i];
    for(int i=0;i<=r-l;++i)y[i]=a[i];
    Multi(x,y,N,x);
    for(int i=mid+1;i<=r;++i)f[i]=(f[i]+x[i-l])%MOD;
    CDQ(mid+1,r);
}
int main()
{
    while(n=read())
    {
        for(int i=1;i<=n;++i)a[i]=read()%MOD;
        int N;for(N=1;N<=n;N<<=1);
        CDQ(1,n);
        printf("%d\n",f[n]);
        memset(f,0,sizeof(f));
        memset(a,0,sizeof(a));
    }
    return 0;
}

转载于:https://www.cnblogs.com/cjyyb/p/8799206.html

hdu5730(FFT+多项式求逆）
KKiseki的博客
09-05 427
（最近在班上好无聊，总感觉同学都好拽，也没有能说很多话的小姐姐，而且都很爱学习。老师也很死板，诸多要求。但我还是我，毕竟看过那么多治愈番。） 终有一天会降临，只属于我的苍蓝天使！ 稳定的题面传送门 题意是这样的，长度分别为1..n的棒子各有a[i]种，每种都有无限条，问组成长度为n的大棒子有多少种方案，mod313。 设f[i]表示长度为i的方案数，根据题意f[0]=1，考虑dp。 枚举
（HDU2011 C++）多项式求和
qq_38780240的博客
08-06 858
多项式求和 Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K (Java/Others) Total Submission(s): 96703 Accepted Submission(s): 54997 Problem Description 多项式的描述如下： 1 - 1/2 ...
shell中实现多种计算方式
qq_46020608的博客
02-22 165
第一种： expr 2 + 7 #注意空格不能少 否则报错；只能计算整数 第二种： echo "2 + 7" |bc #bc用来计算，可以计算小数 echo "scale=3; 8/3" |bc #scale指定小数点位数为3 第三种： echo $(( 2 + 7 )) #一定要空格 第四种： echo $[ 2 + 7 ] #同3 ...
多项式乘积的分治算法
05-14
多项式乘积的分治算法
bzoj 3456 城市规划——分治FFT / 多项式求逆 / 多项式求ln
weixin_34122604的博客
11-30 121
题目：https://www.lydsy.com/JudgeOnline/problem.php?id=3456 分治FFT： 　　设 dp[ i ] 表示 i 个点时连通的方案数。 　　考虑算补集：连通的方案数 == 随便连方案数 - 不连通方案数 　　不连通方案数就和很久之前做过的“地震后的幻想乡”一样，枚举一个连通的点集，其中需要一直包含一个“划分点”保证不重复；其余部分随便连。注意...
HDU 5730 Shell Necklace（CDQ分治+FFT）
ZSQ
07-25 2387
Description 给出长度分别为1~n的珠子，长度为i的珠子有a[i]种，每种珠子有无限个，问用这些珠子串成长度为n的链有多少种方案 Input 多组用例，每组用例首先输入一整数n表示链长，之后n个整数ai表示长度为i的珠子种类数，以n=0结束输入(n<=10^5,0<=ai<=10^7) Output 对每组用例，输出方案数，结果模313 Sample Input 3 1 3
多项式计算-FFT与NTT学习小结
skyword的专栏
01-07 811
【写于2016-10-11，现在从hexo搬运过来】 看了两天FFT 一直没有学fft，直到10天前打beijing online contest时，碰到了一个NTT的裸题，发现必须要看看了 hexo博客配置有一些问题，mathjax公式一直不能正常加载，可以访问我在cmdmarkdown上发布的版本 here 这几天看了一些博客，发现它纯粹在数学上需要功底。FFT,NTT的题目中，简单的那部分...
hdu1007最近点对-分治策略-具体优化
他的博客
03-08 902
1)AC总是让人感到兴奋。 上一篇求最大连续子数组之和了解到分治策略之后，本题先是找到中间界线，依然分成三部分，左边集合、右边集合、两个点分别在左边和右边的集合的情况。前两种情况不断递归，直到该集合内只有两个点，则求两点之间距离并返回，第三种情况如果常规来做，则是左边集合中的每一个点都要对应求到右边集合中每个点的距离，比较出其中最小的，此处O（(n/2)^2），时间复杂度依然是n^2层次，于是继
hdu 1007 平面最近点对 分治
JIBANCANYANG
07-26 350
平面上的分治，把取余分割，然后合并的时候有剪枝。 这里学会了归并合并函数inplace_merge(a, a + m, a + n, cmp) 以mm为分界线分别已经排序好了，只需要按照小于操作符合并就可以了。当然自己实现也只是双指针的问题。#include <cstdio> #include <cstring> #include <cmath> #include <iostream> #inc
分治算法--多项式
gy的博客
06-04 747
from cmath import pi, exp def FFT(A, w): """ 快速傅氏变换 要求将多项式的系数表示长度为2的次方 :param A:多项式 :param w:单位的n次方本原单位根 :return: """ length = len(A) if length == 1: retur...
shell脚本中的数值运算
热门推荐
seteor的专栏
12-17 3万+
shell中的变量没有如c语言中字符型、整型、浮点型的区分，默认将变量赋值当做字符，怎样才能让shell将变量赋值当做数值型呢，答案是用$ 和[], []包含的数学算式可以被shell正确识别。 [root@localhost ~]# cat  test.sh #!/bin/bash var1=100 var2=200 var3=$[$var1+$var2] echo "v
【转】Linux shell 变量 数学 运算
weixin_34248705的博客
03-30 354
源地址：http://blog.chinaunix.net/uid-20671208-id-3552751.html Abstract : 1) Linux shell 中使用 let , [ ] ,(( )) 三种运算符操作 shell 变量进行简单的基本运算；2）Linux shell 中使用 expr 与 bc 两个程序实现高级运算； 1, Linux shell 变量的基本运...
shell中的运算（包括使用多种方法实现加法运算）
ly_qiu的博客
04-27 5965
1.运算符号 符号 作用 + - * / % ** ++ - - < <= > >= = += *= /= %= 2.运算指令 运算指令 作用 (()) ((a=1+2)) let let a=1+2 expr...
linux之shell:数值运算：简单加减乘除
yusiyu的技术专栏
01-06 4835
sh04.sh #!/bin/bash #program #  user inputs 2 integer numbers;program will cross these two numbers. #History: #2014/01/06 yusiyu First release PATH=/bin:/sbin:/usr/bin:/usr/sbin:/usr/local/b
[vulnhub] DarkHole: 1
最新发布
weixin_46099552的博客
11-08 696
vulnhub - DarkHole: 1
day16：shell脚本的表达式与流程控制
qq_62174860的博客
11-05 867
条件判断语句用于检测某个条件是否满足，从而决定是否执行特定的操作。Shell 脚本中常见的条件判断语句主要包括 、 等。if 语句： case 语句： 常见判断的参数 文件判断参数 参数 说明 案例代码 解释 比如 文件的判断 文件/目录存在 判断 是否存在 普通文件 判断 是否为普通文件 是否为目录 判断 是否为目录 文件的权限判断 文件/目录
Fastjson反序列化漏洞利用教程
weixin_38055730的博客
11-05 830
Fastjson反序列化漏洞是一种高危安全漏洞，攻击者通过构造恶意的JSON字符串可以执行任意命令。此漏洞影响广泛，需要引起高度重视。通过本文的学习和实践，读者可以掌握Fastjson反序列化漏洞的基本原理、复现方法和利用技巧。
Shell 脚本入门——变量和字符串
licy__的博客
11-07 264
在 Shell 脚本中，变量的定义非常简单。你只需要将值赋给变量名，不需要声明变量类型。注意：赋值时等号两边不能有空格。
HDU 2000-2099题目源码完整分享
资源摘要信息:"HDU-2000-2099" 知识点概览： HDU 2000-2099是一系列编程题目集，这些题目通常来源于中国一个著名的在线编程评测网站——杭电在线评测系统（HANGZHOU DAILY UNIVERSITY Online Judge，简称HDU OJ）。...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值