QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#579280 | #9309. Graph | Bob_Wang | TL | 1844ms | 32016kb | C++17 | 1.8kb | 2024-09-21 11:23:39 | 2024-09-21 11:23:40 |
Judging History
answer
#include<cstdio>
#include<cmath>
#include<iostream>
#include<unordered_map>
#define ll long long
#define maxn 1000005
using namespace std;
const int M=1000000;
const ll mod=998244353;
ll n,m;
ll tot,prime[maxn],mark[maxn];
ll cnt,id1[maxn],id2[maxn];//将区间映射到连续的数组
ll g[maxn],v[maxn],pre[maxn];
unordered_map<ll,ll>mp;
ll mul(ll x,ll y)
{
if(y<0)
return 1;
x%=mod;
ll ret=1;
while(y)
{
if(y&1)
ret=ret*x%mod;
x=x*x%mod;
y>>=1;
}
return ret;
}
void ready(ll m)
{
for(int i=2;i<=m;++i)
{
if(!mark[i])
prime[++tot]=i;
for(int j=1;j<=tot;++j)
{
ll x=i*prime[j];
if(x>M)
break;
mark[x]=1;
if(i%prime[j]==0)
break;
}
}
}
ll getnum(ll n,ll x)
{
return (x<=m)?id1[x]:id2[n/x];
}
ll f(ll x)//求f(x)
{
if(x==1)
return 0;
else
return 1;
}
ll sumf(ll x)//求f(2)+...+f(x)的和
{
return x-1;
}
void pre_solve(ll n)//预处理g数组
{
for(ll l=1,r;l<=n;l=r+1)
{
r=n/(n/l);
ll val=n/l;
if(val<=m)
id1[val]=++cnt;
else id2[n/val]=++cnt;
v[cnt]=val;
g[cnt]=sumf(val);
}
for(ll j=1;j<=tot;++j)
{
ll x=prime[j];
for(ll i=1;i<=cnt&&x*x<=v[i];++i)
g[i]-=g[getnum(n,v[i]/x)]-(j-1);
}
}
ll getprime(ll n)
{
if(n<M)
return pre[n];
if(mp.count(n)!=0)
return mp[n];
tot=cnt=0;
m=sqrt(n);
ready(m);
pre_solve(n);
return mp[n]=g[getnum(n,n)];
}
int main()
{
ready(M);
for(int i=2;i<=M;++i)
pre[i]=pre[i-1]+(1-mark[i]);
scanf("%lld",&n);
ll ans=1;
for(ll l=1,r;l<=n;l=r+1)
{
r=n/(n/l);
ll cnt=n/l,tot_p=(getprime(cnt)-getprime(cnt>>1))%mod;
ll del;
if(cnt>3)
del=(cnt-tot_p-1)%mod*mul(cnt,tot_p)%mod;
else del=mul(cnt,cnt-2);
ans=ans*mul(del,r-l+1)%mod;
}
printf("%lld",ans);
return 0;
}
详细
Test #1:
score: 100
Accepted
time: 3ms
memory: 21124kb
input:
4
output:
8
result:
ok answer is '8'
Test #2:
score: 0
Accepted
time: 4ms
memory: 21952kb
input:
2
output:
1
result:
ok answer is '1'
Test #3:
score: 0
Accepted
time: 4ms
memory: 23000kb
input:
123
output:
671840470
result:
ok answer is '671840470'
Test #4:
score: 0
Accepted
time: 4ms
memory: 21772kb
input:
233
output:
353738465
result:
ok answer is '353738465'
Test #5:
score: 0
Accepted
time: 9ms
memory: 23192kb
input:
5981
output:
970246821
result:
ok answer is '970246821'
Test #6:
score: 0
Accepted
time: 4ms
memory: 21404kb
input:
86422
output:
897815688
result:
ok answer is '897815688'
Test #7:
score: 0
Accepted
time: 4ms
memory: 23220kb
input:
145444
output:
189843901
result:
ok answer is '189843901'
Test #8:
score: 0
Accepted
time: 4ms
memory: 22256kb
input:
901000
output:
819449452
result:
ok answer is '819449452'
Test #9:
score: 0
Accepted
time: 4ms
memory: 26996kb
input:
1000000
output:
113573943
result:
ok answer is '113573943'
Test #10:
score: 0
Accepted
time: 12ms
memory: 27096kb
input:
23333333
output:
949849384
result:
ok answer is '949849384'
Test #11:
score: 0
Accepted
time: 40ms
memory: 29760kb
input:
102850434
output:
604886751
result:
ok answer is '604886751'
Test #12:
score: 0
Accepted
time: 356ms
memory: 27472kb
input:
998244353
output:
0
result:
ok answer is '0'
Test #13:
score: 0
Accepted
time: 354ms
memory: 28900kb
input:
1000000007
output:
318420284
result:
ok answer is '318420284'
Test #14:
score: 0
Accepted
time: 767ms
memory: 31748kb
input:
2147483547
output:
688759898
result:
ok answer is '688759898'
Test #15:
score: 0
Accepted
time: 1844ms
memory: 32016kb
input:
5120103302
output:
116870489
result:
ok answer is '116870489'
Test #16:
score: -100
Time Limit Exceeded
input:
19834593299