#include <cstdio>
#include <stdio.h>
#include <iostream>
#include <math.h>
#include <vector>
#include <queue>
#include <stack>
#include <deque>
#include <algorithm>
#include <utility>
#include <set>
#include <map>
#include <stdlib.h>
#include <cstring>
#include <string.h>
#include <string>
#include <sstream>
#include <assert.h>
#include <climits>
#include <sstream>
#include <numeric>
#include <time.h>
#include <limits.h>
#include <list>
#include <bitset>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <iomanip>
#include <complex>
#include <chrono>
#include <fstream>
#include <functional>
#include <unistd.h>
// #pragma GCC optimize("Ofast,no-stack-protector")
// #pragma GCC optimize("O3,unroll-loops")
// #pragma GCC target("avx,avx2,bmi,bmi2,lzcnt,popcnt")
#define lowbit(x) ((x) & -(x))
#define ml(a, b) ((1ll * (a) * (b)) % M)
#define tml(a, b) (a) = ((1ll * (a) * (b)) % M)
#define ad(a, b) ((0ll + (a) + (b)) % M)
#define tad(a, b) (a) = ((0ll + (a) + (b)) % M)
#define mi(a, b) ((0ll + M + (a) - (b)) % M)
#define tmi(a, b) (a) = ((0ll + M + (a) - (b)) % M)
#define tmin(a, b) (a) = min((a), (b))
#define tmax(a, b) (a) = max((a), (b))
#define iter(a) (a).begin(), (a).end()
#define riter(a) (a).rbegin(), (a).rend()
#define initi(a, b) memset((a), (b), sizeof(a))
#define cpy(a, b) memcpy((a), (b), sizeof(a))
#define uni(a) a.resize(unique(iter(a)) - a.begin())
#define pb emplace_back
#define mpr make_pair
#define ls(i) ((i) << 1)
#define rs(i) ((i) << 1 | 1)
#define INF 0x3f3f3f3f
#define NIF 0xc0c0c0c0
#define eps 1e-9
#define F first
#define S second
#define AC cin.tie(0)->sync_with_stdio(0)
using namespace std;
typedef long long llt;
typedef pair<int, int> pii;
typedef pair<double, double> pdd;
typedef pair<llt, llt> pll;
typedef complex<double> cd;
// const int M = 998244353;
// random_device rm;
// mt19937 rg(rm());
// default_random_engine rg(rm());
// uniform_int_distribution<int> rd(INT_MIN, INT_MAX);
// uniform_real_distribution<double> rd(0, M_PI);
void db() { cerr << "\n"; }
template <class T, class... U>
void db(T a, U... b) { cerr << a << " ", db(b...); }
inline char gc()
{
const static int SZ = 1 << 16;
static char buf[SZ], *p1, *p2;
if (p1 == p2 && (p2 = buf + fread(p1 = buf, 1, SZ, stdin), p1 == p2))
return -1;
return *p1++;
}
void rd() {}
template <typename T, typename... U>
void rd(T &x, U &...y)
{
x = 0;
bool f = 0;
char c = gc();
while (!isdigit(c))
f ^= !(c ^ 45), c = gc();
while (isdigit(c))
x = (x << 1) + (x << 3) + (c ^ 48), c = gc();
f && (x = -x), rd(y...);
}
template <typename T>
void prt(T x)
{
if (x < 0)
putchar('-'), x = -x;
if (x > 9)
prt(x / 10);
putchar((x % 10) ^ 48);
}
const int N = 4e5 + 5;
int pt[N], dep[N], dst[N], n = 2, m = 2; // head to head
pii ac[20][N]; // po num
inline pii gid(int x)
{
pii res(upper_bound(pt + 1, pt + n + 1, x) - pt - 1, x);
res.S -= pt[res.F];
return res;
}
inline int gdst(const pii &p1, const pii &p2) // up down
{
return dst[p2.F] - dst[p1.F] + p2.S - p1.S;
}
inline pii lca(pii p1, pii p2)
{
if (dep[p1.F] < dep[p2.F])
swap(p1, p2);
for (int i = 19; i >= 0; i--)
if (ac[i][p1.F].F && dep[ac[i][p1.F].F] >= dep[p2.F])
p1 = ac[i][p1.F];
if (p1.F != p2.F)
{
for (int i = 19; i >= 0; i--)
if (ac[i][p1.F].F && ac[i][p2.F].F && ac[i][p1.F].F != ac[i][p2.F].F)
p1 = ac[i][p1.F], p2 = ac[i][p2.F];
p1 = ac[0][p1.F], p2 = ac[0][p2.F];
}
return pii(p1.F, min(p1.S, p2.S));
}
inline int qy(pii p, int d)
{
if (p.S >= d)
return pt[p.F] + p.S - d;
pii cr = p;
for (int i = 19; i >= 0; i--)
if (ac[i][cr.F].F && gdst(mpr(ac[i][cr.F].F, 0), p) < d)
cr = ac[i][cr.F];
d -= gdst(mpr(cr.F, 0), p);
cr = ac[0][cr.F];
return pt[cr.F] + cr.S - d + 1;
}
void init() { pt[1] = 1, pt[2] = 2; }
void path(int x, int y)
{
pii pr = gid(x);
dep[n] = dep[pr.F] + 1;
dst[n] = dst[pr.F] + pr.S + 1;
ac[0][n] = pr;
for (int i = 1; i < 20; i++)
ac[i][n] = ac[i - 1][ac[i - 1][n].F];
m += y, pt[++n] = m;
}
int dig(int x, int y)
{
pii p1 = gid(x), p2 = gid(y);
pii p = lca(p1, p2);
int d1 = gdst(p, p1), d2 = gdst(p, p2);
// db(d1, d2);
if (d1 >= d2)
return qy(p1, (d1 + d2) >> 1);
else
return qy(p2, (d1 + d2 + 1) >> 1);
}
// signed main()
// {
// init();
// int q, x, y;
// char c[10];
// scanf("%d", &q);
// while (q-- && scanf("%s%d%d", c, &x, &y))
// {
// if (c[0] == 'P')
// path(x, y);
// else
// printf("%d\n", dig(x, y)), fflush(stdout);
// }
// }
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
14 ms |
41304 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
9 ms |
41520 KB |
Output is correct |
| 2 |
Correct |
7 ms |
41564 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
296 ms |
67408 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
319 ms |
51656 KB |
Output is correct |
| 2 |
Correct |
245 ms |
67704 KB |
Output is correct |
| 3 |
Correct |
179 ms |
67616 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
313 ms |
67456 KB |
Output is correct |
| 2 |
Correct |
235 ms |
47208 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
472 ms |
67536 KB |
Output is correct |
| 2 |
Correct |
313 ms |
67724 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
129 ms |
51520 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
277 ms |
67800 KB |
Output is correct |
| 2 |
Correct |
161 ms |
79344 KB |
Output is correct |
| 3 |
Correct |
147 ms |
79444 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
296 ms |
67696 KB |
Output is correct |
| # |
Verdict |
Execution time |
Memory |
Grader output |
| 1 |
Correct |
303 ms |
67664 KB |
Output is correct |
