QOJ.ac
QOJ
ID | 题目 | 提交者 | 结果 | 用时 | 内存 | 语言 | 文件大小 | 提交时间 | 测评时间 |
---|---|---|---|---|---|---|---|---|---|
#244687 | #4478. Jo loves counting | ucup-team203 | TL | 0ms | 0kb | C++20 | 2.4kb | 2023-11-09 14:28:47 | 2023-11-09 14:28:48 |
answer
#include <bits/stdc++.h>
const int N = 2e6 + 10;
using namespace std;
#define ll __int128_t
const ll Mod = (1ll << 57) * 29ll + 1;
ll Pow(ll n, ll k)
{
ll res = 1;
for (ll x = n; k; x = x * x % Mod, k >>= 1)
if (k & 1)
res = res * x % Mod;
return res;
}
ll isp[N], p[N], tot;
ll M;
ll ans = 0,inv2;
struct node
{
/* data */
ll now, pos, h;
};
void calc()
{
stack<node> st;
st.push((node){1, 0, 1});
while (!st.empty())
{
node tmp=st.top();
st.pop();
ll now=tmp.now,pos=tmp.pos,h=tmp.h;
// printf("%lld %lld %lld\n",(long long) pos,(long long) M,(long long) p[pos]*p[pos]);
if (p[pos] * p[pos] > M / now)
{
if (now == 1)
{
ans += (M * (M - 1)%Mod *inv2) % Mod;
}
else
{
ans += (h > 0 ? 1 : -1) * now * Pow((h > 0 ? h : -h), Mod - 2) % Mod * (M / now) % Mod * (M / now + 1) * inv2 % Mod;
ans = (ans % Mod + Mod) % Mod;
// printf("%lld %lld\n",(long long)now,(long long)(h>0?1:-1)*now*Pow(abs(h),Mod-2)%Mod);
}
// printf("%lld\n",(long long)ans);
continue;
}
// calc(now, pos + 1, h);
// if (p[pos] * p[pos] > M / now)continue;
st.push((node){now, pos + 1, h});
// if(p[pos]==0)
// {
// printf("%lld"p[pos-1]);
// }
for (ll i = 2, pj = p[pos] * p[pos]; now <= M / pj; i++, pj *= p[pos])
{
// calc(now * pj, pos + 1, h * (-i * (i - 1)));
st.push((node){now * pj, pos + 1, h * (-i * (i - 1))});
}
}
}
void solve()
{
long long _;
cin >> _;
M = _;
calc();
ans = ans * inv2 % Mod;
printf("%lld\n", (long long)ans);
}
int main()
{
inv2=Pow(2,Mod-2);
for (int i = 2; i <= N - 10; i++)
{
if (!isp[i])
p[tot++] = i;
for (int j = 0; i * p[j] < N && j < tot; j++)
{
isp[i * p[j]] = 1;
if (i % p[j])
;
else
break;
}
}
int T;
cin >> T;
while (T--)
{
solve();
}
}
/*
12
4
959139
9859
125987
209411
965585325539
213058376259
451941492387
690824608515
934002691939
1000000000000
1
*/
详细
Test #1:
score: 0
Time Limit Exceeded
input:
12 4 959139 9859 125987 209411 965585325539 213058376259 451941492387 690824608515 934002691939 1000000000000 1
output:
2 1063521647235962441 1654756389743847357 28123286989126918 186817287120753871 509981696528294966 224839384346900823 2717574651366567155 1474379426507740960 330832676794761953