高精度求幂 public static char[] exponentiation(string a,int r) { ]; string b = ""; string c = a; ; i < r-; i++) { aa = acm.Quadrature(c, a); b = ""; foreach (var item in aa) { b += item; } c = b; } return aa; }…

A * B Problem Plus 链接:http://acm.hdu.edu.cn/showproblem.php?pid=1402 Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Submission(s): 26449    Accepted Submission(s): 6917 Problem Description Calculate A * B.  …

用递推的方式写的写挂了,而如果用组合数又不会高精度除法,偶然看到了别人的只用高精度乘低精度求组合数的方法,记录一下 #include<bits/stdc++.h> using namespace std; const int maxn=60010; const long long M=1000000000; typedef long long LL; LL num[maxn]; int cnt[maxn*2];//记录素数因子的个数 int len=1; void get_cnt(int x,…

一.多项式求逆 给定一个多项式 \(F(x)\),请求出一个多项式 \(G(x)\), 满足 \(F(x) * G(x) \equiv 1 ( \mathrm{mod\:} x^n )\).系数对 \(998244353\)取模. 考虑递归求解,当\(F\)的最高次为\(0\)时,\(G_0=F_0^{-1}\) 假设我们知道了\(F(x)\)在模\(x^{\left \lceil \frac{n}{2}\right \rceil}\)意义下的逆元\(G'\) 那么\(F∗G′≡1(\mathr…

多项式 代码 const int nsz=(int)4e5+50; const ll nmod=998244353,g=3,ginv=332748118ll; //basic math ll qp(ll a,ll b){ ll res=1; for(;b;a=a*a%nmod,b>>=1)if(b&1)res=res*a%nmod; return res; } ll inv(ll n){ return qp(n,nmod-2); } //polynomial operations //…

https://nanti.jisuanke.com/t/31716 题意 n颗糖果n个人,按顺序给每个人任意数目(至少一个)糖果,问分配方案有多少. 分析 插板法或者暴力打表后发现答案就为2^(n-1),只是这个n有点大.于是马上用java.然而现实相当残酷,超时. 然后想到降幂,即(a^b)%m=a^(b%phi(m))%m,当gcd(b,m)==1.这里显然互质,于是降幂后仍然用java写,还是tle. 而后还尝试了C++大数来写,可能是使用姿势错误,也t了. 到了最后一小时,没错,我们队…

城市规划 Time Limit: 40 Sec  Memory Limit: 256 MBSubmit: 1091  Solved: 629[Submit][Status][Discuss] Description 刚刚解决完电力网络的问题, 阿狸又被领导的任务给难住了. 刚才说过, 阿狸的国家有n个城市, 现在国家需要在某些城市对之间建立一些贸易路线, 使得整个国家的任意两个城市都直接或间接的连通. 为了省钱, 每两个城市之间最多只能有一条直接的贸易路径. 对于两个建立路线的方案, 如果存在一…

时间限制: 1 s  空间限制: 256000 KB  题目等级 : 白银 Silver 题目描述 Description 用高精度计算出S=n! 其中"!"表示阶乘,例如:5!=5*4*3*2*1 输入描述 Input Description 输入正整数N 输出描述 Output Description 输出计算结果S. 样例输入 Sample Input 3 样例输出 Sample Output 6 数据范围及提示 Data Size & Hint n<=100 这题…

hdu 5730 Shell Necklace 题意:求递推式\(f_n = \sum_{i=1}^n a_i f_{n-i}\),模313 多么优秀的模板题 可以用分治fft,也可以多项式求逆 分治fft 注意过程中把r-l+1当做次数界就可以了,因为其中一个向量是[l,mid],我们只需要[mid+1,r]的结果. 多项式求逆 变成了 \[ A(x) = \frac{f_0}{1-B(x)} \] 的形式 要用拆系数fft,直接把之前的代码复制上就可以啦 #include <iostream…

    关于FFT,咱们都会迫不及待地 @  .....(大雾)(貌似被玩坏了...)    .....0.0学习FFT前先orz FFT君.         首先先是更详细的链接(手写版题解点赞0v0) FFT的资料    其实众所周知的最详细的算法解释在<算法导论>上...然后咱就是边看着那个边码理解的...    首先来看看多项式乘法和快速FFT的关系,然后咱们再来看能否聊到卷积什么的东西...    其实觉得还是去看算法导论最好.   [一.多项式及其表达方式.]      首先什么是…