Submission #798895

#	Time	Username	Problem	Language	Result	Execution time	Memory
798895		QwertyPi	Game (IOI13_game)	C++14	63 / 100	1275 ms	256000 KiB

game

#include "game.h"
#include <bits/stdc++.h>

using namespace std;

long long gcd2(long long X, long long Y) {
	return __gcd(X, Y);
}

struct node{
	long long a = 0;
	node *px1 = 0, *px2 = 0, *py1 = 0, *py2 = 0;
};

int R, C;

long long get(node *v){
	if(!v) return 0;
	return v->a; 
}

node root;
void update_y(int y, long long K, node *v, int y1 = 0, int y2 = C - 1){
	if(y1 == y2){
		v->a = K;
	}else{
		int ym = (y1 + y2) / 2;
		if(y <= ym){
			if(!v->py1) v->py1 = new node();
			update_y(y, K, v->py1, y1, ym);
		}else{
			if(!v->py2) v->py2 = new node();
			update_y(y, K, v->py2, ym + 1, y2);
		}
		v->a = gcd2(get(v->py1), get(v->py2));
	}
}

long long calc_y(int qy1, int qy2, node *v, int y1 = 0, int y2 = C - 1){
	if(!v || qy1 > y2 || y1 > qy2) return 0;
	if(qy1 <= y1 && y2 <= qy2) return v->a;
	int ym = (y1 + y2) / 2;
	long long a1 = 0, a2 = 0;
	if(v->py1) a1 = calc_y(qy1, qy2, v->py1, y1, ym);
	if(v->py2) a2 = calc_y(qy1, qy2, v->py2, ym + 1, y2);
	return gcd2(a1, a2);
}

void update_x(int x, int y, long long K, node *v, int x1 = 0, int x2 = R - 1){
	if(x1 == x2){
		update_y(y, K, v);
		return;
	}
	int xm = (x1 + x2) / 2;
	if(x <= xm){
		if(!v->px1) v->px1 = new node();
		update_x(x, y, K, v->px1, x1, xm);
	}else{
		if(!v->px2) v->px2 = new node();
		update_x(x, y, K, v->px2, xm + 1, x2);
	}
	long long A = gcd2(calc_y(y, y, v->px1), calc_y(y, y, v->px2));
	update_y(y, A, v);
}

long long calc_x(int qx1, int qx2, int qy1, int qy2, node *v, int x1 = 0, int x2 = R - 1){
	if(!v || qx1 > x2 || x1 > qx2) return 0;
	if(qx1 <= x1 && x2 <= qx2) return calc_y(qy1, qy2, v);
	int xm = (x1 + x2) / 2;
	long long a1 = 0, a2 = 0;
	if(v->px1) a1 = calc_x(qx1, qx2, qy1, qy2, v->px1, x1, xm);
	if(v->px2) a2 = calc_x(qx1, qx2, qy1, qy2, v->px2, xm + 1, x2);
	return gcd2(a1, a2);
}

void init(int R, int C) {
	root = node(); ::R = R, ::C = C;
}

void update(int P, int Q, long long K) {
    update_x(P, Q, K, &root);
}

long long calculate(int P, int Q, int U, int V) {
    return calc_x(P, U, Q, V, &root);
}

• Subtask 1: 10 / 10
• Subtask 2: 27 / 27
• Subtask 3: 26 / 26
• Subtask 4: 0 / 17
• Subtask 5: 0 / 20