Submission #795267

#	Time	Username	Problem	Language	Result	Execution time	Memory
795267		Josia	Horses (IOI15_horses)	C++17	0 / 100	1563 ms	51348 KiB

horses

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

using namespace std;

#define ll long long

const ll mod = 1e9+7;

struct seg {
	vector<int> tree;

	void init(int n) {
		tree.assign(n*4, 0);
	}

	void update(int v, int rl, int rr, int pos, int x) {
		if (rl == rr) {
			tree[v] = x;
			return;
		}

		int rm = (rl + rr)/2;

		if (pos <= rm) update(v*2, rl, rm, pos, x);
		else update(v*2+1, rm+1, rr, pos, x);

		tree[v] = max(tree[v*2], tree[v*2+1]);
	}

	int query(int v, int rl, int rr, int ql, int qr) {
		if (ql > qr) return 0;

		if (rl == ql && rr == qr) return tree[v];

		int rm = (rl + rr)/2;

		return max(query(v*2, rl, rm, ql, min(qr, rm)), query(v*2+1, rm+1, rr, max(ql, rm+1), qr));
	}
};


struct segY {
	vector<pair<int, int>> tree;

	void init(int n) {
		tree.assign(n*4, {0, 0});
	}

	void update(int v, int rl, int rr, int pos, int x) {
		if (rl == rr) {
			tree[v] = {x, pos};
			return;
		}

		int rm = (rl + rr)/2;

		if (pos <= rm) update(v*2, rl, rm, pos, x);
		else update(v*2+1, rm+1, rr, pos, x);

		tree[v] = max(tree[v*2], tree[v*2+1]);
	}

	pair<int, int> query(int v, int rl, int rr, int ql, int qr) {
		if (ql > qr) return {0, -1};

		if (rl == ql && rr == qr) return tree[v];

		int rm = (rl + rr)/2;

		return max(query(v*2, rl, rm, ql, min(qr, rm)), query(v*2+1, rm+1, rr, max(ql, rm+1), qr));
	}
};


ll mult(ll a, ll b) {
	return (a*b)%mod;
}

struct seg2 {
	vector<ll> tree;

	void init(int n) {
		tree.assign(n*4, 0);
	}

	void update(int v, int rl, int rr, int pos, int x) {
		if (rl == rr) {
			tree[v] = x;
			return;
		}

		int rm = (rl + rr)/2;

		if (pos <= rm) update(v*2, rl, rm, pos, x);
		else update(v*2+1, rm+1, rr, pos, x);

		tree[v] = mult(tree[v*2], tree[v*2+1]);
	}

	ll query(int v, int rl, int rr, int ql, int qr) {
		if (ql > qr) return 1;

		if (rl == ql && rr == qr) return tree[v];

		int rm = (rl + rr)/2;

		return mult(query(v*2, rl, rm, ql, min(qr, rm)), query(v*2+1, rm+1, rr, max(ql, rm+1), qr));
	}
};


seg tree;
segY treeY;
seg2 tree2;
int N;
vector<ll> X, Y;

struct getSol {
	int sol = 0;

	int BS(int pos) {
		int l = 0, r = pos-1;
		while (l<r) {
			int m = (l+r+1)/2;
			// cerr << m << " " << pos-1 << " " << tree.query(1, 0, N-1, m, pos-1) << "\n";
			if (tree.query(1, 0, N-1, m, pos-1) > 1) {
				l = m;
			} else {
				r = m-1;
			}
		}

		if (tree.query(1, 0, N-1, l, pos-1) <= 1) return -1;

		return l;
	}

	int BSinv(int pos) {
		if (pos == N-1) return N-1;
		int l = pos+1, r = N-1;
		while (l<r) {
			int m = (l+r+1)/2;
			// cerr << m << " " << pos-1 << " " << tree.query(1, 0, N-1, m, pos-1) << "\n";
			if (tree.query(1, 0, N-1, pos+1, m) > 1) {
				r = m-1;
			} else {
				l = m;
			}
		}

		return l;
	}


	int findInd() {
		int currInd = N-1;
		ll currProd = 1;
		for (int i = N-1; i>=0; i = BS(i)) {
			// cerr << i << " " << currProd << " " << currInd << " " << BSinv(i) << "\n";
			if (currProd*Y[currInd] <= treeY.query(1, 0, N-1, i, BSinv(i)).first) {currInd = treeY.query(1, 0, N-1, i, BSinv(i)).second; currProd=1;}


			currProd *= X[i];
			if (currProd > 1e9) break;
		}
		return currInd;
	}


	getSol() {
		int ind = findInd();
		// cerr << ind << " " << tree2.query(1, 0, N-1, 0, ind) << " " << Y[ind] << "\n\n";

		sol = (tree2.query(1, 0, N-1, 0, ind) * Y[ind])%mod;
	}
};




int init(int _N, int _X[], int _Y[]) {
	N = _N;
	X.resize(N);
	Y.resize(N);
	for (int i = 0; i<N; i++) X[i] = _X[i];
	for (int i = 0; i<N; i++) Y[i] = _Y[i];

	tree.init(N);
	treeY.init(N);
	tree2.init(N);

	for (int i = 0; i<N; i++) {
		tree.update(1, 0, N-1, i, X[i]);
		treeY.update(1, 0, N-1, i, Y[i]);
		tree2.update(1, 0, N-1, i, X[i]);
	}
	return getSol().sol;
}

int updateX(int pos, int val) {
	tree.update(1, 0, N-1, pos, val);
	tree2.update(1, 0, N-1, pos, val);
	X[pos] = val;

	return getSol().sol;
}

int updateY(int pos, int val) {
	treeY.update(1, 0, N-1, pos, val);
	Y[pos] = val;

	return getSol().sol;
}

Compilation message (stderr)

horses.cpp: In member function 'int getSol::findInd()':
horses.cpp:165:8: warning: conversion from 'long long int' to 'double' may change value [-Wconversion]
  165 |    if (currProd > 1e9) break;
      |        ^~~~~~~~
horses.cpp: In constructor 'getSol::getSol()':
horses.cpp:175:50: warning: conversion from 'long long int' to 'int' may change value [-Wconversion]
  175 |   sol = (tree2.query(1, 0, N-1, 0, ind) * Y[ind])%mod;
      |         ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~^~~~
horses.cpp: In function 'int init(int, int*, int*)':
horses.cpp:194:33: warning: conversion from '__gnu_cxx::__alloc_traits<std::allocator<long long int>, long long int>::value_type' {aka 'long long int'} to 'int' may change value [-Wconversion]
  194 |   tree.update(1, 0, N-1, i, X[i]);
      |                                 ^
horses.cpp:195:34: warning: conversion from '__gnu_cxx::__alloc_traits<std::allocator<long long int>, long long int>::value_type' {aka 'long long int'} to 'int' may change value [-Wconversion]
  195 |   treeY.update(1, 0, N-1, i, Y[i]);
      |                                  ^
horses.cpp:196:34: warning: conversion from '__gnu_cxx::__alloc_traits<std::allocator<long long int>, long long int>::value_type' {aka 'long long int'} to 'int' may change value [-Wconversion]
  196 |   tree2.update(1, 0, N-1, i, X[i]);
      |                                  ^

• Subtask 1: 0 / 17
• Subtask 2: 0 / 17
• Subtask 3: 0 / 20
• Subtask 4: 0 / 23
• Subtask 5: 0 / 23