結果
問題 |
No.3030 Kruskal-Katona
|
ユーザー |
![]() |
提出日時 | 2025-02-21 22:57:08 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 44 ms / 2,000 ms |
コード長 | 2,652 bytes |
コンパイル時間 | 2,262 ms |
コンパイル使用メモリ | 199,716 KB |
実行使用メモリ | 18,924 KB |
最終ジャッジ日時 | 2025-02-21 22:57:13 |
合計ジャッジ時間 | 4,640 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 27 |
ソースコード
#include<bits/stdc++.h> // #include<atcoder/all> // #include<boost/multiprecision/cpp_int.hpp> using namespace std; // using namespace atcoder; // using bint = boost::multiprecision::cpp_int; using ll = long long; using ull = unsigned long long; using ld = long double; using P = pair<ll,ll>; using vi = vector<ll>; using vvi = vector<vi>; using vvvi = vector<vvi>; using ve = vector<vector<int>>; using vb = vector<bool>; using vvb = vector<vb>; #define rep(i,s,n) for(ll i = (ll)s;i < (ll)n;i++) #define rrep(i,l,r) for(ll i = r-1;i >= l;i--) #define ALL(x) (x).begin(),(x).end() #define sz(c) ((ll)(c).size()) #define LB(A,x) (int)(lower_bound(A.begin(),A.end(),x)-A.begin()) #define UB(A,x) (int)(upper_bound(A.begin(),A.end(),x)-A.begin()) // #define MOD 1000000007 #define MOD 998244353 template<typename T>using min_priority_queue=priority_queue<T,vector<T>,greater<T>>; template<typename T>ostream&operator<<(ostream&os,vector<T>&v){for(int i = 0;i < v.size();i++)os<<v[i]<<(i+1!=v.size()?" ":"");return os;} template<typename T>istream&operator>>(istream&is,vector<T>&v){for(T&in:v)is>>in;return is;} template<typename T1,typename T2>ostream&operator<<(ostream&os,pair<T1,T2>&p){os<<p.first<<" "<<p.second;return os;} template<typename T1,typename T2>istream&operator>>(istream&is,pair<T1,T2>&p){is>>p.first>>p.second;return is;} template<typename T> inline bool chmax(T &a,T b){if(a < b){a = b;return true;}return false;} template<typename T> inline bool chmin(T &a,T b){if(a > b){a = b;return true;}return false;} ld dist(ld x1,ld y1,ld x2, ld y2){return sqrtl((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2));} class Comb{ public: const int n; // const int mod = 1000000007; const int mod = 998244353; vector<long long> fac,refac; Comb(int N) : n(N),fac(N+1,1),refac(N+1,1) { for(int i = 0;i < n;i++)fac[i+1] = fac[i]*(i+1)%mod; refac[n] = modpow(fac[n],mod-2); for(int i = n;i > 0;i--)refac[i-1] = refac[i]*i%mod; } long long comb(long long n,long long r){ if(n < r)return 0; long long res = fac[n]*refac[r]%mod; res = res*refac[n-r]%mod; return res; } long long modpow(long long n,long long r){ long long res = 1; while(r){ if(r & 1)res = res*n%mod; n = n*n%mod; r >>= 1; } return res; } }; int main(){ ios_base::sync_with_stdio(0), cin.tie(0); ll n,i;cin >> n >> i; Comb cb(1000000); vi res; while(n){ if(i == 1){ res.emplace_back(n);break; } rep(k,i,1000000){ if(cb.comb(k,i) > n){ res.emplace_back(k-1); n -= cb.comb(k-1,i); i--; break; } } } cout << res << endl; return 0; }