BraketBN

【题目链接】 同【BZOJ3260的题解】 /* Pigonometry */ #include #include using namespace std; typedef unsigned long long ULL; typedef long long LL; const ULL p = 1000000007; void exgcd(ULL a, ULL b, LL &x,

【题目链接】 公式为： 设 ans = ∑(n! / ((d!)^(n/d)*(n/d)!)) 则答案为m ^ ans 证明： 考虑现在有d * k个点，d代表每个团的点数，那么k就是个数了，记方案数为Ak。 然后现在又来了d个点，记方案数为Ak+1。（即现在有n = d * (k + 1)个点） 我们选择一个点，让这个点与其他的d - 1个点组成团，方案数为C(d *

【题目链接】 答案为C(n + R - L + 1, n) - 1。 大体推导过程是固定长度，然后隔板法，然后对长度求和。 暴力算组合数那里，要判断n /* Pigonometry */ #include #include #include using namespace std; typedef unsigned long long ULL; const int

【题目链接】 【thy_asdf的题解】 /* Telekinetic Forest Guard */ #include #include #include using namespace std; typedef long long LL; LL n, k, m, p; inline LL qpow(LL a, LL n) { LL res = 1; for(; n

因为新数列中的数线性无关，所以答案就是C(n, k) % p，因为p比较小，写个扩展Lucas就好了。 证明一下线性无关： 首先有两个引理： 定理1：一个n×n的矩阵A是非奇异的充要条件为 A的行列式不等于0。 定理2：令x1, x2, ..., xn为R^n中的n个向量，并令X = (x1, ..., xn)。向量x1, x2, ..., xn线性无关的充要条件是X

【题目链接】 原来那个板有问题...重新写了下。 把分解质因数写挂了，WA了一次... /* Telekinetic Forest Guard */ #include #include #include using namespace std; typedef long long LL; const int maxn = 10; LL w[maxn]; inlin

【题目链接】 太神了【将狼踩尽的题解】 /* Telekinetic Forest Guard */ #include #include using namespace std; typedef long long LL; const int maxn = 2000001, p = 20100403; int n, m, fact[maxn]; inline int mul(

【题目链接】 还是python好用 【大爷的题解】 def A(n): res = 1 for i in range(1, n + 1): res *= i return res def C(n, m): if n < m: return 0 return A(n) // A(m) // A(n - m) n, m = [int(i) for i in raw_inpu

【题目链接】 同【CodeChef-CSEQ的题解】 /* Telekinetic Forest Guard */ #include #include #include using namespace std; typedef unsigned long long ULL; const int maxn = 1000505; const ULL p = 1000003; ULL

设m 所以求出组合数就好了... 发现n * m最大为10^12，那么m最大为10^6。 将组合数用阶乘形式展开，为 (n + m + 1)! / n! / (m + 1)!。 将前两项约掉，得到 (n + 1) * (n + 2) * ... * (n + m + 1) / (m + 1)! 分子有m项，分母有m项，所以可以O(m)求出来。 因为有inv(

【题目链接】 定理题... /* Think Thank Thunk */ #include #include #include using namespace std; typedef long long LL; const int maxn = 105, p = 1000003; int n, fact[200005], A[maxn][maxn], du[maxn];