QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#883420 | #8593. Coin | arnur2937 | 0 | 0ms | 12016kb | C++23 | 6.1kb | 2025-02-05 16:17:42 | 2025-02-05 16:17:59 |
Judging History
answer
#include <bits/stdc++.h>
// #pragma GCC optimize("Ofast,unroll-loops,no-stack-protector,fast-math,inline")
// #pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,lzcnt,mmx,abm,avx,avx2,fma")
using namespace std;
#define static_assert(...);
#define speed ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#define bigInt __int128
#define int long long
#define dl double long
#define fl float
#define all(s) s.begin(), s.end()
#define pub push_back
#define puf push_front
#define pob pop_back
#define pof pob_front
#define ins insert
#define F first
#define S second
#define len length
const int N = 100005;
const int M = 1005;
const int LN = 131072;
const int MOD = 1e9 + 7;//998244353;
const int BLOCK = 500;
int binpow(int a, int b) {//, int MOD) {
int res = 1;
a %= MOD;
while (b > 0) {
if (b & 1) {
res = res * a;
res %= MOD;
}
a = a * a;
a %= MOD;
b >>= 1;
}
return res;
}
int mdiv(int a, int b) {
int ret = (a % MOD) * binpow(b, MOD - 2);
ret %= MOD;
return ret;
}
int mmul(int a, int b) {
int ret = (a % MOD) * (b % MOD);
ret %= MOD;
return ret;
}
int madd(int a, int b) {
int ret = (a % MOD) + (b % MOD);
ret %= MOD;
return ret;
}
int msub(int a, int b) {
int ret = (a % MOD) - (b % MOD);
ret = (ret + MOD) % MOD;
return ret;
}
int GCD(int a, int b) {
if (b == 0) {
return a;
}
return GCD(b, a % b);
}
struct pqComp
{
bool operator()(const pair<int, int>& p1, const pair<int, int>& p2) const
{
return (p1.F < p2.F) || (p1.F == p2.F && p1.S < p2.S);
}
};
bool pCompF(pair<int, int>& p1, pair<int, int>& p2)
{
return p1.F < p2.F;
}
bool pCompS(const pair<int, int>& p1, const pair<int, int>& p2)
{
return p1.S < p2.S;
}
bool pCompFS(pair<int, int>& p1, pair<int, int>& p2)
{
return (p1.S < p2.S) || (p1.S == p2.S && p1.F < p2.F);
}
vector <vector<int>> DS {{0, -1}, {0, 1}, {-1, 0}, {1, 0}};
int n, m, ln = 1, x[2*N], y[2*N], mn[2*N], mx[2*N], lf[2*N], rg[2*N], vis[2*N];
vector<int> tp, g[2*N];
pair<int,int> tmn[8*N], tmx[8*N];
void dfs(int v) {
vis[v] = 1;
for (int u: g[v]) {
if (vis[u]) continue;
dfs(u);
}
tp.pub(v);
}
void updmn(int pos, int x) {
int v = pos + ln - 1;
tmn[v] = {x, pos};
v /= 2;
while (v) {
tmn[v] = (tmn[v * 2].F < tmn[v * 2 + 1].F ? tmn[v * 2] : tmn[v * 2 + 1]);
v /= 2;
}
}
void updmx(int pos, int x) {
int v = pos + ln - 1;
tmx[v] = {x, pos};
v /= 2;
while (v) {
tmx[v] = (tmx[v * 2].F > tmx[v * 2 + 1].F ? tmx[v * 2] : tmx[v * 2 + 1]);
v /= 2;
}
}
pair<int,int> getmn(int l, int r, int v = 1, int tl = 1, int tr = ln) {
if (r < tl || tr < l) {
return {1e9, 0};
}
if (l <= tl && tr <= r) {
return tmn[v];
}
int m = (tl + tr) / 2;
auto xl = getmn(l, r, v * 2, tl, m), xr = getmn(l, r, v * 2 + 1, m + 1, tr);
return (xl.F < xr.F ? xl : xr);
}
pair<int,int> getmx(int l, int r, int v = 1, int tl = 1, int tr = ln) {
if (r < tl || tr < l) {
return {0, 0};
}
if (l <= tl && tr <= r) {
return tmx[v];
}
int m = (tl + tr) / 2;
auto xl = getmx(l, r, v * 2, tl, m), xr = getmx(l, r, v * 2 + 1, m + 1, tr);
return (xl.F > xr.F ? xl : xr);
}
int upd(int i) {
int ok = 0;
while (lf[i] >= 1 && mn[lf[i]] <= i) {
lf[i] = lf[lf[i]];
ok = 1;
}
while (rg[i] <= n && mx[rg[i]] >= i) {
rg[i] = rg[rg[i]];
ok = 1;
}
return ok;
}
void solve() {
cin >> n >> m;
while (ln < n) {
ln *= 2;
}
for (int i = 1; i < ln * 2; ++i) {
tmn[i].F = 1e9;
}
for (int i = 1; i <= m; ++i) {
cin >> x[i] >> y[i];
g[x[i]].pub(y[i]);
}
for (int i = 1; i <= n; ++i) {
if (!vis[i]) {
dfs(i);
}
}
tp.pub(0);
reverse(all(tp));
// for (int i = 1; i <= n; ++i) {
// cout << tp[i] << ' ';
// }
// cout << '\n';
vector<int> pos(n + 1);
for (int i = 1; i <= n; ++i) {
pos[tp[i]] = i;
mn[i] = 1e9;
lf[i] = i - 1;
rg[i] = i + 1;
updmn(i, lf[i]);
updmx(i, rg[i]);
}
vector<int> ans(n + 1, -1);
for (int i = 1; i <= m; ++i) {
mx[pos[y[i]]] = max(mx[pos[y[i]]], pos[x[i]]);
mn[pos[x[i]]] = min(mn[pos[x[i]]], pos[y[i]]);
vector<array<int, 3>> add;
while (true) {
auto j = getmn(pos[y[i]], n);
//if (i == 8)
//cout << "! " << j.S << ' ' << j.F << ' ' << lf[j.S] << ' ' << rg[j.S] << '\n';
if (j.F >= pos[y[i]] || j.F <= pos[x[i]]) break;
upd(j.S);
add.pub({0, j.S, lf[j.S]});
updmn(j.S, 1e9);
}
//return;
while (true) {
auto j = getmx(1, pos[x[i]]);
if (j.F <= pos[x[i]] || j.F >= pos[y[i]]) break;
upd(j.S);
add.pub({1, j.S, rg[j.S]});
updmx(j.S, 0);
}
// for (int j = 1; j <= n; ++j) {
// cout << lf[j] << ' ' << rg[j] << '\n';
// }
// cout << '\n';
for (auto [t, x, y]: add) {
//cout << t << ' ' << x << ' ' << y << '\n';
if (t) {
if (y < n + 1) {
updmx(x, y);
}
}
else {
if (y > 0) {
updmn(x, y);
}
}
if (lf[x] == 0 && rg[x] == n + 1) {
ans[tp[x]] = i;
}
}
}
for (int i = 1; i <= n; ++i) {
cout << ans[i] << ' ';
}
}
signed main() {
speed;
int T = 1;
//cin >> T;
while (T--) {
solve();
}
}
/*
НЕ ЗАХОДИТ РЕШЕНИЕ?
1) ПРОВЕРЬ НА ОЧЕВИДНЫЕ ОШИБКИ В КОДЕ
2) ПРОВЕРЬ НА ПЕРЕПОЛНЕНИЯ
3) УБЕДИСЬ В ПРАВИЛЬНОСТИ АЛГОРИТМА
*/
详细
Subtask #1:
score: 0
Wrong Answer
Test #1:
score: 0
Wrong Answer
time: 0ms
memory: 12016kb
input:
4 4 2 4 3 1 4 1 2 3
output:
-1 -1 -1 -1
result:
wrong answer wa
Subtask #2:
score: 0
Skipped
Dependency #1:
0%
Subtask #3:
score: 0
Skipped
Dependency #1:
0%
Subtask #4:
score: 0
Skipped
Dependency #1:
0%