LOJ 2409「THUPC 2017」小 L 的计算题 / Sum
思路
和玩游戏一题类似
定义\(A_k(x)=\sum_{i=0}^\infty a_k^ix^i=\frac{1}{1-a_kx}\)
用\(\ln 'x\)代替\(\frac{1}{x}\),
所以就是求
\]
这样没法快速计算
所以再设\(G(x)=\sum _{i=1}^n (ln(1-a_ix))'\)
所以
\]
所以
\]
\]
然后上分治+NTT就可以在\(O(n\log^2n)\)的时间内解决了
代码
#include <cstdio>
#include <algorithm>
#include <cstring>
#define int long long
using namespace std;
const int MOD = 998244353;
const int G = 3;
const int invG = 332748118;
const int MAXN = 2000000;
int pow(int a,int b){
int ans=1;
while(b){
if(b&1)
ans=(1LL*ans*a)%MOD;
a=(1LL*a*a)%MOD;
b>>=1;
}
return ans;
}
int rev[MAXN],inv_val[MAXN];
void cal_rev(int n,int lim){
for(int i=0;i<n;i++)
rev[i]=(rev[i>>1]>>1)|((i&1)<<(lim-1));
}
void cal_inv(int n){
inv_val[0]=0;
inv_val[1]=1;
for(int i=2;i<=n;i++)
inv_val[i]=(1LL*(MOD-MOD/i)*inv_val[MOD%i])%MOD;
}
void NTT(int *a,int opt,int n,int lim){
for(int i=0;i<n;i++)
if(i<rev[i])
swap(a[i],a[rev[i]]);
for(int i=2;i<=n;i<<=1){
int len=i/2,tmp=pow((opt)?G:invG,(MOD-1)/i);
for(int j=0;j<n;j+=i){
int arr=1;
for(int k=j;k<j+len;k++){
int t=(1LL*a[k+len]*arr)%MOD;
a[k+len]=(a[k]-t+MOD)%MOD;
a[k]=(a[k]+t)%MOD;
arr=(1LL*arr*tmp)%MOD;
}
}
}
if(!opt){
int invN = pow(n,MOD-2);
for(int i=0;i<n;i++)
a[i]=(1LL*a[i]*invN)%MOD;
}
}
void mul(int *a,int *b,int &at,int bt){
int midlen=1,midlim=0;
while((midlen)<(at+bt+2))
midlen<<=1,midlim++;
cal_rev(midlen,midlim);
static int tmp[MAXN];
for(int i=0;i<midlen;i++)
tmp[i]=b[i];
NTT(a,1,midlen,midlim);
NTT(tmp,1,midlen,midlim);
for(int i=0;i<midlen;i++)
a[i]=(1LL*a[i]*tmp[i])%MOD;
NTT(a,0,midlen,midlim);
at+=bt;
for(int i=0;i<=at;i++)
tmp[i]=0;
for(int i=at+1;i<midlen;i++)
a[i]=0,tmp[i]=0;
}
void inv(int *a,int *b,int dep,int &midlen,int &midlim){
if(dep==1){
b[0]=pow(a[0],MOD-2);
return;
}
inv(a,b,(dep+1)>>1,midlen,midlim);
static int tmp[MAXN];
while((dep<<1)>midlen)
midlen<<=1,midlim++;
for(int i=0;i<dep;i++)
tmp[i]=a[i];
for(int i=dep;i<midlen;i++)
tmp[i]=0;
cal_rev(midlen,midlim);
NTT(tmp,1,midlen,midlim);
NTT(b,1,midlen,midlim);
for(int i=0;i<midlen;i++)
b[i]=(1LL*b[i]*(2-1LL*tmp[i]*b[i]%MOD+MOD)%MOD)%MOD;
NTT(b,0,midlen,midlim);
for(int i=dep;i<midlen;i++)
b[i]=0;
}
void qd(int *a,int &t){
for(int i=0;i<t;i++)
a[i]=1LL*a[i+1]*(i+1)%MOD;
a[t]=0;
t--;
}
void jf(int *a,int &t){
t++;
for(int i=t;i>0;i--){
a[i]=1LL*a[i-1]*inv_val[i]%MOD;
}
a[0]=0;
}
void ln(int *a,int *b,int n){
static int tmp[MAXN];
for(int i=0;i<n;i++)
tmp[i]=0,b[i]=a[i];
int midlen=1,midlim=0;
inv(a,tmp,n,midlen,midlim);
int t=n;
qd(b,t);
mul(b,tmp,t,n);
jf(b,t);
for(int i=0;i<n;i++)
tmp[i]=0;
for(int i=n;i<midlen;i++)
b[i]=tmp[i]=0;
}
int n,a[MAXN],b[MAXN],c[MAXN];
int cnt=0,barrel[40][MAXN];
void solve(int l,int r,int *val,int &len){
if(l==r){
val[0]=1;
val[1]=MOD-a[l];
len=1;
return;
}
int *la=barrel[cnt++],*ra=barrel[cnt++];
int num=cnt,lenl,lenr;
// printf("num=%lld\n",cnt);
int mid=(l+r)>>1;
solve(l,mid,la,lenl);
solve(mid+1,r,ra,lenr);
int midlen=1,midlim=0;
while(midlen<(lenl+lenr+2))
midlen<<=1,midlim++;
cal_rev(midlen,midlim);
NTT(la,1,midlen,midlim);
NTT(ra,1,midlen,midlim);
for(int i=0;i<midlen;i++)
val[i]=(1LL*la[i]*ra[i])%MOD;
NTT(val,0,midlen,midlim);
for(int i=0;i<midlen;i++)
la[i]=ra[i]=0;
len=lenl+lenr;
cnt=num-2;
}
signed main(){
// freopen("test.in","r",stdin);
// freopen("test.out","w",stdout);
int T;
scanf("%lld",&T);
while(T--){
// printf("ok %lld\n",T);
scanf("%lld",&n);
cal_inv(n+10);
for(int i=1;i<=n;i++)
scanf("%lld",&a[i]);
int len=0;
solve(1,n,b,len);
// for(int i=0;i<=n;i++)
// printf("%lld ",b[i]);
// printf("\n");
ln(b,c,n+1);
// for(int i=0;i<=n;i++)
// printf("%lld ",c[i]);
// printf("\n");
int t=n;
qd(c,t);
for(int i=n;i>=1;i--)
c[i]=MOD-c[i-1];
c[0]=n;
// for(int i=0;i<=n;i++)
// printf("%lld ",c[i]);
// printf("\n");
// for(int i=1;i<=n;i++)
// printf("%lld ",c[i]);
// printf("\n");
int ans=0;
for(int i=1;i<=n;i++)
ans^=c[i];
printf("%lld\n",ans);
for(int i=0;i<=n;i++)
a[i]=b[i]=c[i]=0;
}
return 0;
}
LOJ 2409「THUPC 2017」小 L 的计算题 / Sum的更多相关文章
- LOJ#2409. 「THUPC 2017」小 L 的计算题 / Sum(生成函数)
题意 给定一个长为 \(n\) 的序列 \(\{a_i\}\) 对于 \(k \in [1, n]\) 求 \[ f_k = \sum_{i = 1}^{n} a_i^k \pmod {9982443 ...
- 题解 「THUPC 2017」小 L 的计算题 / Sum
题目传送门 题目大意 给出 \(a_{1,2,...,n}\),对于 \(\forall k\in [1,n]\) ,求出: \[\sum_{i=1}^{n}a_i^k \] \(n\le 2\tim ...
- LOJ 2288「THUWC 2017」大葱的神力
LOJ 2288「THUWC 2017」大葱的神力 Link Solution 比较水的提交答案题了吧 第一个点爆搜 第二个点爆搜+剪枝,我的剪枝就是先算出 \(mx[i]\) 表示选取第 \(i \ ...
- @loj - 2977@ 「THUSCH 2017」巧克力
目录 @description@ @solution@ @accepted code@ @details@ @description@ 「人生就像一盒巧克力,你永远不知道吃到的下一块是什么味道.」 明 ...
- LOJ #2978「THUSCH 2017」杜老师
听说LOJ传了THUSC题赶紧上去看一波 随便点了一题都不会做想了好久才会写暴力爆了一发过了... LOJ #2978 题意 $ T$次询问,每次询问$ L,R$,问有多少种选取区间中数的方案使得选出 ...
- LOJ 2980 「THUSCH 2017」大魔法师——线段树
题目:https://loj.ac/problem/2980 线段树维护矩阵. 然后是 30 分.似乎是被卡常了?…… #include<cstdio> #include<cstri ...
- LOJ 2979 「THUSCH 2017」换桌——多路增广费用流
题目:https://loj.ac/problem/2979 原来的思路: 优化连边.一看就是同一个桌子相邻座位之间连边.相邻桌子对应座位之间连边. 每个座位向它所属的桌子连边.然后每个人建一个点,向 ...
- LOJ 2978 「THUSCH 2017」杜老师——bitset+线性基+结论
题目:https://loj.ac/problem/2978 题解:https://www.cnblogs.com/Paul-Guderian/p/10248782.html 第 i 个数的 bits ...
- LOJ 2997 「THUSCH 2017」巧克力——思路+随机化+斯坦纳树
题目:https://loj.ac/problem/2977 想到斯坦纳树.但以为只能做 “包含一些点” 而不是 “包含一些颜色” .而且不太会处理中位数. 其实 “包含一些颜色” 用斯坦纳树做也和普 ...
随机推荐
- ASCII Unicode UTF-8 之间的关系
转载请标明:https://i.cnblogs.com/EditPosts.aspx?opt=1 1. ASCII ASCII 只有127个字符,表示英文字母的大小写.数字和一些符号,但由于其他语言用 ...
- spring-data-redis2.0以上配置redis连接
1.所需的maven依赖 <dependency> <groupId>org.springframework.data</groupId> <artifact ...
- jacoco统计server端功能测试覆盖率
jacoco可以统计,功能测试时,server代码调用的覆盖情况.这里对服务器端的java代码进行统计. 操作步骤如下: 第一步:更改server的启动脚本,使用jacocoagent.jar ...
- 循环结构 : while .. for
# ###循环结构 : while .. for ''' while 循环 可以提高代码的效率,减少代码的冗余 while 条件表达式: code1 code2 如果条件表达式成立,返回True , ...
- 使用@Autowired时,取值为null
如果取不到,可以考虑其他方式 场景: @Autowired private StringRedisTemplate redisTemplate; 想使用redisTemplate,但是使用时为null ...
- P3803 【模板】多项式乘法(FFT)
传送门: 参考博客 1:大佬 attack 参考博客 2:大佬 胡小兔 在这里再膜拜一下这两位大佬 Orz%%% #include<iostream> #include<cstd ...
- git加速和只下载部分目录
浅复制 工作要用到的.git有1.8G太大了.下载过程要好几个小时,太慢了.可以这样操作 git clone 默认会下载项目的完整历史版本,如果你只关心最新版的代码,而不关心之前的历史信息,可以使用 ...
- PHP----------线程安全和非线程安全的介绍
1.Linux下的PHP,没有线程安全版和非线程安全版之分.从2000年10月20日发布的第一个Windows版的PHP3.0.17开始的都是线程安全的版本,直至5.2.1版本开始有Thread Sa ...
- UGUI背包系统
在Unity3d中,UGUI提供了Scroll Rect.Grid Layout Group.Mask这三个组件,下面就给大家介绍下如何用这个三个组件来实现滚动视图. 首先放置好背包的背景图 在矩形线 ...
- php魔术变量以及命名空间
魔术变量: PHP 向它运行的任何脚本提供了大量的预定义常量. 不过很多常量都是由不同的扩展库定义的,只有在加载了这些扩展库时才会出现,或者动态加载后,或者在编译时已经包括进去了. 有八个魔术常量它们 ...