CSP-S day1

T1 临时换题（huan）
1.1 题目描述
现在都已经八点半了， pufanyi 现在因为 std 出锅很自闭，于是决定临时换题。原
来的 Day1T1 将会被移到 Day2，现在你看到的是原来的 Day2T1。
本题均使用十进制计数。
和一般的女孩子一样，小小迪喜欢教钧钧数数。
给定两个正整数 l 和 r，对于任意 x，满足 l ≤ x ≤ r，把 x 的所有约数全部写下
来。对于每个写下来的数，只保留最高位的那个数码。求 [1, 9] 中每个数码出现的次数。
这里的约数是指：对于正整数 y，如果正整数 x 使得存在正整数 k 有 kx = y， x 是
y 的约数。
1.2 输入输出格式
1.2.1 输入格式
从文件 huan.in 中读入数据。
输入一行两个正整数 l 和 r。
1.2.2 输出格式
输出到文件 huan.out 中。
输出一行 9 个非负整数，每两个数字之间用一个空格分开，第 i 个数表示数码 i 出
现的次数。
1.3 样例
1.3.1 样例 1 输入
1 4
1.3.2 样例 1 输出
4 2 1 1 0 0 0 0 0
1.3.3 样例 1 解释
数字 约数 最高位
1 1 1
2 1, 2 1, 2
3 1, 3 1, 3
4 1, 2, 4 1, 2, 4
1.3.4 样例 2 输入
233 238
1.3.5 样例 2 输出
13 10 4 2 2 1 3 0 1
1.3.6 样例 2 解释
数字 约数 最高位
233 1, 233 1, 2
234 1, 2, 3, 6, 9, 13, 18, 1, 2, 3, 6, 9, 1, 1, 2, 3, 7, 1, 2
26, 39, 78, 117, 234
235 1, 5, 47, 235 1, 5, 4, 2
236 1, 2, 4, 59, 118, 236 1, 2, 4, 5, 1, 2
237 1, 3, 79, 237 1, 3, 7, 2
238 1, 2, 7, 14, 17, 1, 2, 7, 1, 1, 3, 1, 2
34, 119, 238
1.4 数据范围与提示
本题共 10 个测试点，每个测试点 10 分。
对于所有测试点，均有 l ≤ r。
每个测试点的具体限制见下表。
测试点编号 r 特殊性质
1 ∼ 2 ≤ 10000 无
3 ∼ 5 ≤ 107
6 ∼ 7 ≤ 109 r − l ≤ 1000
8 ∼ 10 无
本题答案有 32 位整型溢出风险，故统计答案时请务必使用 64 位整型。

思路：
由题意可得 ans[getmax(i)]+=n/i
暴力的话1e7可以过直接50分，1e9会TLE
因为其实只需要考虑最高位对答案的贡献，并且我们发现在某一段区间内，对答案的贡献是一样的
比如（50,100]=1
而这个区间为 (n/k+1,n/k]，贡献值为k;
枚举首位遍历不同的区间
1->(10,19)->(100-199)
再在这个区间内根据 (n/k+1,n/k]划分区间
最后答案累加a(当前区间的值得贡献）*d(区间长度）
例：100
30-39的区间
a=3，区间33 （30,33）100/34=2,100/2=50,50和39进行比较（34，39）
所以划分为（30,33）（34,39）

#include <bits/stdc++.h>
#define ll long long
using namespace std;
inline int read(){
    int s=0,w=1;
    char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')w=-1,ch=getchar();}
    while(ch>='0'&&ch<='9'){s=s*10+ch-'0',ch=getchar();}
    return s*w;
}
ll ans[2][10],l,r;
inline void solve(int tot,int n){
    for(int i=1;i<=9;i++){
        for(int j=1;j*i<=n;j*=10){
            ll d,l=1ll*i*j,r=min(1ll*i*j+j-1,(ll)n);//d区间长度，l,r左右边界
            for(ll k=l;k<=r;k+=d){
                ll a=n/k,b=n%k;//a为区间的统一贡献值，b=n-a*k;
                d=1+min(b/a,r-k);//b/a=n/a-k,(n/a)即为断点的值，r-k防止越界
                ans[tot][i]+=a*d; 
            } 
        }
    }
}
int main(){
    freopen("huan.in", "r", stdin);
    freopen("huan.out", "w", stdout);
    scanf("%d%d", &l, &r);
    solve(0, r), solve(1, l - 1);
    for (int i = 1; i <= 9; ++i) {
        printf("%lld ", ans[0][i] - ans[1][i]);
    }
    return 0;
}