| # | 제출 시각 | 아이디 | 문제 | 언어 | 결과 | 실행 시간 | 메모리 |
|---|---|---|---|---|---|---|---|
| 140322 | edenooo | Game (IOI13_game) | C++17 | 0 ms | 0 KiB |
이 제출은 이전 버전의 oj.uz에서 채점하였습니다. 현재는 제출 당시와는 다른 서버에서 채점을 하기 때문에, 다시 제출하면 결과가 달라질 수도 있습니다.
#include<stdio.h>
#include<iostream>
#include<vector>
#include<string>
#include<queue>
#include<deque>
#include<set>
#include<map>
#include<algorithm>
#include<math.h>
//#include "game.h"
using namespace std;
#define INF 1234567890
#define ll long long
struct SEG {
int l, r, root;
};
struct Node {
int l, r, s, e;
ll val;
};
int Y, X, Q;
vector<SEG> segy(1);
vector<Node> seg(1);
ll Gcd(ll a, ll b)
{
return a ? Gcd(b%a, a) : b;
}
ll QueryX(int L, int R, int n)
{
int l = seg[n].s, r = seg[n].e;
if (r < L || R < l) return 0LL;
if (L <= l && r <= R) return seg[n].val;
ll left = 0, right = 0;
if (seg[n].l)
left = QueryX(L, R, seg[n].l);
if (seg[n].r)
right = QueryX(L, R, seg[n].r);
return Gcd(left, right);
}
ll QueryY(int sy, int sx, int ey, int ex, int n, int l, int r)
{
if (r < sy || ey < l) return 0LL;
if (sy <= l && r <= ey) return QueryX(sx, ex, segy[n].root);
int mid = (int)((ll)l + r >> 1);
ll left = 0, right = 0;
if (segy[n].l)
left = QueryY(sy, sx, ey, ex, segy[n].l, l, mid);
if (segy[n].r)
right = QueryY(sy, sx, ey, ex, segy[n].r, mid + 1, r);
return Gcd(left, right);
}
void UpdateX(int idx, ll val, int n, int l, int r)
{
if (r < idx || idx < l) return;
if (l == r)
{
seg[n].val = val;
return;
}
int mid = (int)((ll)l + r >> 1);
if (idx <= mid)
{
if (seg[n].l == 0)
{
seg[n].l = seg.size();
seg.push_back({ 0,0,idx,idx,val });
}
else if (seg[seg[n].l].s <= idx && idx <= seg[seg[n].l].e)
UpdateX(idx, val, seg[n].l, seg[seg[n].l].s, seg[seg[n].l].e);
else // 새로운 LCA 만들고 부모 자식 관계 재설정하기
{
int tmp = seg[n].l, tmp2 = seg.size();
seg.push_back({ 0,0,idx,idx,val });
if (seg[tmp].s > seg[tmp2].s) swap(tmp, tmp2);
int s = 0, e = X - 1; //LCA의 구간 찾기
while (1)
{
int m = s + e >> 1;
if (s <= seg[tmp].s && seg[tmp2].e <= m)
e = m;
else if (m + 1 <= seg[tmp].s && seg[tmp2].e <= e)
s = m + 1;
else
break;
}
seg[n].l = seg.size();
seg.push_back({ tmp,tmp2,s,e,Gcd(seg[tmp].val,seg[tmp2].val) });
}
}
else
{
if (seg[n].r == 0)
{
seg[n].r = seg.size();
seg.push_back({ 0,0,idx,idx,val });
}
else if (seg[seg[n].r].s <= idx && idx <= seg[seg[n].r].e)
UpdateX(idx, val, seg[n].r, seg[seg[n].r].s, seg[seg[n].r].e);
else // 새로운 LCA 만들고 부모 자식 관계 재설정하기
{
int tmp = seg[n].r, tmp2 = seg.size();
seg.push_back({ 0,0,idx,idx,val });
if (seg[tmp].s > seg[tmp2].s) swap(tmp, tmp2);
int s = 0, e = X - 1; //LCA의 구간 찾기
while (1)
{
int m = s + e >> 1;
if (s <= seg[tmp].s && seg[tmp2].e <= m)
e = m;
else if (m + 1 <= seg[tmp].s && seg[tmp2].e <= e)
s = m + 1;
else
break;
}
seg[n].r = seg.size();
seg.push_back({ tmp,tmp2,s,e,Gcd(seg[tmp].val,seg[tmp2].val) });
}
}
seg[n].val = Gcd(seg[seg[n].l].val, seg[seg[n].r].val);
}
void UpdateY(int y, int x, ll val, int n, int l, int r)
{
if (r < y || y < l) return;
if (l == r)
{
//역원이 없으면 업데이트를 밖으로 빼면 안 됨!
UpdateX(x, val, segy[n].root, 0, X - 1);
return;
}
int mid = (int)((ll)l + r >> 1);
if (y <= mid)
{
if (segy[n].l == 0)
{
segy[n].l = segy.size();
segy.push_back({ 0,0,(int)seg.size() });
seg.push_back({ 0,0,0,X - 1,0 });
}
UpdateY(y, x, val, segy[n].l, l, mid);
}
else
{
if (segy[n].r == 0)
{
segy[n].r = segy.size();
segy.push_back({ 0,0,(int)seg.size() });
seg.push_back({ 0,0,0,X - 1,0 });
}
UpdateY(y, x, val, segy[n].r, mid + 1, r);
}
//고친 부분 (중요)
ll updL = 0, updR = 0;
if (segy[n].l)
updL = QueryX(x, x, segy[segy[n].l].root);
if (segy[n].r)
updR = QueryX(x, x, segy[segy[n].r].root);
UpdateX(x, Gcd(updL, updR), segy[n].root, 0, X - 1);
}
void init(int R, int C)
{
Y = R, X = C;
segy.push_back({ 0,0,(int)seg.size() });
seg.push_back({ 0,0,0,X-1,0 });
}
void update(int P, int Q, ll K)
{
UpdateY(P, Q, K, 1, 0, Y - 1);
}
ll calculate(int P, int Q, int U, int V)
{
return QueryY(P, Q, U, V, 1, 0, Y - 1);
}
int main()
{
scanf("%d %d %d", &Y, &X, &Q);
init(Y, X);
while (Q--)
{
int t, ay, ax, by, bx;
ll k;
scanf("%d %d %d", &t, &ay, &ax);
if (t == 1)
{
scanf("%lld", &k);
update(ay, ax, k);
}
else
{
scanf("%d %d", &by, &bx);
printf("%lld\n", calculate(ay, ax, by, bx));
}
}
return 0;
}