You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
270 lines
6.2 KiB
270 lines
6.2 KiB
// https://www.hackerrank.com/challenges/arithmetic-progressions
|
|
|
|
#include <stdio.h>
|
|
#include <stdlib.h>
|
|
#include <string.h>
|
|
#include <ctype.h>
|
|
|
|
#define INTS_START_SIZE 8
|
|
#define INTS_DELIM " \t\n"
|
|
#define BUF_SIZE 1 << 20
|
|
#define LINE_SIZE 16
|
|
#define BIG_MOD 1000003
|
|
|
|
struct adp {
|
|
long a;
|
|
long d;
|
|
long p;
|
|
};
|
|
|
|
long pow_mod(long a, long b, long mod, long acc) {
|
|
long powa = a % mod;
|
|
while (b > 0) {
|
|
if (b & 1) {
|
|
acc = (acc * powa) % mod;
|
|
}
|
|
powa = (powa * powa) % mod;
|
|
b >>= 1;
|
|
}
|
|
return acc;
|
|
}
|
|
|
|
long factorial_mod(long n, long mod, long acc) {
|
|
int i;
|
|
for (i = 2; i <= n; i++) {
|
|
acc = (acc * i) % mod;
|
|
}
|
|
return acc;
|
|
}
|
|
|
|
long *fenwick_make(size_t n) {
|
|
long *fenwick = calloc(n, sizeof *fenwick);
|
|
return fenwick;
|
|
}
|
|
|
|
long *fenwick_prodmake(size_t n) {
|
|
long *fenwick = malloc(n * sizeof *fenwick);
|
|
if (fenwick == NULL) {
|
|
return NULL;
|
|
}
|
|
size_t i;
|
|
for (i = 0; i < n; i++) {
|
|
fenwick[i] = 1;
|
|
}
|
|
return fenwick;
|
|
}
|
|
|
|
// Add val to element i in fenwick tree
|
|
void fenwick_update(long *fenwick, size_t n, size_t i, long val) {
|
|
i++;
|
|
if (i < 1 || i > n) {
|
|
return;
|
|
}
|
|
while (i <= n) {
|
|
fenwick[i-1] += val;
|
|
i += i & -i; // parent
|
|
}
|
|
}
|
|
|
|
// Sum of elements in [0, i[
|
|
long fenwick_prefixsum(long *fenwick, size_t n, size_t i) {
|
|
long sum = 0;
|
|
while (i > 0) {
|
|
sum += fenwick[i-1];
|
|
i -= i & -i; // predecessor
|
|
}
|
|
return sum;
|
|
}
|
|
|
|
// Sum of elements in [i, j[
|
|
long fenwick_intervalsum(long *fenwick, size_t n, size_t i, size_t j) {
|
|
return fenwick_prefixsum(fenwick, n, j) - fenwick_prefixsum(fenwick, n, i);
|
|
}
|
|
|
|
// Multiply element i by val in fenwick tree of products
|
|
void fenwick_produpdate(long *fenwick, size_t n, size_t i, long val, long mod) {
|
|
i++;
|
|
if (i < 1 || i > n) {
|
|
return;
|
|
}
|
|
while (i <= n) {
|
|
fenwick[i-1] = (fenwick[i-1] * val) % mod;
|
|
i += i & -i; // parent
|
|
}
|
|
}
|
|
|
|
// Product of elements in [0, i[
|
|
long fenwick_prefixprod(long *fenwick, size_t n, size_t i, long mod) {
|
|
long prod = 1;
|
|
while (i > 0) {
|
|
prod = (prod * fenwick[i-1]) % mod;
|
|
i -= i & -i; // predecessor
|
|
}
|
|
return prod;
|
|
}
|
|
|
|
// Product of elements in [i, j[
|
|
long fenwick_intervalprod(long *fenwick, size_t n, size_t i, size_t j, long mod) {
|
|
long p_j = fenwick_prefixprod(fenwick, n, j, mod);
|
|
long p_i = fenwick_prefixprod(fenwick, n, i, mod);
|
|
// p_i ^ -1 = p_i ^ BIG_MOD - 2 in this field
|
|
return pow_mod(p_i, BIG_MOD - 2, BIG_MOD, p_j);
|
|
}
|
|
|
|
size_t read_ints(char *line, int **ints) {
|
|
if (line == NULL) {
|
|
return 0;
|
|
}
|
|
size_t size = INTS_START_SIZE;
|
|
*ints = malloc(INTS_START_SIZE * sizeof *ints);
|
|
if (*ints == NULL) {
|
|
return 0;
|
|
}
|
|
|
|
char *tok = line;
|
|
char *numend;
|
|
size_t i = 0;
|
|
long tmp;
|
|
while (*tok != '\n' && *tok != '\0') {
|
|
if (isspace(*tok)) {
|
|
tok++;
|
|
} else if (isdigit(*tok)) {
|
|
tmp = strtol(tok, &numend, 10);
|
|
if (i == size) {
|
|
// If we have filled array, double its size
|
|
size *= 2;
|
|
*ints = realloc(*ints, size * sizeof *ints);
|
|
if (*ints == NULL) {
|
|
return 0;
|
|
}
|
|
}
|
|
// store in array
|
|
(*ints)[i] = tmp;
|
|
i++;
|
|
tok = numend;
|
|
} else {
|
|
free(*ints);
|
|
*ints = NULL;
|
|
return 0;
|
|
}
|
|
}
|
|
return i;
|
|
}
|
|
|
|
int read_adp(char *buf, int n, struct adp *adp, long *ptree, long *vtree) {
|
|
int *ints;
|
|
size_t ints_len;
|
|
|
|
int i;
|
|
char *line;
|
|
for (i = 0; i < n; i++) {
|
|
line = strtok(NULL, "\n");
|
|
ints_len = read_ints(line, &ints);
|
|
if (ints_len != 3) {
|
|
return 1;
|
|
}
|
|
adp[i].a = ints[0];
|
|
adp[i].d = ints[1];
|
|
adp[i].p = ints[2];
|
|
long v = pow_mod(adp[i].d, adp[i].p, BIG_MOD, 1);
|
|
fenwick_produpdate(vtree, n, i, v, BIG_MOD);
|
|
fenwick_update(ptree, n, i, ints[2]);
|
|
free(ints);
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
int min_const_diff(
|
|
int n,
|
|
struct adp *adp,
|
|
long *ptree,
|
|
long *vtree,
|
|
int i,
|
|
int j)
|
|
{
|
|
if (i < 1 || j > n) {
|
|
return 1;
|
|
}
|
|
long v = fenwick_intervalprod(vtree, n, i-1, j, BIG_MOD);
|
|
long k = fenwick_intervalsum(ptree, n, i-1, j);
|
|
v = factorial_mod(k, BIG_MOD, v);
|
|
printf("%ld %ld\n", k, v);
|
|
return 0;
|
|
}
|
|
|
|
int add_powers(int n, struct adp *adp, long *ptree, long *vtree, int i, int j, int v) {
|
|
if (i < 1 || j > n) {
|
|
return 1;
|
|
}
|
|
int k;
|
|
for (k = i-1; k < j; k++) {
|
|
long factor = pow_mod(adp[k].d, v, BIG_MOD, 1);
|
|
fenwick_produpdate(vtree, n, k, factor, BIG_MOD);
|
|
fenwick_update(ptree, n, k, v);
|
|
adp[k].p += v;
|
|
}
|
|
return 0;
|
|
}
|
|
|
|
int handle_query(char *str, int n, struct adp *adp, long *ptree, long *vtree) {
|
|
int *ints;
|
|
int ints_len = read_ints(str, &ints);
|
|
if (ints_len < 1) {
|
|
return 1;
|
|
}
|
|
int query_type = ints[0];
|
|
if (query_type == 0) {
|
|
if (ints_len != 3) {
|
|
return 1;
|
|
}
|
|
return min_const_diff(n, adp, ptree, vtree, ints[1], ints[2]);
|
|
}
|
|
if (query_type == 1) {
|
|
if (ints_len != 4) {
|
|
return 1;
|
|
}
|
|
return add_powers(n, adp, ptree, vtree, ints[1], ints[2], ints[3]);
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
int main(int argc, char **argv) {
|
|
char buf[BUF_SIZE];
|
|
fread(buf, BUF_SIZE, 1, stdin);
|
|
|
|
char *line = strtok(buf, "\n");
|
|
|
|
char *end;
|
|
int n = strtol(line, &end, 10);
|
|
if ((!isspace(*end) && *end != '\0') || n < 1) {
|
|
return 0;
|
|
}
|
|
|
|
struct adp *adp = malloc(n * sizeof *adp);
|
|
long *ptree = fenwick_make(n);
|
|
long *vtree = fenwick_prodmake(n);
|
|
if (adp == NULL || ptree == NULL || vtree == NULL) {
|
|
return 0;
|
|
}
|
|
int err = read_adp(buf, n, adp, ptree, vtree);
|
|
if (err) {
|
|
return 0;
|
|
}
|
|
|
|
line = strtok(NULL, "\n");
|
|
int q = strtol(line, &end, 10);
|
|
if ((!isspace(*end) && *end != '\0') || q < 1) {
|
|
return 0;
|
|
}
|
|
|
|
int i;
|
|
for (i = 0; i < q; i++) {
|
|
line = strtok(NULL, "\n");
|
|
err = handle_query(line, n, adp, ptree, vtree);
|
|
if (err) {
|
|
return 0;
|
|
}
|
|
}
|
|
|
|
return 0;
|
|
}
|