//和以前写的fft不太一样,可能是因为要取模??
#include<iostream>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;
const int mod=,maxn=;
int mx,n,m,inv[maxn],a[maxn],b[maxn],c[maxn],na[maxn],w[][maxn],pos[maxn];
int qmi(int x,int y){
int t=;
for(;y;y>>=,x=(ll)x*x%mod)if(y&)t=(ll)t*x%mod;
return t;
}
void pre(int n){
int i,x=qmi(,(mod-)/n);//以前这里的取值都和mod无关啊,取了模了不一样了?
w[][]=w[][]=;
for(int i=;i<n;++i)w[][i]=w[][n-i]=(ll)w[][i-]*x%mod;
for(int i=;i<n;++i){
pos[i]=pos[i>>]>>;
if(i&)pos[i]|=n>>;
}
}
void fnt(int *a,int n,int flag){
if(n>mx)mx=n;
int i,j,k,l,x,u,v;
for(i=;i<n;++i)na[pos[i]]=a[i];
memcpy(a,na,sizeof(int)*n);
for(k=;k<n;k<<=){
for(i=,x=n/k>>;i<n;i+=k<<)
for(j=i,l=;j<i+k;++j,l+=x){
u=a[j];v=(ll)a[j+k]*w[flag][l]%mod;
a[j]=(u+v)%mod;a[j+k]=(u-v+mod)%mod;
}
}
if(flag){
x=qmi(n,mod-);
for(i=;i<n;++i)a[i]=(ll)a[i]*x%mod;
}
}
void solve_inv(int *a,int *b,int n){
if(n==){b[]=qmi(a[],mod-);return;}
int i;solve_inv(a,b,n>>);
memcpy(c,a,sizeof(int)*n);memset(c+n,,sizeof(int)*n);
pre(n<<);
fnt(b,n<<,);fnt(c,n<<,);
for(i=;i<(n<<);++i)b[i]=(-(ll)b[i]*c[i]%mod+mod)*b[i]%mod;
fnt(b,n<<,);memset(b+n,,sizeof(int)*n);
}
int main(){
int i,n;scanf("%d",&n);
inv[]=inv[]=a[]=m=;
while(m<=n)m<<=;
for(int i=;i<=n;++i)inv[i]=mod-(ll)inv[mod%i]*(mod/i)%mod;
for(int i=;i<=n;++i)inv[i]=(ll)inv[i-]*inv[i]%mod;
for(int i=;i<=n;++i)a[i]=((mod-inv[i])<<)%mod;
solve_inv(a,b,m);
int ans=b[n];
for(int i=n;i;--i)ans=((ll)ans*i+b[i-])%mod;
printf("%d\n",ans);
return ;
}

学习地址:http://blog.csdn.net/lych_cys/article/details/51512278

bzoj4555(多项式求逆解法)的更多相关文章

  1. 【BZOJ】4555: [Tjoi2016&Heoi2016]求和 排列组合+多项式求逆 或 斯特林数+NTT

    [题意]给定n,求Σi=0~nΣj=1~i s(i,j)*2^j*j!,n<=10^5. [算法]生成函数+排列组合+多项式求逆 [题解]参考: [BZOJ4555][Tjoi2016& ...

  2. luogu P4725 多项式对数函数 (模板题、FFT、多项式求逆、求导和积分)

    手动博客搬家: 本文发表于20181125 13:25:03, 原地址https://blog.csdn.net/suncongbo/article/details/84487306 题目链接: ht ...

  3. hdu 5730 Shell Necklace [分治fft | 多项式求逆]

    hdu 5730 Shell Necklace 题意:求递推式\(f_n = \sum_{i=1}^n a_i f_{n-i}\),模313 多么优秀的模板题 可以用分治fft,也可以多项式求逆 分治 ...

  4. BZOJ 4555: [Tjoi2016&Heoi2016]求和 [分治FFT 组合计数 | 多项式求逆]

    4555: [Tjoi2016&Heoi2016]求和 题意:求\[ \sum_{i=0}^n \sum_{j=0}^i S(i,j)\cdot 2^j\cdot j! \\ S是第二类斯特林 ...

  5. NTT+多项式求逆+多项式开方(BZOJ3625)

    定义多项式$h(x)$的每一项系数$h_i$,为i在c[1]~c[n]中的出现次数. 定义多项式$f(x)$的每一项系数$f_i$,为权值为i的方案数. 通过简单的分析我们可以发现:$f(x)=\fr ...

  6. Re.多项式求逆

    前言 emmm暂无 多项式求逆目的 顾名思义 就是求出一个多项式的摸xn时的逆 给定一个多项式F(x),请求出一个多项式G(x),满足F(x)∗G(x)≡1(modxn),系数对998244353取模 ...

  7. BZOJ 3456: 城市规划 与 多项式求逆算法介绍(多项式求逆, dp)

    题面 求有 \(n\) 个点的无向有标号连通图个数 . \((1 \le n \le 1.3 * 10^5)\) 题解 首先考虑 dp ... 直接算可行的方案数 , 容易算重复 . 我们用总方案数减 ...

  8. 洛谷P4841 城市规划(生成函数 多项式求逆)

    题意 链接 Sol Orz yyb 一开始想的是直接设\(f_i\)表示\(i\)个点的无向联通图个数,枚举最后一个联通块转移,发现有一种情况转移不到... 正解是先设\(g(n)\)表示\(n\)个 ...

  9. LOJ2527 HAOI2018 染色 容斥、生成函数、多项式求逆

    传送门 调了1h竟然是因为1004535809写成了998244353 "恰好有\(K\)种颜色出现了\(S\)次"的限制似乎并不容易达到,考虑容斥计算. 令\(c_j\)表示强制 ...

随机推荐

  1. Vue中table表头合并的用法

    <div class="panel-container"> <div> <table class="table-head" wid ...

  2. day31网络编程

    网络编程1. 目标:编写一个C/S架构的软件    C/S: Client(用户端)--------基于网络----------Server(服务端)    B/S: Browser-------基于 ...

  3. BZOJ 2173 luoguo P4451 [国家集训队]整数的lqp拆分

    整数的lqp拆分 [问题描述] lqp在为出题而烦恼,他完全没有头绪,好烦啊… 他首先想到了整数拆分.整数拆分是个很有趣的问题.给你一个正整数N,对于N的一个整数拆分就是满足任意m>0,a1 , ...

  4. DeepLearning初窥门径

    说明: 最近在看Ng的DL课程,感觉说的非常好,浅显易懂! 本来打算记录一下自己的学习过程,网上几个大神总结的太完美了,根本没必要自己去写了,而且浪费时间~~ 网易地址:http://mooc.stu ...

  5. File mapping

    文件映射的三个功能: 1.File mapping allows the process to use both random input and output (I/O) and sequentia ...

  6. leetcode42

    class Solution: def calLeft(self,height,rightval,left,right): if left>=right: return 0 sums = 0 m ...

  7. Linux命令:mapfile

    mapfile [-n 计数] [-O 起始序号] [-s 计数] [-t] [-u fd] [-C 回调] [-c 量子] [数组] 两个作用: 一是给index数组读取标准输入来赋值:二是文件描述 ...

  8. cmake find_package 中,include_directories,target_link_libraries 的值怎么知道?

    拿Sophus库为例: find_package(Sophus REQUIRED) include_directories(${Sophus_INCLUDE_DIRS}) target_link_li ...

  9. Python学习—基础篇之文件操作

    文件操作 文件操作也是编程中需要熟练掌握的技能,尤其是在后台接口编写和数据分析过程中,对各种类型的文件进行操作,获取文件信息或者对信息进行存储是十分重要的.本篇博客中将主要对常见的文本格式文件和Exc ...

  10. React createRef:引用

    一 代码 import React, { Component } from 'react'; class Box extends Component { render() { return <b ...