結果
問題 | No.737 PopCount |
ユーザー | hayaten |
提出日時 | 2022-09-23 23:14:27 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 2 ms / 1,000 ms |
コード長 | 3,121 bytes |
コンパイル時間 | 5,964 ms |
コンパイル使用メモリ | 287,904 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2024-12-22 05:43:02 |
合計ジャッジ時間 | 6,254 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,820 KB |
testcase_01 | AC | 2 ms
6,820 KB |
testcase_02 | AC | 2 ms
6,816 KB |
testcase_03 | AC | 2 ms
6,820 KB |
testcase_04 | AC | 2 ms
6,820 KB |
testcase_05 | AC | 2 ms
6,816 KB |
testcase_06 | AC | 2 ms
6,820 KB |
testcase_07 | AC | 2 ms
6,816 KB |
testcase_08 | AC | 2 ms
6,816 KB |
testcase_09 | AC | 2 ms
6,820 KB |
testcase_10 | AC | 2 ms
6,816 KB |
testcase_11 | AC | 2 ms
6,816 KB |
testcase_12 | AC | 2 ms
6,820 KB |
testcase_13 | AC | 2 ms
6,820 KB |
testcase_14 | AC | 2 ms
6,820 KB |
testcase_15 | AC | 2 ms
6,816 KB |
testcase_16 | AC | 2 ms
6,816 KB |
testcase_17 | AC | 2 ms
6,820 KB |
ソースコード
#include <bits/stdc++.h> #include <atcoder/all> #include <ext/pb_ds/assoc_container.hpp> #include <ext/pb_ds/tree_policy.hpp> using namespace std; using namespace __gnu_pbds; template <typename T> using ordered_set = tree<T, null_type, less<T>, rb_tree_tag, tree_order_statistics_node_update>; struct __INIT{__INIT(){cin.tie(0);ios::sync_with_stdio(false);cout<<fixed<<setprecision(15);}} __init; #define REP(i,n) for(ll i=0;i<ll(n);i++) #define REPD(i,n) for(ll i=n-1;i>=0;i--) #define FOR(i,a,b) for(ll i=a;i<=ll(b);i++) #define FORD(i,a,b) for(ll i=a;i>=ll(b);i--) #define overload4(_1,_2,_3,_4,name,...) name #define overload3(_1,_2,_3,name,...) name #define rep1(n) for(ll i=0;i<n;++i) #define rep2(i,n) for(ll i=0;i<n;++i) #define rep3(i,a,b) for(ll i=a;i<b;++i) #define rep4(i,a,b,c) for(ll i=a;i<b;i+=c) #define rep(...) overload4(__VA_ARGS__,rep4,rep3,rep2,rep1)(__VA_ARGS__) #define rrep1(n) for(ll i=n;i--;) #define rrep2(i,n) for(ll i=n;i--;) #define rrep3(i,a,b) for(ll i=b;i-->(a);) #define rrep4(i,a,b,c) for(ll i=(a)+((b)-(a)-1)/(c)*(c);i>=(a);i-=c) #define rrep(...) overload4(__VA_ARGS__,rrep4,rrep3,rrep2,rrep1)(__VA_ARGS__) #define ALL(v) v.begin(), v.end() #define endl "\n" #define fi first #define se second #define popcount(bit) __builtin_popcount(bit) #define popcountll(bit) __builtin_popcountll(bit) #define pb push_back #define eb emplace_back #define sz(x) ((ll)(x).size()) using namespace atcoder; using P = pair<int, int>; using PL = pair<long long, long long>; using Graph = vector<vector<int>>; typedef long long ll; template<class T> inline bool chmax(T& a, T b) { if (a < b) { a = b; return 1; } return 0; } template<class T> inline bool chmin(T& a, T b) { if (a > b) { a = b; return 1; } return 0; } inline void Yes(bool b = true) { cout << (b ? "Yes" : "No") << '\n'; } inline void YES(bool b = true) { cout << (b ? "YES" : "NO") << '\n'; } inline void OKNG(bool b = true) { cout << (b ? "OK" : "NG") << '\n'; } const int dx[4] = {1, 0, -1, 0}; const int dy[4] = {0, 1, 0, -1}; const int fx[8] = {0, 1, 1, 1, 0, -1, -1, -1}; const int fy[8] = {1, 1, 0, -1, -1, -1, 0, 1}; template <typename T> const auto INF = numeric_limits<T>::max()/2; using mint = modint1000000007; string to_binary(ll N, int digit = 0) { string res = ""; while (N) { if (N & 1) res = "1" + res; else res = "0" + res; N >>= 1; } if (0 < digit) { int n = digit - res.length(); rep(i, 0, n) res = "0" + res; } return res; } int main() { ll n; cin >> n; vector<ll> s; while(n > 0){ s.push_back(n % 2); n /= 2; } reverse(ALL(s)); n = s.size(); mint com[n + 1][n + 2][2] = {0}; mint sum[n + 1][n + 2][2] = {0}; com[0][0][0] = 1; rep(i, n)rep(j, n+1)rep(k, 2){ ll nd = s[i]; rep(d, 2){ ll nk = k; ll nj = j; if(k == 0 && d > nd)continue; if(d < nd)nk = 1; if(d == 1)nj++; com[i+1][nj][nk] += com[i][j][k]; sum[i+1][nj][nk] += sum[i][j][k] * 2 + com[i][j][k] * d; } } mint ans = 0; rep(i, n+2)rep(j, 2)ans += sum[n][i][j] * i; cout << ans.val() << endl; }