[BZOJ4026]dC Loves Number Theory

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cmath>
#include <stack>
#include <vector>
#include <queue>
#include <cstring>
#include <string>
#include <map>
#include <set>
using namespace std;

const int BufferSize = 1 << 16;
char buffer[BufferSize], *Head, *Tail;
inline char Getchar() {
	if(Head == Tail) {
		int l = fread(buffer, 1, BufferSize, stdin);
		Tail = (Head = buffer) + l;
	}
	return *Head++;
}
int read() {
	int x = 0, f = 1; char c = Getchar();
	while(!isdigit(c)){ if(c == ‘-‘) f = -1; c = Getchar(); }
	while(isdigit(c)){ x = x * 10 + c - ‘0‘; c = Getchar(); }
	return x * f;
}

#define maxn 800010
#define maxnode 12800010
#define maxp 1000010
#define maxlog 20
#define MOD 1000777
#define LL long long

int A[maxn], nn, pos[maxn], lstp[maxp], Mulsum[maxn], inv[MOD];
void gcd(int a, int b, int& x, int& y) {
	if(!b){ x = 1; y = 0; return ; }
	gcd(b, a % b, y, x); y -= a / b * x;
	return ;
}
int Inv(int a) {
	int x, y;
	gcd(a, MOD, x, y);
	return (x % MOD + MOD) % MOD;
}
int prime[maxp], cnt, val[maxp], nxt[maxp], tval[maxlog];
bool vis[maxp];
void prime_table() {
	for(int i = 2; i <= maxp - 10; i++) {
		if(!vis[i]) prime[++cnt] = i, val[i] = prime[cnt], nxt[i] = i;
		for(int j = 1; j <= cnt && i * prime[j] <= maxp - 10; j++) {
			vis[i*prime[j]] = 1;
			val[i*prime[j]] = prime[j]; nxt[i*prime[j]] = i;
			if(i % prime[j] == 0) break;
		}
	}
	return ;
}

int ToT, rt[maxn], Mul[maxnode], lc[maxnode], rc[maxnode];
void update(int& y, int x, int l, int r, int p, int v) {
	if(v > 1) Mul[y = ++ToT] = ((LL)Mul[x] * (v - 1) % MOD) * inv[v] % MOD;
	else Mul[y = ++ToT] = Mul[x];
//	printf("%d %d %d: %d(%d %d %d)\n", y, l, r, Mul[y], x, v, inv[v]);
	if(l == r) return ;
	int mid = l + r >> 1; lc[y] = lc[x]; rc[y] = rc[x];
	if(p <= mid) update(lc[y], lc[x], l, mid, p, v);
	else update(rc[y], rc[x], mid + 1, r, p, v);
	return ;
}
LL query(int o, int l, int r, int qr) {
	if(!o) return 1;
	if(r <= qr) return Mul[o];
	int mid = l + r >> 1;
	LL ans = query(lc[o], l, mid, qr);
	if(qr > mid) (ans *= query(rc[o], mid + 1, r, qr)) %= MOD;
	return ans;
}

int ANS[maxn], len;
char Out[maxn];
int main() {
//	freopen("data.in", "r", stdin);
//	freopen("data.out", "w", stdout);
	for(int i = 0; i < MOD; i++) inv[i] = Inv(i);
	prime_table();
	int n = read(), q = read();
	Mulsum[0] = 1;
	for(int i = 1; i <= n; i++) {
		int tmp = read();
		Mulsum[i] = (LL)Mulsum[i-1] * tmp % MOD;
		pos[i] = nn + 1;
		if(tmp == 1) A[++nn] = 1;
		else {
			int j = tmp, tc = 0;
			for(; nxt[j] != j; j = nxt[j]) tval[++tc] = val[j];
			tval[++tc] = val[j];
			tc = unique(tval + 1, tval + tc + 1) - tval - 1;
			for(j = 1; j <= tc; j++) A[++nn] = tval[j];
		}
	}
	pos[n+1] = nn + 1;
	n = nn; Mul[0] = 1;
	for(int i = 1; i <= n; i++) {
		update(rt[i], rt[i-1], 0, n, lstp[A[i]], A[i]);
		lstp[A[i]] = i;
	}
	/*printf("n: %d\n", n);
	for(int i = 1; i <= n; i++) printf("%d%c", A[i], i < n ? ‘ ‘ : ‘\n‘);
	for(int i = 1; i <= n; i++) printf("%d ", pos[i]); putchar(‘\n‘);*/
	int lst = 0;
	for(int i = 1; i <= q; i++) {
		int ql = read() ^ lst, qr = read() ^ lst;
		int l = pos[ql], r = pos[qr+1] - 1;
		int Multi = (LL)Mulsum[qr] * inv[Mulsum[ql-1]] % MOD;
//		printf("Multi: %d %d %d %d %d\n", Multi, rt[r], rt[l-1], query(rt[r], 0, n, l - 1), query(rt[l-1], 0, n, l - 1));
		lst = (query(rt[r], 0, n, l - 1) * inv[query(rt[l-1],0,n,l-1)] % MOD) * Multi % MOD;
		ANS[i] = lst;
//		lst = 0;
	}
	int num[10], cntn;
	for(int i = 1; i <= q; i++) {
		int tmp = ANS[i];
//		printf("%d\n", ANS[i]);
		if(!tmp) Out[len++] = ‘0‘;
		cntn = 0;
		while(tmp) num[++cntn] = tmp % 10, tmp /= 10;
		for(int j = cntn; j; j--) Out[len++] = num[j] + ‘0‘;
		if(i < q) Out[len++] = ‘\n‘; else Out[len++] = ‘\0‘;
	}
	puts(Out);
	
	return 0;
}