牛客练习赛89E-牛牛小数点【数论】

1|0正题

题目链接:https://ac.nowcoder.com/acm/contest/11179/E

---

1|1题目大意

定义$f(x)$表示$\frac{1}{x}$的混循环节长度（如果没有循环节就是$0$），$T$组询问给出$l,r$求

$$\sum_{i=l}^rf(i)$$

$1\leq T\leq 100,1\leq l\leq r\leq 10^{15}$

---

1|2解题思路

设$a_i$表示$i$位之后的余数，那么出现循环节当且仅当$a_i$出现重复。

$$a_i=a_{i-1}\times 10\% n\Rightarrow a_i=10^i\%n$$

那么出现循环节一定满足存在一个正整数$k$使得

$$a_i\times 10^k\equiv a_i(mod\ n)$$

$$\Rightarrow 10^k\equiv 1(mod\ \frac{n}{gcd(a_i,n)})$$

我们知道$10^k\equiv 1(mod\ n)$有解当且仅当$gcd(10,n)=1$。

也就是说我们要找到第一个$i$使得$gcd(10,\frac{n}{gcd(a_i,n)})=1$。

而$a_i$每次乘十，所以相当于$n$每次在质因数中去掉一个$2$和$5$，直到和$10$互质。

但是这样还没有结束，因为如果没有循环节就是$0$，而这里则会统计小数的长度，得减去这些情况，不难发现有循环节的话当且仅当存在某个$10^k\%n=0$，也就是说$n$只由$2$和$5$构成，暴力枚举这些数就好了。

时间复杂度：$O(T\log_2n\log_5 n)$

---

1|3code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll P=998244353;
ll T,l,r;
int Ask(ll n){
    ll ans=n;
    for(ll i=1,pw=2;pw<=n;i++,pw=pw*2ll)ans=(ans+n/pw)%P;
    for(ll i=1,pw=5;pw<=n;i++,pw=pw*5ll)ans=(ans+n/pw)%P;
    for(ll i=1,pw=10;pw<=n;i++,pw=pw*10ll)ans=(ans-n/pw)%P;
    for(ll i=1,pw=1;pw<=n;i++,pw=pw*2ll)
        for(ll j=1,qw=1;pw*qw<=n;j++,qw=qw*5ll)
            (ans-=max(i,j))%=P;
    return (ans+P)%P;
}
signed main()
{
    scanf("%lld",&T);
    while(T--){
        scanf("%lld%lld",&l,&r);
        printf("%lld\n",(Ask(r)-Ask(l-1)+P)%P);
    }
    return 0;
}

__EOF__

本文作者：QuantAsk
本文链接：https://www.cnblogs.com/QuantAsk/p/15333340.html
关于博主：退役OIer，GD划水选手
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！