-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtemplate.cpp
117 lines (93 loc) · 2.11 KB
/
template.cpp
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
#include <bits/stdc++.h>
#define ll long long
#define int ll
#define pii pair<ll, ll>
#define ff first
#define ss second
#define mp make_pair
#define pb push_back
using namespace std;
ll MOD = 998244353;
ll pwr(ll x, ll y, ll mod = MOD) {
ll res = 1;
x = x%mod;
while (y > 0) {
if (y&1) res = (res*x)%mod;
y = y>>1;
x = (x*x)%mod;
}
return res;
}
inline ll addmod(ll a, ll b){
return ((a+b)%MOD);
}
inline ll mulmod(ll a, ll b){
return ((a*b)%MOD);
}
bool check_composite(int n, int a, int d, int s) {
int x = pwr(a, d, n);
if (x == 1 || x == n - 1)
return false;
for (int r = 1; r < s; r++) {
x = x * x % n;
if (x == n - 1)
return false;
}
return true;
};
const int first12prime[12] = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37};
bool isPrime(int n) { // O(log n) time complexity
if (n < 2)
return false;
int r = 0;
int d = n - 1;
while ((d & 1) == 0) {
d >>= 1;
r++;
}
for (int i = 0; i < 12; i++) {
int a = first12prime[i];
if (n == a)
return true;
if (check_composite(n, a, d, r))
return false;
}
return true;
}
const int PM = 2e5+5;
vector<int> primes;
vector<bool> is_comp;
void sieve () {
is_comp.resize(PM, false);
for (int i = 2; i < PM; i++) {
if (!is_comp[i]) primes.pb(i);
for (int j = 0; j < primes.size() && primes[j]*i < PM; j++) {
is_comp[primes[j]*i] = true;
if (i % primes[j] == 0) break;
}
}
}
inline int sfunc(int a, int b) {
return a+b;
}
int n; // array size
int *stree; // stree = new int[2 * n]
void build() { // build the tree
for (int i = n - 1; i > 0; --i) stree[i] = sfunc(stree[i<<1], stree[i<<1|1]);
}
void modify(int p, int value) { // set value at position p
for (stree[p += n] = value; p > 1; p >>= 1) stree[p>>1] = sfunc(stree[p], stree[p^1]);
}
int query(int l, int r) { // query on interval [l, r)
int res = 0; // DEFAULT
for (l += n, r += n; l < r; l >>= 1, r >>= 1) {
if (l&1) res = sfunc(res, stree[l++]);
if (r&1) res = sfunc(res, stree[--r]);
}
return res;
}
signed main () {
ios_base::sync_with_stdio(false); cin.tie(NULL); cout.tie(NULL);
int t; cin >> t;
while (t--) {}
}