【题解】Atcoder ARC#90 F-Number of Digits

最新推荐文章于 2020-12-01 00:57:32 发布

转载最新推荐文章于 2020-12-01 00:57:32 发布 · 134 阅读

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/twilight-sx/p/9739078.html

本文深入探讨了AtCoder平台上的难题解决策略，通过分析特定条件下的数据处理方法，提出了有效的two-pointer算法和遍历优化方案。针对(f(l)>7)的情况，文章详细解释了如何通过设定(t=x+y)来简化问题，并通过实例代码展示了具体的实现过程。

　　Atcoder刷不动的每日一题...

　　首先注意到一个事实：随着 \(l, r\) 的增大，\(f(r) - f(l)\) 会越来越小。考虑暴力处理出小数据的情况，我们可以发现对于左端点 \(f(l) <= 7\) 的情况下，右端点的最大限度为 \(\frac{10^8}{8} + 10^7\) 。这个范围并不大，可以直接用 two-pointer 处理出来。

　　那么这部分的数据和后面的数据有什么不同呢？当 \(f(l) > 7\) 的时候，\(f(r) - f(l) <= 1\)。那么设 \(l = f(l)\)，我们有：

　　\( x * l + y * (l + 1) = S \)

令 \( t = x + y \)

　原式等于 \( t * l + y = S \)；

　　我们令 \(x + y > y\) ，即 \( x != 0 \)，那么最后一个式子实际上表达的是 \( y = S \ mod \ t \)。也就是说，对于任何的一个确定的 \(t\)，（当 \(l\) 的范围 \(> 7\)）我们都可以找到唯一对应的 \(x, y\) 与之对应。那么我们就可以扫一遍所有的 \(t\)以求得答案。注意当 \( y = 0 \) 时，对答案的贡献为这个数字范围内所有长为 \(t\) 的区间。

#include <bits/stdc++.h>
using namespace std;
#define maxn 23000000
#define int long long
#define mod 1000000007
int digit[maxn], ans, S;

int read()
{
    int x = 0, k = 1;
    char c; c = getchar();
    while(c < '0' || c > '9') { if(c == '-') k = -1; c = getchar(); }
    while(c >= '0' && c <= '9') x = x * 10 + c - '0', c = getchar();
    return x * k;
}

int Qpow(int times)
{
    int x = 10, base = 1;
    for(; times; times >>= 1, x = x * x % mod)
        if(times & 1) base = base * x % mod;
    return base;
}

void Up(int &x, int y) { x = (x + y) % mod; }
signed main()
{
    S = read(); int lim = S / 8, now = 0;
    for(int i = 1; i < maxn; i ++) digit[i] = digit[i / 10] + 1;
    for(int i = 1, now = 0, tem = 0; i < 1e7; i ++)
    {
        while(tem < S) { tem += digit[now]; now ++; }
        if(tem == S) { ans += 1; if(ans >= mod) ans -= mod; } 
        tem -= digit[i];
    }
    for(int i = 1; i <= lim; i ++)
    {
        if(!(S % i)) 
        {
            int l = S / i, sum = Qpow(l - 1) * 9 % mod;
            Up(ans, (sum - i + 1 + mod) % mod);
        }
        else 
        { 
            int l = (S - S % i) / i;
            ans += 1; if(ans >= mod) ans -= mod; 
        }
    }
    printf("%lld\n", ans);
    return 0;
}

转载于:https://www.cnblogs.com/twilight-sx/p/9739078.html