一开始以为暴力搞,后来看了数据范围还以为要FFT,各种被虐,然后Orz Seter大神!!!
我只想到了前三位:a * b <=> 10^(log(a) + log(b)),于是把乘的数都先log了最后再变回去就可以了。。。
然后后九位的方法:
Seter:"对于素数a,在N!中出现了N / a + N / a ^ 2+...次"
于是C(M, N) = N! / M! / (N - M)!就可以求每个质因数出现的次数,搞定!
/**************************************************************
Problem: 1300
User: rausen
Language: C++
Result: Accepted
Time:164 ms
Memory:1792 kb
****************************************************************/ #include <cstdio>
#include <cmath>
#include <algorithm> using namespace std;
typedef long long LL; int n, m;
LL ans = ;
bool p[]; int main(){
scanf("%d%d", &n, &m);
m = min(m, n - m); int i, j, x, y, c, s;
double k = , f = ;
for (i = ; i <= n; ++i)
if (!p[i]){
for (j = n, x = n - m, y = m; j;){
c = (j /= i) - (x /= i) - (y /= i), s = i;
for (c <<= ; c >>= ; s *= s)
if (c & ) ans *= s, ans %= (LL) 1e12;
}
if (i <= )
for (j = i * i; j <= n; j += i)
p[j] = ;
}
for (i = ; i <= m; ++i){
if (f > 1e7)
k += log10(f), f = ;
f = f * (n - m + i) / i;
}
k += log10(f); if (k < ) printf("%lld\n", ans);
else printf("%d...%09lld\n", (int) (pow(, k - floor(k) + ) + 1e-), ans % (int) 1e9);
return ;
}
p.s.话说Seter竟然可以C、C++混用,真是碉堡了。。。