[数论][Lucas定理] 51Nod 1778 小Q的集合
2018-02-25 11:37
316 查看
SolutionSolution
通过枚举⌊im⌋⌊im⌋和imodmimodm。由Lucas定理是可以独立计算贡献的。ans====∑i=0n(ni)(ik−(n−i)k)2∑i=0n(⌊nm⌋⌊im⌋)(nmodmimodm)(ik−(n−i)k)2∑i=0⌊nm⌋(⌊nm⌋i)∑j=0nmodm(nmodmj)(ik−(n−i)k)22⌊nm⌋∑j=0nmodm(nmodmj)(ik−(n−i)k)2(57)(58)(59)(60)(57)ans=∑i=0n(ni)(ik−(n−i)k)2(58)=∑i=0n(⌊nm⌋⌊im⌋)(nmodmimodm)(ik−(n−i)k)2(59)=∑i=0⌊nm⌋(⌊nm⌋i)∑j=0nmodm(nmodmj)(ik−(n−i)k)2(60)=2⌊nm⌋∑j=0nmodm(nmodmj)(ik−(n−i)k)2预处理ikik好像可以O(m)O(m)做:
首先这是一个完全积性函数,可以线性筛,需要计算快速幂的只有质数,大概复杂度就是O(π(m)logm)O(π(m)logm),就O(m)O(m)了。
#include <bits/stdc++.h> #define show(x) cerr << #x << " = " << x << endl using namespace std; typedef long long ll; typedef pair<int, int> Pairs; const int N = 1010101; inline char get(void) { static char buf[100000], *S = buf, *T = buf; if (S == T) { T = (S = buf) + fread(buf, 1, 100000, stdin); if (S == T) return EOF; } return *S++; } template<typename T> inline void read(T &x) { static char c; x = 0; int sgn = 0; for (c = get(); c < '0' || c > '9'; c = get()) if (c == '-') sgn = 1; for (; c >= '0' && c <= '9'; c = get()) x = x * 10 + c - '0'; if (sgn) x = -x; } int k, m, pcnt; int fk , n ; int ndm, nmm, ans; int s ; int fac , ifac ; int prime , vis ; inline int pwr(int a, int b) { int c = 1; while (b) { if (b & 1) c = (ll)c * a % m; b >>= 1; a = (ll)a * a % m; } return c; } inline void pre(int n) { ifac[1] = 1; for (int i = 2; i <= n; i++) ifac[i] = (ll)(m - m / i) * ifac[m % i] % m; fac[0] = ifac[0] = 1; for (int i = 1; i <= n; i++) { fac[i] = (ll)fac[i - 1] * i % m; ifac[i] = (ll)ifac[i - 1] * ifac[i] % m; } fk[1] = 1; for (int i = 2; i <= n; i++) { if (!vis[i]) { prime[++pcnt] = i; fk[i] = pwr(i, k); } for (int j = 1, x; j <= pcnt && (x = prime[j] * i) <= n; j++) { vis[x] = 1; fk[x] = (ll)fk[i] * fk[prime[j]] % m; if (i % prime[j] == 0) break; } } } inline int C(int n, int m) { return (ll)fac * ifac[m] % ::m * ifac[n - m] % ::m; } int main(void) { freopen("1.in", "r", stdin); freopen("1.out", "w", stdout); for (char c = get(); c >= '0' && c <= '9'; c = get()) n[++*n] = c - '0'; read(k); read(m); pre(m - 1); int cur = 0; for (int i = 1; i <= *n; i++) { cur = cur * 10 + n[i]; nmm = (nmm * 10 + n[i]) % m; ndm = (ndm * 10 + cur / m) % (m - 1); cur %= m; } for (int j = 0; j <= nmm; j++) ans = (ans + (ll)pwr(fk[j] - fk[nmm - j] + m, 2) * C(nmm, j) % m) % m; ans = (ll)ans * pwr(2, ndm) % m; cout << ans << endl; return 0; }
相关文章推荐
- [方差+lucas定理] 51nod 算法马拉松25 D. 小Q的集合
- 51Nod 1778 小Q的集合
- [期望与方差 Lucas定理] 算法马拉松25 D 小Q的集合
- 51Nod1778 小Q的集合 【组合数】【Lucas定理】
- 51nod 1120 机器人走方格 V3 卡特兰数 lucas定理
- 教你一波Lucas(卢卡斯)定理在数论解题中的应用
- 数论 Lucas定理 hdu3037
- 51nod 1616 最小集合(数论)(枚举)
- 【数论】Lucas定理与拓展Lucas定理
- BZOJ 1951 SDOI2010 古代猪文 数论 Lucas定理
- 51nod 1616 最小集合(数论)(枚举)
- 【BZOJ 1951】 [Sdoi2010]古代猪文|数论|中国剩余定理|Lucas
- 51nod 1616 最小集合(数论)(枚举)
- ACM 数论 HDU 3037 Saving Beans (Lucas定理,大数取魔)
- 【数论】Lucas定理推导过程+模板
- 数论-Lucas(卢卡斯定理)
- 51nod 1616 最小集合(数论)(枚举)
- Bzoj 4403: 序列统计 Lucas定理,组合数学,数论
- lucas (数论定理)
- CPC23-4-K. 喵喵的神数 (数论 Lucas定理)