QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#570103 | #9309. Graph | u2x1 | AC ✓ | 856ms | 76188kb | C++17 | 3.3kb | 2024-09-17 13:49:58 | 2024-09-17 13:50:06 |
Judging History
answer
#include<bits/stdc++.h>
#include<ext/pb_ds/assoc_container.hpp>
#include<ext/pb_ds/hash_policy.hpp>
#define ll long long
#define all(v) v.begin(),v.end()
#define sz(v) ((ll)v.size())
#define V vector
#define vi V<int>
#define vll V<ll>
#define eb emplace_back
#define pii pair<int,int>
#define pll pair<ll,ll>
#define fi first
#define se second
#define A array
#define pb push_back
#define mset multiset
#define umap unordered_map
#define ui unsigned int
#define ull unsigned ll
#define i128 __int128
#define cerr if (test) cerr
#define freopen if (test) freopen
#define whtest if (test)
using namespace std;
const int test = 0;
namespace jhsy {
constexpr int P = 998244353,N = 1e6,BUC = 999983;
template<class T,class U,int N,int BUC,T(*hash)(T)=nullptr>
struct HashMap {
struct Node {
int nxt; T key; U val;
};
A<Node,N> a;
A<int,BUC> head,col{};
int sz{},ver{1};
HashMap() {}
void clear() {
sz = 0; ver++;
}
auto end() {return a.end();}
auto find(const T &key) {
const int id = hash(key);
if (col[id] < ver) {
return a.end();
}
for (int i = head[id]; i != -1; i = a[i].nxt) {
if (key == a[i].key) {
return a.begin()+i;
}
}
return a.end();
}
auto ins(const T &key,const U &val = {}) {
const int id = hash(key);
if (col[id] < ver) {
head[id] = -1; col[id] = ver;
}
int x = sz++;
a[x] = {head[id],key,val};
head[id] = x;
return a.begin()+x;
}
auto& operator[] (const T key) {
auto it = find(key);
if (it != end()) {
return it->val;
}
return ins(key)->val;
}
};
vi p;
void init(const int n = N) {
A<bool,N+1> isp; isp.fill(1);
for (int i = 2; i <= n; i++) {
if (isp[i]) {
p.eb(i);
}
for (auto x:p) {
if (i*x > n) {
break;
}
isp[i*x] = 0;
if (!(i%x)) {
break;
}
}
}
}
ll power(ll b,ll k) {
ll a = 1;
b = (b%P+P)%P;
k = (k%(P-1)+(P-1))%(P-1);
for (; k; b = b*b%P, k >>= 1) {
if (k&1) {
a = a*b%P;
}
}
return a;
}
void main() {
ll n;
cin >> n;
ll ans = 1;
auto hash = [](ll x) {
return x % BUC;
};
HashMap<ll,ll,N*2,N*2,hash> mp;
{
vll vec;
for (ll lst = 0,cur; lst < n; lst = cur) {
ll k = n/(lst+1); cur = n/k;
vec.eb(cur);
}
reverse(all(vec));
for (auto x:vec) {
mp[x] = x-1;
}
for (int k = 0; k < sz(p); k++) {
ll t = 1ll*p[k]*p[k];
if (t > n) {
break;
}
for (auto x:vec) {
if (t > x) {
break;
}
mp[x] -= mp[x/p[k]]-k;
}
}
}
auto get_p = [&](ll x) -> ll {
return mp[x];
};
auto get = [&](ll x) {
ll res = get_p(x)-get_p(x/2);
if (x/2 < 2 && x >= 2) {
res--;
}
return res;
};
auto calc = [&](ll x) -> ll {
if (x == 1) {
return 1;
}
ll p = get(x)+1,q = x-p;
return q%P*power(p+q,p-1)%P;
};
for (ll lst = 0,cur; lst < n; lst = cur) {
ll k = n/(lst+1); cur = n/k;
ans = ans*power(calc(k),cur-lst)%P;
}
cout << ans << '\n';
}
}
int main() {
freopen("test.in", "r", stdin);
freopen("test.out", "w", stdout);
ios::sync_with_stdio(false);
cin.tie(nullptr);
cout << fixed << setprecision(20);
jhsy::init();
int T = 1;
// cin >> T;
while (T--) {
jhsy::main();
}
return 0;
}
这程序好像有点Bug,我给组数据试试?
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 66504kb
input:
4
output:
8
result:
ok answer is '8'
Test #2:
score: 0
Accepted
time: 4ms
memory: 66508kb
input:
2
output:
1
result:
ok answer is '1'
Test #3:
score: 0
Accepted
time: 0ms
memory: 66276kb
input:
123
output:
671840470
result:
ok answer is '671840470'
Test #4:
score: 0
Accepted
time: 13ms
memory: 66192kb
input:
233
output:
353738465
result:
ok answer is '353738465'
Test #5:
score: 0
Accepted
time: 11ms
memory: 66204kb
input:
5981
output:
970246821
result:
ok answer is '970246821'
Test #6:
score: 0
Accepted
time: 11ms
memory: 66208kb
input:
86422
output:
897815688
result:
ok answer is '897815688'
Test #7:
score: 0
Accepted
time: 4ms
memory: 66208kb
input:
145444
output:
189843901
result:
ok answer is '189843901'
Test #8:
score: 0
Accepted
time: 12ms
memory: 66276kb
input:
901000
output:
819449452
result:
ok answer is '819449452'
Test #9:
score: 0
Accepted
time: 7ms
memory: 66208kb
input:
1000000
output:
113573943
result:
ok answer is '113573943'
Test #10:
score: 0
Accepted
time: 3ms
memory: 66324kb
input:
23333333
output:
949849384
result:
ok answer is '949849384'
Test #11:
score: 0
Accepted
time: 11ms
memory: 66436kb
input:
102850434
output:
604886751
result:
ok answer is '604886751'
Test #12:
score: 0
Accepted
time: 27ms
memory: 66620kb
input:
998244353
output:
0
result:
ok answer is '0'
Test #13:
score: 0
Accepted
time: 27ms
memory: 66700kb
input:
1000000007
output:
318420284
result:
ok answer is '318420284'
Test #14:
score: 0
Accepted
time: 44ms
memory: 67196kb
input:
2147483547
output:
688759898
result:
ok answer is '688759898'
Test #15:
score: 0
Accepted
time: 85ms
memory: 68076kb
input:
5120103302
output:
116870489
result:
ok answer is '116870489'
Test #16:
score: 0
Accepted
time: 213ms
memory: 71328kb
input:
19834593299
output:
523663743
result:
ok answer is '523663743'
Test #17:
score: 0
Accepted
time: 494ms
memory: 71712kb
input:
52500109238
output:
195086665
result:
ok answer is '195086665'
Test #18:
score: 0
Accepted
time: 736ms
memory: 76124kb
input:
84848352911
output:
107959260
result:
ok answer is '107959260'
Test #19:
score: 0
Accepted
time: 856ms
memory: 75608kb
input:
99824435322
output:
0
result:
ok answer is '0'
Test #20:
score: 0
Accepted
time: 830ms
memory: 75068kb
input:
99999999354
output:
316301711
result:
ok answer is '316301711'
Test #21:
score: 0
Accepted
time: 836ms
memory: 76188kb
input:
100000000000
output:
396843576
result:
ok answer is '396843576'
Extra Test:
score: 0
Extra Test Passed