QOJ.ac
QOJ
ID | Problem | Submitter | Result | Time | Memory | Language | File size | Submit time | Judge time |
---|---|---|---|---|---|---|---|---|---|
#574862 | #9309. Graph | NTT | AC ✓ | 453ms | 19872kb | C++23 | 1.2kb | 2024-09-19 03:18:06 | 2024-09-19 03:18:07 |
Judging History
answer
#pragma GCC optimize("Ofast,inline,unroll-loops")
#include<bits/stdc++.h>
using ll=long long;
constexpr int B=320000,M=998244353;
auto qpow(ll a,auto b){int res=1;for(;b;a=a*a%M,b>>=1)if(b&1)res=res*a%M;return res;}
int v[B];
std::vector<int>prime;
ll n;
using Info=std::array<ll,2*B>;
auto&access(Info&a,ll x){
if(x<=B)return a[x+B];
return a[n/x];
}
signed main(){
v[1]=1;
prime.reserve(27608);
for(int i=2;i<B;++i){
if(!v[i]){
v[i]=i;
prime.emplace_back(i);
}
for(auto&p:prime){
if(i*p>=B)break;
v[i*p]=p;
if(v[i]==p)break;
}
}
std::cin>>n;
std::vector<ll>b;
Info f={};
for(ll l=1,r;l<=n;l=r+1){
ll q=n/l;
r=n/q;
b.emplace_back(q);
access(f,q)=q-1;
}
int cnt=0;
for(auto&p:prime){
for(auto&x:b){
if(p*(ll)p>x)break;
access(f,x)-=access(f,x/p)-cnt;
}
cnt++;
}
ll ans=1;
for(ll l=1,r;l<=n;l=r+1){
ll q=n/l;
r=n/q;
// std::cerr<<l<<' '<<r<<' '<<q<<'\n';
ans=(ans*qpow([&]->ll{
if(q==1)return 1;
if(q==2)return 1;
if(q==3)return 3;
// int c=pi[q]-pi[q/2];
ll c=access(f,q)-access(f,q/2);
return (q-c-1)%M*qpow(q%M,c)%M;
}()%M,r-l+1))%M;
}
printf("%lld\n",ans);
}
这程序好像有点Bug,我给组数据试试?
Details
Tip: Click on the bar to expand more detailed information
Test #1:
score: 100
Accepted
time: 0ms
memory: 10184kb
input:
4
output:
8
result:
ok answer is '8'
Test #2:
score: 0
Accepted
time: 0ms
memory: 10208kb
input:
2
output:
1
result:
ok answer is '1'
Test #3:
score: 0
Accepted
time: 0ms
memory: 10124kb
input:
123
output:
671840470
result:
ok answer is '671840470'
Test #4:
score: 0
Accepted
time: 4ms
memory: 10112kb
input:
233
output:
353738465
result:
ok answer is '353738465'
Test #5:
score: 0
Accepted
time: 0ms
memory: 10212kb
input:
5981
output:
970246821
result:
ok answer is '970246821'
Test #6:
score: 0
Accepted
time: 0ms
memory: 10124kb
input:
86422
output:
897815688
result:
ok answer is '897815688'
Test #7:
score: 0
Accepted
time: 0ms
memory: 10152kb
input:
145444
output:
189843901
result:
ok answer is '189843901'
Test #8:
score: 0
Accepted
time: 0ms
memory: 10428kb
input:
901000
output:
819449452
result:
ok answer is '819449452'
Test #9:
score: 0
Accepted
time: 0ms
memory: 10140kb
input:
1000000
output:
113573943
result:
ok answer is '113573943'
Test #10:
score: 0
Accepted
time: 0ms
memory: 10216kb
input:
23333333
output:
949849384
result:
ok answer is '949849384'
Test #11:
score: 0
Accepted
time: 8ms
memory: 10264kb
input:
102850434
output:
604886751
result:
ok answer is '604886751'
Test #12:
score: 0
Accepted
time: 15ms
memory: 10236kb
input:
998244353
output:
0
result:
ok answer is '0'
Test #13:
score: 0
Accepted
time: 15ms
memory: 10528kb
input:
1000000007
output:
318420284
result:
ok answer is '318420284'
Test #14:
score: 0
Accepted
time: 26ms
memory: 10580kb
input:
2147483547
output:
688759898
result:
ok answer is '688759898'
Test #15:
score: 0
Accepted
time: 52ms
memory: 11584kb
input:
5120103302
output:
116870489
result:
ok answer is '116870489'
Test #16:
score: 0
Accepted
time: 133ms
memory: 14424kb
input:
19834593299
output:
523663743
result:
ok answer is '523663743'
Test #17:
score: 0
Accepted
time: 267ms
memory: 14520kb
input:
52500109238
output:
195086665
result:
ok answer is '195086665'
Test #18:
score: 0
Accepted
time: 394ms
memory: 19364kb
input:
84848352911
output:
107959260
result:
ok answer is '107959260'
Test #19:
score: 0
Accepted
time: 453ms
memory: 19504kb
input:
99824435322
output:
0
result:
ok answer is '0'
Test #20:
score: 0
Accepted
time: 453ms
memory: 18272kb
input:
99999999354
output:
316301711
result:
ok answer is '316301711'
Test #21:
score: 0
Accepted
time: 437ms
memory: 19872kb
input:
100000000000
output:
396843576
result:
ok answer is '396843576'
Extra Test:
score: 0
Extra Test Passed