多项式
c20181220_xiang_m_y
学习他人,提升自己;
提升自己,帮助他人。
展开
-
UOJ#551【UNR #4】校园闲逛 【多项式好题(雾)】
题目描述Linknnn 个点 mmm 条边的无向带权图,1≤1\le1≤ 边权 ≤V\le V≤V,QQQ 次询问,从 xxx 到 yyy 有多少条权值和为 www 的路径,路径可以经过重复边和重复点。1≤n≤8,0≤m≤3∗105,1≤V≤65000,0≤Q≤100001\le n\le 8, 0\le m\le 3*10^5,1\le V\le65000,0\le Q\le 100001≤n≤8,0≤m≤3∗105,1≤V≤65000,0≤Q≤10000题目分析暴力 DP:dp[i][j][v原创 2020-08-13 08:01:42 · 367 阅读 · 0 评论 -
20200802 T3 我永远喜欢【生成函数容斥,拉格朗日反演】
题目描述有 nnn 种颜色的石子,每种 cic_ici 个,记一个石子序列首尾相接后极长连续段的长度为 lil_ili,求所有石子序列的 1∏li!\frac 1{\prod l_i!}∏li!1 的和。n≤105,∑ci≤2∗105n\le10^5,\sum c_i\le2*10^5n≤105,∑ci≤2∗105题目分析先考虑去掉首尾相接的情况怎么做。因为限制了极长,容易想到把每种颜色分成几段,然后合并,但是不好保证相同颜色不被并在一起。先不谈容斥做法,题解给出了一种用生成函数解决的原创 2020-08-09 20:49:41 · 304 阅读 · 0 评论 -
20200803 T3 我的朋友们【分治NTT优化递推】
题意描述Lynn 要买新发卡,拉上了君儿和薇帮她参考。可供选择的发卡总共有 nnn 个。Lynn 会问很多次君儿某个发卡是否适合 Lynn。由于君儿沉迷于她的数竞习题,所以某一次 Lynn 问她第 iii 个发卡是否适合 Lynn 的时候,她有 1−pi(0<pi<1)1-p_i(0<p_i<1)1−pi(0<pi<1) 的概率会告诉 Lynn “这个发卡很适合你”,有 pip_ipi 的概率告诉 Lynn “这个发卡一点也不适合你”。因为君儿太沉迷了,所以两原创 2020-08-03 23:01:31 · 423 阅读 · 0 评论 -
20200713 T1 集合划分(JOISCD1T1)【多项式优化DP】
题目描述长度为 2n2n2n 的 {Ai},{Bi}\{A_i\},\{B_i\}{Ai},{Bi} 序列,现在要选出一个新序列 {Ci}\{C_i\}{Ci},其中 Ci=AiC_i=A_iCi=Ai 或 Ci=BiC_i=B_iCi=Bi要求 Ci≤Ci+1C_i\le C_{i+1}Ci≤Ci+1,且选择的 AAA 的个数恰为 nnn 个,问方案数,方案不同当且仅当某个位置的选择不同。n≤5∗104,mod 998244353n\le 5*10^4,\mod 99824435原创 2020-07-22 16:40:37 · 263 阅读 · 0 评论 -
20200721 T3 s3mple【DP,点值优化卷积】
题目描述n≤200,x≤109n\le200,x\le 10^9n≤200,x≤109题目分析容易想到枚举最大值的位置,然后递归两边,于是有这样的DP:上面式子里面有个地方要改成 f(n−1−i,m−min(i,n−1−i)−1−j)f(n-1-i,m-min(i,n-1-i)-1-j)f(n−1−i,m−min(i,n−1−i)−1−j)可以DP求出区间长度为 nnn 的 xxx 的最大值:mx[n]=min(n+12,n+1−n+12)mx[n]=min(\frac {n+1}2,n+1原创 2020-07-21 20:05:23 · 251 阅读 · 0 评论 -
20200630 多项式难题训练
LOJ#6247. 九个太阳(单位根反演)给定 n≤1015,k≤220n\le10^{15},k\le2^{20}n≤1015,k≤220,kkk 是 2 的幂,求:∑k∣i,0≤i≤n(ni)\sum_{k|i,0\le i\le n}\binom nik∣i,0≤i≤n∑(in) 单位根反演相当于就是求 f(x)=∑i=0n(ni)xi=(x+1)nf(x)=\sum_{i=0}^n \binom nix^i=(x+1)^nf(x)=∑i=0n(in)xi=(x+1)n 的 kkk 的倍原创 2020-07-02 14:49:57 · 2418 阅读 · 0 评论 -
LOJ#2320. 「清华集训 2017」生成树计数
题目描述:图中有 nnn 个连通块,每个连通块有 aia_iai 个点,需要再连 n−1n-1n−1 条边,使其变成一棵树。对一种连边方案,设原图中第 iii 个连通块连出了 did_idi 条边,那么这棵树 TTT 的价值为val(T)=(∏i=1ndim)∗(∑i=1ndim)\text{val}(T)=\left(\prod_{i=1}^nd_i^m\right)*\left(\sum_{i=1}^nd_i^m\right)val(T)=(i=1∏ndim)∗(i=1∑ndim)求所原创 2020-06-22 20:07:41 · 171 阅读 · 0 评论 -
多项式运算
多项式乘法多项式求逆多项式开根前置知识:多项式牛顿迭代:用于寻找多项式零点,已知多项式GGG,求满足G(F(x))≡0(modxn)G(F(x))\equiv0\pmod {x^n}G(F(x))≡0(modxn)的F(x)F(x)F(x)假设已经求出了G(F0(x))≡0(modxn2)G(F_0(x))\equiv0\pmod {x^{\frac n2}}G(F0(x))≡0(modx2n)将G(F(x))G(F(x))G(F(x))在F0(x)F_0(x)F0(x)处泰勒展开:G原创 2020-06-22 20:14:42 · 1272 阅读 · 0 评论 -
FWT学习笔记
拓展:给出序列a1,...,ana_1,...,a_na1,...,an以及ppp,求i∈[1,n]i\in[1,n]i∈[1,n]的(1+pxai)(1+px^{a_i})(1+pxai)的异或卷积。分治FWT复杂度太高。考虑(1+pxa)(1+px^{a})(1+pxa)的对应的多项式bixib_ix^ibixi的FWTFWTFWT变换:ci=∑j(−1)d(i&j)bjc_i=\sum\limits_{j}(-1)^{d(i\&j)}b_jci=j∑(−1)d(i&a原创 2020-05-14 22:10:56 · 685 阅读 · 0 评论 -
20200501 zr T1【k次多项式的n次幂的前x项系数】
题目描述:题目分析:第一种是链式法则,第二种是(uv)′=u′v+uv′(uv)'=u'v+uv'(uv)′=u′v+uv′令fi=Pn(x)[i]f_i=P^n(x)[i]fi=Pn(x)[i],gi=P(x)[i]g_i=P(x)[i]gi=P(x)[i],已知ggg且最高项次数为kkk,我们要求fff的前XXX项首先将第二个式子的第二项移到左边,可以得到:nPn(x)P′...原创 2020-05-05 15:46:25 · 612 阅读 · 0 评论 -
20200330 T1 烷烃计数
题目描述:计算化学式为CnH2n+2C_nH_{2n+2}CnH2n+2的烷烃的同分异构体个数。等价于求 nnn 个点的无标号无根树并且满足每个点的度数 ≤4\le4≤4 的树的个数。多组数据,对998244353取模。题解:1. T=1, n≤20001.~~~T=1,~n\le20001. T=1,&n...原创 2020-04-04 12:15:03 · 998 阅读 · 1 评论 -
模拟赛20200221(zjx)【带删除线性基(线段树分治),不经过某条边的最短路(最短路树),嵌套for循环循环次数(拉格朗日插值)】
T1:题意:给出mmm个值在[0,232)[0,2^{32})[0,232)内的数,kkk次操作添加一个数或删除一个数,或询问某个xxx是否能被这些数异或得到。m,k≤105m,k\le10^5m,k≤105题解:线性基套个线段树分治即可。网上好像有基于线性基性质的神奇做法,不过复杂度似乎都是两个log。Code:#include<bits/stdc++.h>#defi...原创 2020-02-21 22:43:43 · 442 阅读 · 1 评论 -
BZOJ3601 一个人的数论【莫比乌斯反演+拉格朗日插值求自然数幂和多项式系数】
题目描述:题目分析:这似乎是官方题解:题解1这是比较详细的题解:题解2拉格朗日插值求系数调错调了1h+…#include<bits/stdc++.h>#define maxn 105using namespace std;const int mod = 1e9+7;int d,n,a[maxn],s[maxn],f[maxn],len;inline int Pow...原创 2019-10-19 17:31:08 · 157 阅读 · 0 评论 -
BZOJ2655: calc【容斥/插值】
题目描述:一个序列a1,…,an是合法的,当且仅当:长度为给定的n。a1,…,an都是[1,A]中的整数。a1,…,an互不相等。一个序列的值定义为它里面所有数的乘积,即a1a2…an。求所有不同合法序列的值的和。两个序列不同当且仅当他们任意一位不一样。输出答案对一个数mod取余的结果。1≤n≤A≤1091\le n\le A\le10^91≤n≤A≤109,A<mo...原创 2019-09-16 23:07:46 · 203 阅读 · 0 评论 -
拉格朗日插值模板
贴dalao博客学习一波顺便打个板:#include<cstdio>#include<cctype>#define maxn 2005inline void read(int &a){ char c;bool f=0; while(!isdigit(c=getchar())) if(c=='-') f=1; for(a=c-'0';...原创 2019-03-18 09:56:31 · 122 阅读 · 0 评论 -
51nod 1850 抽卡大赛 (十二省联考模测) 【"动态"(多项式) DP】
题目传送门题目分析:如果只求第i个人的收入,那么只需要枚举他选抽到哪一张卡,然后O(n2)DP算出他排名为k的概率即可。求一个人是O(n3),求n个人就是O(n4)的。打完这个暴力发现居然有63分然后就溜了。。O(n3)的做法:对于一个权值A,显然只有大于它的值会对它的排名造成影响。所以我们将所有的A值放在一起排序,考虑如何求出Ai的期望排名。如果不同的A由不同的人选择,那么这是个很...原创 2019-03-24 21:00:04 · 171 阅读 · 0 评论 -
多项式求逆元模板【NTT/拆系数FFT】
一. 洛谷4238:模998244353,NTT:#include&amp;amp;lt;cstdio&amp;amp;gt;#include&amp;amp;lt;algorithm&amp;amp;gt;#define maxn 400005#define LL long longusing namespace std;const int mod = 998244353, G = 3;LL wn,w;inline LL k原创 2018-12-28 15:00:09 · 424 阅读 · 0 评论