#include #include #include using namespace std; using i64 = long long; using u64 = unsigned long long; using i32 = int; using u32 = unsigned int; #define rep(i,n) for(int i=0; i<(n); i++) vector> factorize(u64 n){ vector> res; for(u64 i=2; i*i<=n; i++){ if(n%i) continue; res.push_back({i,0}); while(n%i==0){ n/=i; res.back().second++; } } if(n!=1) res.push_back({n,1}); return res; } vector divs(vector> F){ vector res = {1}; for(auto f:F){ int t=res.size() * f.second; rep(i,t) res.push_back(res[i] * f.first); } return move(res); } vector> gcd_frequency(u64 n){ vector> ans; ans.push_back(make_pair((u64)1, (u64)1)); for(auto f : factorize(n)){ vector> buf; u64 multiplier = 1; u64 inv_multiplier = 1; for(u64 p_idx=0; p_idx> N >> K >> M; auto F = factorize(M); u64 ans = 5001001001001001001; for(auto[ p,e ] : F){ u64 fpres = fp(N,p) - fp(K,p) - fp(N-K,p); fpres /= e; ans = min(ans, fpres); } cout << ans << endl; return 0; } struct ios_do_not_sync{ ios_do_not_sync(){ std::ios::sync_with_stdio(false); std::cin.tie(nullptr); } } ios_do_not_sync_instance;