#include <bits/stdc++.h>
using ll = long long;
using namespace std;
#define forn(a, b) for (int a = 0; a < (b); a++)
#define rofn(a, b) for (int a = (b) - 1; a >= 0; a--)
#define LOG(a) (63 - __builtin_clzll(a))
#define gc() (pt1==pt2&&(pt2=(pt1=buf)+fread(buf,1,1<<21,stdin),pt1==pt2)?EOF:*pt1++)
char buf[1<<21],*pt1,*pt2;
inline ll read() {
ll x = 0;char c;
while(!isdigit(c = gc()));
while(isdigit(c)) x = (x<<3)+(x<<1)+(c^48), c = gc();
return x;
}
const ll INFL = 1ll << 60;
const int INF = 1 << 29;
const int N = 61;
const int M = 2 * N;
const int MOD = 1e9 + 7;
const int X = 2.6e8;
ll d0[M][M], d1[M][M];
ll n;
int ptr, ans, q;
int pw2[X + 1];
//unordered_map<ll, int> calced;
inline ll pw(ll a, int b) {
ll re = 1;
while (b) {
if (b & 1)
re = (re * a) % MOD;
a = (a * a) % MOD;
b >>= 1;
}
return re;
}
inline int p2(ll b) {
b %= (MOD - 1);
ll a = 1;
while (b > X) {
a = (a * pw2[X]) % MOD;
b -= X;
}
return a * pw2[b] % MOD;
}
void calc_dist(ll v) {
for (int k = LOG(v) + 1; k < N; k++) {
int x = k - LOG(v) - 1;
if ((v << x) > n)
d1[ptr][x] = 0;
else
d1[ptr][x] = min(n - (v << x) + 1, ((ll)1 << x));
}
if (v > 1) {
for (int k = 2; k < N; k++)
d0[ptr][k] -= d1[ptr][k - 2];
}
}
inline void pass() {
for (int k = 1; k < M; k++)
d0[ptr + 1][k] = d0[ptr][k - 1] + d1[ptr][k - 1];
ptr++;
}
void Complete(ll v) {
//cout << "C: " << v << ' ' << sub << " PTR: " << ptr << '\n';
calc_dist(v);
if ((v << 1) <= n) {
pass();
Complete(v << 1);
q <<= 1;
if (q >= MOD)
q -= MOD;
}
ll s = 0;
for (int k = 0; s != n; k++) {
s += d0[ptr][k] + d1[ptr][k];
q += p2(n - s) - 1;
if (q >= MOD)
q -= MOD;
}
ptr--;
}
inline int solve(ll a, ll b) {
int x = LOG(b) - LOG(a);
a <<= x;
return x + (a <= b);
}
void F(ll v, ll sub) {
if (((sub) & (sub + 1)) == 0) {
q = 0;
Complete(v);
ans += q;
if (ans >= MOD)
ans -= MOD;
return;
}
//cout << "F: " << v << ' ' << sub << " PTR: " << ptr << '\n';
calc_dist(v);
ll r = (v << 1) | 1;
ll l = (v << 1);
int cntr = solve(r, n), cntl = solve(l, n);
if (cntr == cntl) {
q = 0;
if ((v << 1) <= n) {
pass();
Complete(v << 1);
ans += q;
if (ans >= MOD)
ans -= MOD;
}
if (((v << 1) | 1) <= n) {
pass();
F((v << 1) | 1, sub - (1ll << cntl));
}
}
else {
if (((v << 1) | 1) <= n) {
q = 0;
pass();
Complete((v << 1) | 1);
ans += q;
if (ans >= MOD)
ans -= MOD;
}
if ((v << 1) <= n) {
pass();
F((v << 1), sub - (1ll << cntr));
}
}
ll s = 0;
for (int k = 0; s != n; k++) {
s += d0[ptr][k] + d1[ptr][k];
ans += p2(n - s) - 1;
if (ans >= MOD)
ans -= MOD;
}
ptr--;
}
signed main() {
ios::sync_with_stdio(0), cin.tie(0);
pw2[0] = 1;
for (int i = 1; i <= X; i++) {
pw2[i] = pw2[i - 1] * 2;
if (pw2[i] >= MOD)
pw2[i] -= MOD;
}
int query = read();
//cin >> query;
while (query--) {
//cin >> n;
n = read();
F(1, n);
ll re = p2(n) - 1;
ans = (ans * pw(re, MOD - 2)) % MOD;
cout << ans << '\n';
ans = q = ptr = 0;
}
return 0;
}
//a
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1280 ms |
1017880 KB |
Output is correct |
2 |
Correct |
1287 ms |
1017840 KB |
Output is correct |
3 |
Correct |
1274 ms |
1017864 KB |
Output is correct |
4 |
Correct |
1287 ms |
1017872 KB |
Output is correct |
5 |
Correct |
1278 ms |
1017984 KB |
Output is correct |
6 |
Correct |
1289 ms |
1018088 KB |
Output is correct |
7 |
Correct |
1255 ms |
1017900 KB |
Output is correct |
8 |
Correct |
1266 ms |
1017884 KB |
Output is correct |
9 |
Correct |
1278 ms |
1017988 KB |
Output is correct |
10 |
Correct |
1321 ms |
1018000 KB |
Output is correct |
11 |
Correct |
1306 ms |
1018032 KB |
Output is correct |
12 |
Correct |
1271 ms |
1017860 KB |
Output is correct |
13 |
Correct |
1676 ms |
1017932 KB |
Output is correct |
14 |
Correct |
1297 ms |
1017928 KB |
Output is correct |
15 |
Correct |
1279 ms |
1017940 KB |
Output is correct |
16 |
Correct |
1286 ms |
1017960 KB |
Output is correct |
17 |
Correct |
1302 ms |
1017940 KB |
Output is correct |
18 |
Correct |
1299 ms |
1017948 KB |
Output is correct |
19 |
Correct |
1294 ms |
1018060 KB |
Output is correct |
20 |
Correct |
1319 ms |
1017976 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1316 ms |
1018044 KB |
Output is correct |
2 |
Correct |
1322 ms |
1017984 KB |
Output is correct |
3 |
Correct |
1306 ms |
1018024 KB |
Output is correct |
4 |
Correct |
1323 ms |
1018160 KB |
Output is correct |
5 |
Correct |
1306 ms |
1017952 KB |
Output is correct |
6 |
Correct |
1316 ms |
1018056 KB |
Output is correct |
7 |
Correct |
1327 ms |
1018044 KB |
Output is correct |
8 |
Correct |
1316 ms |
1017964 KB |
Output is correct |
9 |
Correct |
1305 ms |
1018124 KB |
Output is correct |
10 |
Correct |
1358 ms |
1017972 KB |
Output is correct |
11 |
Correct |
1445 ms |
1018112 KB |
Output is correct |
12 |
Correct |
1442 ms |
1018096 KB |
Output is correct |
13 |
Correct |
1424 ms |
1018084 KB |
Output is correct |
14 |
Correct |
1441 ms |
1018060 KB |
Output is correct |
15 |
Correct |
1445 ms |
1018076 KB |
Output is correct |
16 |
Correct |
1864 ms |
1017988 KB |
Output is correct |
17 |
Correct |
1434 ms |
1018068 KB |
Output is correct |
18 |
Correct |
1419 ms |
1018304 KB |
Output is correct |
19 |
Correct |
1472 ms |
1018004 KB |
Output is correct |
20 |
Correct |
1415 ms |
1018084 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1383 ms |
1018096 KB |
Output is correct |
2 |
Correct |
1415 ms |
1018048 KB |
Output is correct |
3 |
Correct |
1404 ms |
1018052 KB |
Output is correct |
4 |
Correct |
1420 ms |
1017976 KB |
Output is correct |
5 |
Correct |
1402 ms |
1018044 KB |
Output is correct |
6 |
Correct |
1404 ms |
1017936 KB |
Output is correct |
7 |
Correct |
1416 ms |
1018112 KB |
Output is correct |
8 |
Correct |
1424 ms |
1018052 KB |
Output is correct |
9 |
Correct |
1301 ms |
1017936 KB |
Output is correct |
10 |
Correct |
1295 ms |
1018080 KB |
Output is correct |
11 |
Correct |
1319 ms |
1017952 KB |
Output is correct |
12 |
Correct |
1316 ms |
1018048 KB |
Output is correct |
13 |
Correct |
1691 ms |
1018120 KB |
Output is correct |
14 |
Correct |
1674 ms |
1018148 KB |
Output is correct |
15 |
Correct |
1820 ms |
1017976 KB |
Output is correct |
16 |
Correct |
1771 ms |
1018120 KB |
Output is correct |
17 |
Correct |
1670 ms |
1017928 KB |
Output is correct |
18 |
Correct |
1471 ms |
1017924 KB |
Output is correct |
19 |
Correct |
1428 ms |
1018004 KB |
Output is correct |
20 |
Correct |
1454 ms |
1017940 KB |
Output is correct |
# |
결과 |
실행 시간 |
메모리 |
Grader output |
1 |
Correct |
1280 ms |
1017880 KB |
Output is correct |
2 |
Correct |
1287 ms |
1017840 KB |
Output is correct |
3 |
Correct |
1274 ms |
1017864 KB |
Output is correct |
4 |
Correct |
1287 ms |
1017872 KB |
Output is correct |
5 |
Correct |
1278 ms |
1017984 KB |
Output is correct |
6 |
Correct |
1289 ms |
1018088 KB |
Output is correct |
7 |
Correct |
1255 ms |
1017900 KB |
Output is correct |
8 |
Correct |
1266 ms |
1017884 KB |
Output is correct |
9 |
Correct |
1278 ms |
1017988 KB |
Output is correct |
10 |
Correct |
1321 ms |
1018000 KB |
Output is correct |
11 |
Correct |
1306 ms |
1018032 KB |
Output is correct |
12 |
Correct |
1271 ms |
1017860 KB |
Output is correct |
13 |
Correct |
1676 ms |
1017932 KB |
Output is correct |
14 |
Correct |
1297 ms |
1017928 KB |
Output is correct |
15 |
Correct |
1279 ms |
1017940 KB |
Output is correct |
16 |
Correct |
1286 ms |
1017960 KB |
Output is correct |
17 |
Correct |
1302 ms |
1017940 KB |
Output is correct |
18 |
Correct |
1299 ms |
1017948 KB |
Output is correct |
19 |
Correct |
1294 ms |
1018060 KB |
Output is correct |
20 |
Correct |
1319 ms |
1017976 KB |
Output is correct |
21 |
Correct |
1316 ms |
1018044 KB |
Output is correct |
22 |
Correct |
1322 ms |
1017984 KB |
Output is correct |
23 |
Correct |
1306 ms |
1018024 KB |
Output is correct |
24 |
Correct |
1323 ms |
1018160 KB |
Output is correct |
25 |
Correct |
1306 ms |
1017952 KB |
Output is correct |
26 |
Correct |
1316 ms |
1018056 KB |
Output is correct |
27 |
Correct |
1327 ms |
1018044 KB |
Output is correct |
28 |
Correct |
1316 ms |
1017964 KB |
Output is correct |
29 |
Correct |
1305 ms |
1018124 KB |
Output is correct |
30 |
Correct |
1358 ms |
1017972 KB |
Output is correct |
31 |
Correct |
1445 ms |
1018112 KB |
Output is correct |
32 |
Correct |
1442 ms |
1018096 KB |
Output is correct |
33 |
Correct |
1424 ms |
1018084 KB |
Output is correct |
34 |
Correct |
1441 ms |
1018060 KB |
Output is correct |
35 |
Correct |
1445 ms |
1018076 KB |
Output is correct |
36 |
Correct |
1864 ms |
1017988 KB |
Output is correct |
37 |
Correct |
1434 ms |
1018068 KB |
Output is correct |
38 |
Correct |
1419 ms |
1018304 KB |
Output is correct |
39 |
Correct |
1472 ms |
1018004 KB |
Output is correct |
40 |
Correct |
1415 ms |
1018084 KB |
Output is correct |
41 |
Correct |
1383 ms |
1018096 KB |
Output is correct |
42 |
Correct |
1415 ms |
1018048 KB |
Output is correct |
43 |
Correct |
1404 ms |
1018052 KB |
Output is correct |
44 |
Correct |
1420 ms |
1017976 KB |
Output is correct |
45 |
Correct |
1402 ms |
1018044 KB |
Output is correct |
46 |
Correct |
1404 ms |
1017936 KB |
Output is correct |
47 |
Correct |
1416 ms |
1018112 KB |
Output is correct |
48 |
Correct |
1424 ms |
1018052 KB |
Output is correct |
49 |
Correct |
1301 ms |
1017936 KB |
Output is correct |
50 |
Correct |
1295 ms |
1018080 KB |
Output is correct |
51 |
Correct |
1319 ms |
1017952 KB |
Output is correct |
52 |
Correct |
1316 ms |
1018048 KB |
Output is correct |
53 |
Correct |
1691 ms |
1018120 KB |
Output is correct |
54 |
Correct |
1674 ms |
1018148 KB |
Output is correct |
55 |
Correct |
1820 ms |
1017976 KB |
Output is correct |
56 |
Correct |
1771 ms |
1018120 KB |
Output is correct |
57 |
Correct |
1670 ms |
1017928 KB |
Output is correct |
58 |
Correct |
1471 ms |
1017924 KB |
Output is correct |
59 |
Correct |
1428 ms |
1018004 KB |
Output is correct |
60 |
Correct |
1454 ms |
1017940 KB |
Output is correct |
61 |
Correct |
1445 ms |
1017820 KB |
Output is correct |
62 |
Correct |
1465 ms |
1017956 KB |
Output is correct |
63 |
Correct |
2507 ms |
1018064 KB |
Output is correct |
64 |
Correct |
2458 ms |
1017992 KB |
Output is correct |
65 |
Correct |
2417 ms |
1018060 KB |
Output is correct |
66 |
Correct |
2470 ms |
1018068 KB |
Output is correct |
67 |
Correct |
2497 ms |
1018076 KB |
Output is correct |
68 |
Correct |
2545 ms |
1018072 KB |
Output is correct |
69 |
Correct |
2329 ms |
1018252 KB |
Output is correct |
70 |
Correct |
2311 ms |
1018076 KB |
Output is correct |
71 |
Correct |
1333 ms |
1018016 KB |
Output is correct |
72 |
Correct |
1316 ms |
1018308 KB |
Output is correct |
73 |
Correct |
1412 ms |
1018024 KB |
Output is correct |
74 |
Correct |
1689 ms |
1018144 KB |
Output is correct |
75 |
Execution timed out |
3088 ms |
1018088 KB |
Time limit exceeded |
76 |
Halted |
0 ms |
0 KB |
- |