2024HBCPC：J Iris‘ Food

题目描述

猫猫 Iris 又在玩 oql 的键盘了。
在 T 天的时间里，每一天 Iris 会在玩 oql 键盘的时候敲下若干个字符，而这些字符恰好全都是 $0\sim 9$ 这 $10$ 个阿拉伯数字。经过统计，数字 $i$ 有 $a_i$个。oql 每天可能会给 Iris 喂一定数量的猫粮。他决定，在这 ${\sum }_{i=0}^9{a}_i$ 个数字里选择 $m$ 个，经过重新排列后形成一个十进制下 $m$ 位、且不包含前导 0 的非负整数 x，然后 oql 将在这天给 Iris 投喂 $x$ 克的猫粮。
然而，Iris 还是一只小猫，不宜食用太多的猫粮。因此 oql 想让这个数字尽可能小。请你帮助 Iris 计算出她每天能得到多少猫粮。由于答案可能很大，你只需要输出答案对 $10^9+7$ 取模的结果。请注意，你需要输出的是最小答案模 $10^9+7$ 后的结果，而不是 $x$ $mod$ $10^9+7$ 的最小值。

Input

第一行一个正整数 $T(1\le T\le 10^4)$，表示天数。
第 $2\sim (T+1)$ 行，每行 $11$ 个由空格分隔的非负整数 $m,a_0,a_1,\cdots ,a_9(1\le m\le 10^9,0\le a_i\le 10^9)$，表示第 $i$ 天的情况。
数据保证有解，即至少能形成一个 $m$ 位不包含前导 $0$ 的非负整数。

Output

对于每一天，输出一行一个整数，表示 Iris 这一天能得到的猫粮对 $10^9+7$ 取模后的结果。

Sample Input 1

3
3 1 0 0 0 3 0 0 0 0 0
1 2 0 0 0 0 0 0 0 0 0
4 0 1 1 1 3 0 0 0 0 0

Sample Output 1

404
0
1234

解题思路

很简单的一个贪心，因为不能有前导零，所以可以先从 $1\sim 9$ 中拿出一个最小的数作为第一个数，之后就可以从 $0\sim 9$ 从小往大依次选择 $m-1$ 个数即可

接下来的问题就是高精度取余了，由同余定理可知：

假设一个数$x=abcd=a\ast 1000+b\ast 100+c\ast 10+d$
$=(((a\ast 10+b)\ast 10+c)\ast 10+d)$
那么 $x\%p=((((((a\%p)\ast 10+b)\%p)\ast 10+c)\%p)\ast 10+d)\%p$

所以，对于一个数字字符串 $S$，求它取余 $p$，就是如下实现

long long res = 0;
for (char c : S)
{
    int x = c - '0';
    res = (res * 10 + x) % p;
}
return res;

那么，如何计算答案呢？
假设 $a=xxxxxx$，如何求$b=xxxxxxyyyyyyyymodp$ 的值
我们可得到下列式子

$b=a\ast 10^8+yyyyyyyy$
$b\%p=(a\ast 10^8\%p+yyyyyyyy\%p)\%p$

现在，问题都解决了，我们只需要预处理出一个数组 $cnt$，$cnt[i]$表示长度为 $i$ 的全 $1$ 数字 $mod10^9+7$ 的值，上面的式子就可以得到以下变换:

$b\%p=(a\ast 10^8\%p+(11111111\ast y)\%p)\%p$

最后，注意特判 $m=1$ 且 $a[0]\ge 1$时，答案为 $0$

实现代码

#pragma GCC optimize(3, "Ofast", "inline")
#include<bits/stdc++.h>
#define int long long
#define x first
#define y second
#define debug(x) cerr << #x" = " << x << '\n';
using namespace std;

const int mod = 1e9 + 7;
int cnt[10000010];

int qmi(int a, int b)
{
    int res = 1;
    while (b)
    {
        if (b & 1) res = res * a % mod;
        b >>= 1;
        a = a * a % mod;
    }
    return res;
}

void init()
{
    for (int i = 1; i <= 10000000; i++) // 预处理出10^7个1 % p的值
    {
        cnt[i] = cnt[i - 1] * 10 + 1;
        cnt[i] %= mod;
    }
}

void solve()
{
    int m;
    cin >> m;
    vector<int> a(10);
    for (int i = 0; i < 10; i++) cin >> a[i];
    int res = 0;
    if (m == 1 && a[0]) return cout << 0 << '\n', void();
    for (int i = 0; i < 10; i++)
    {
        if (a[i] && i)
        {
            res = res * 10 + i;
            res %= mod;
            m --;
            a[i] --;
            break;
        }
    }
    // 如果发现只能取0，则直接输出0
    if (res == 0) return cout << 0 << '\n', void();
    for (int i = 0; i < 10; i++)
    {
        int x = min(m, a[i]);
        m -= x;
        while (x >= 10000000)
        {
            res = (res * qmi(10, 10000000) % mod + cnt[10000000] * i % mod) % mod;
            x -= 10000000;
        }
        if (x) res = (res * qmi(10, x) % mod + cnt[x] * i % mod) % mod;
    }
    cout << res << '\n';
}

signed main()
{
    // freopen("Sample.in", "r", stdin);
    ios::sync_with_stdio(false);
    cin.tie(nullptr);
    init();
    int T = 1;
    cin >> T;
    while (T--) solve();
    return 0;
}