Tuyển sinh 🎓 Mình nhận đào tạo lập trình giải thuật thi đấu: thi vào lớp 10 chuyên tin · thi học sinh giỏi tin học cấp tỉnh/thành · thi vào đội tuyển quốc gia · Olympic tin học sinh viên / ICPC · phỏng vấn thuật toán. Bấm vào đây để xem thông tin khóa học → Xem khóa học

Đếm Chữ Số

Đếm Chữ Số
Nguồn: Codeforces Group Độ khó: Dễ Mảng đếm Tính trước

Đề bài

Bạn An đánh số trang cho một quyển sách có $N$ trang. Sau khi đánh số xong, bạn ấy hỏi các bạn các câu hỏi sau:

  • Có bao nhiêu chữ số 0 mà bạn ấy đã sử dụng?
  • Có bao nhiêu chữ số 1 mà bạn ấy đã sử dụng?
  • …
  • Có bao nhiêu chữ số 9 mà bạn ấy đã sử dụng?

Input

  • Dòng đầu tiên là $T$ – số test.
  • $T$ dòng tiếp theo, mỗi dòng là một số nguyên $N$.
  • Chúng ta không quan trọng câu chuyện $N$ chẵn hay $N$ lẻ nhé.

Output

  • In ra $T$ dòng, mỗi dòng gồm 10 số là số lượng chữ số 0, chữ số 1, …, chữ số 9 được sử dụng.

Giới hạn

  • $1 \le T \le 1000$
  • $1 \le N \le 10^6$

Subtask:

  • 50% số test đầu tiên có $T \le 100$, $N \le 10^4$.
  • 50% số test còn lại không có giới hạn gì thêm.
Ví dụ – input
1
1
Ví dụ – output
0 1 0 0 0 0 0 0 0 0

Lời mở đầu

Bài này mình được một bạn học viên hỏi nên sẵn tiện mình viết bài giải chi tiết để nhiều bạn được tham khảo luôn. Hy vọng mang đến được giá trị cho các bạn 😇. Bạn nào muốn thử nộp bài để hệ thống chấm xem có đúng hay không thì link bài đây nhé:

Dành cho học viên · Ưu tiênNộp bài trên hệ thống CLUE OJ

Link này chỉ học viên đã tham gia tổ chức mới vào được. Học viên nên ưu tiên nộp ở đây vì sẽ được chấm nhanh hơn.

Nộp bài trên CLUE OJ
Dành cho tất cả mọi ngườiNộp bài trên Codeforces

Bạn nào muốn thử nộp bài để hệ thống chấm xem có đúng hay không thì dùng link này nhé. Chưa có tài khoản thì tạo tài khoản rồi đăng nhập nộp bài thử hen.

Nộp bài trên Codeforces

Lộ trình lời giải, bấm vào từng phần để nhảy tới:

Đề bàiHiểu đúng đề: có T bộ test

Trước tiên mình giải thích rõ đề bài vì mình thấy nhiều bạn đọc không kỹ cứ làm nộp sai, đó là ta có $T$ bộ test các bạn nhé, và kết quả sẽ xuất ra $T$ dòng, mỗi dòng có 10 số là số lượng các chữ số từ 0 đến 9 nha các bạn. Ví dụ với bộ test input như sau thì kết quả sẽ xuất ra là:

Ví dụ – input
3
5
10
15
Ví dụ – output
0 1 1 1 1 1 0 0 0 0
1 2 1 1 1 1 1 1 1 1
1 8 2 2 2 2 1 1 1 1
Chữ số0123456789
N = 50111110000
N = 101211111111
N = 151822221111

Giải thích: Dòng đầu tiên $N = 5$ thì từ 1 tới 5 ta sẽ có các chữ số 1, 2, 3, 4, 5 nên tại vị trí các chữ số này sẽ có số lượng là 1, còn các chữ số còn lại sẽ có số lượng là 0, nên ta có dòng kết quả tương ứng là: 0 1 1 1 1 1 0 0 0 0.

Dòng tiếp theo $N = 10$ thì từ 1 tới 10 sẽ là: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10. Từ 1 tới 9 sẽ là các chữ số từ 1 tới 9 đều xuất hiện 1 lần, đến số 10 thì chữ số 1 được xuất hiện thêm 1 lần nữa thành ra 2 lần, và có thêm xuất hiện của chữ số 0. Nên tóm lại với $N = 10$ thì có đủ các chữ số từ 0 đến 9 xuất hiện 1 lần, riêng chữ số 1 xuất hiện 2 lần, nên ta có dòng kết quả tương ứng là: 1 2 1 1 1 1 1 1 1 1.

Dòng tiếp theo $N = 15$ thì tương tự như mình đã phân tích ở 2 trường hợp $N = 5$, $N = 10$ ở trên các bạn tự suy ra hen (gợi ý: các số 11, 12, 13, 14, 15 góp thêm 6 chữ số 1 và mỗi chữ số 2, 3, 4, 5 thêm 1 lần).

Tại vì đề bài đưa ra ví dụ với $T = 1$ test case nên nhiều bạn chưa rành về vụ có nhiều bộ test thế này, vì thế mình thấy các bạn nộp sai nhiều.

Cách 1Vét cạn từng bộ test – O(T · N · log N)

Giờ ta bàn về cách giải quyết cho bài này. Thì một cách đơn giản nhất mà nhiều bạn sẽ nhìn thấy và làm đó là cứ cho vòng lặp chạy $T$ bộ test, với mỗi bộ test đọc số $N$ vào và chạy vòng lặp với giá trị $i$ từ 1 tới $N$, mỗi giá trị $i$ sẽ đi phân tích và đếm số lần xuất hiện của các chữ số từ 0 đến 9 của giá trị $i$ đó, và cứ thế cộng dồn hết lại; sau cùng khi vòng lặp $i$ kết thúc ta in ra kết quả số lượng của từng chữ số. Rồi cứ thế quy trình ở trên lặp lại qua bộ test tiếp theo.

Mảng đếm 10 phần tử thay cho 10 biến

Nói sâu hơn về việc mỗi giá trị $i$ cần đếm xem có bao nhiêu chữ số từ 0 đến 9 xuất hiện trong đó. Thì chỗ này bạn nào mà không có kinh nghiệm thì sẽ làm cách tạo ra 10 biến đếm tương ứng với 10 chữ số từ 0 đến 9, rồi sau đó phải xét 10 trường hợp if else (vd: if(digit == 0) cnt0++; else if(digit == 1) cnt1++; …). Như thế vẫn sẽ đúng nhưng rất dài dòng.

Một cách gọn nhất trong trường hợp này là ta sẽ tạo ra mảng có 10 phần tử để index tương ứng từ 0 đến 9, và nó cũng tương ứng là các chữ số từ 0 đến 9 luôn. Như thế a[i] hiểu là số lần xuất hiện của chữ số $i$. Ví dụ a[3] = 5 thì ta hiểu chữ số 3 xuất hiện 5 lần. Ban đầu mảng 10 phần tử khởi tạo tất cả giá trị là 0. Rồi cứ mỗi chữ số phân tích ra được từ số $i$ thì sẽ gọi a[digit]++;, tức là tăng thêm 1 tại vị trí a[digit]. Sau cùng thì ta duyệt vòng lặp từ 0 đến 9 và in ra các giá trị của mảng tại các vị trí đó thì đó là dãy đáp án.

Tách từng chữ số của một số

Về việc với từng số $i$ tiến hành đi phân tích ra từng chữ số của nó thì cái này quá cơ bản nhưng mình cũng nói luôn kẻo bạn nào không biết: với 1 số thì ta chia lấy dư cho 10 để lấy ra được chữ số cuối cùng của nó, ví dụ: 123 % 10 = 3. Rồi để lấy ra được chữ số tiếp theo thì ta phải loại bỏ đi chữ số mà ta vừa lấy ra bằng cách lấy số đó chia nguyên cho 10, ví dụ: 123 / 10 = 12 (vì cả 2 đều đang là số nguyên nên phép chia 2 số nguyên cho nhau kết quả sẽ ra số nguyên, chứ nếu 1 trong 2 mà là số thực thì kết quả phép chia sẽ là 12.3).

Tách chữ số của 123 (lấy dần từ cuối lên)
123% 10 = 3/ 10 = 12
12% 10 = 2/ 10 = 1
1% 10 = 1/ 10 = 0
0x = 0 → dừng
Lấy được các chữ số 3, 2, 1 → a[3]++, a[2]++, a[1]++.

Rồi số còn lại là 12 thì ta lại lặp lại quy trình trên: 12 % 10 = 2, 12 / 10 = 1, 1 % 10 = 1, 1 / 10 = 0. Đến khi số là 0 rồi thì ta dừng quy trình vì từ bây giờ nó sẽ cứ lặp vô tận 0 % 10 = 0, 0 / 10 = 0, …

Source code để các bạn tham khảo với hướng làm mình trình bày ở trên:

#include <bits/stdc++.h>
using namespace std;

int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    int t;
    cin >> t;
    while(t--){
        int n;
        cin >> n;
        int a[10] = {0};
        for(int i = 1; i <= n; ++i){
            int x = i;
            while(x != 0){
                a[x % 10]++;
                x /= 10;
            }
        }
        for(int i = 0; i < 10; ++i){
            cout << a[i] << " ";
        }
        cout << endl;
    }
    return 0;
}
import sys
input = sys.stdin.readline

t = int(input())
for _ in range(t):
    n = int(input())
    a = [0] * 10
    for i in range(1, n + 1):
        x = i
        while x != 0:
            a[x % 10] += 1
            x //= 10
    print(*a)
C++▶ Chạy tay từng dòng codeCách 1 (vét cạn): mỗi bộ test đếm lại chữ số của 1..N bằng % 10 và / 10, O(T · N · log N).
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3
4int main(){
5 ios::sync_with_stdio(0);
6 cin.tie(0);
7 int t;
8 cin >> t;
9 while(t--){
10 int n;
11 cin >> n;
12 int a[10] = {0};
13 for(int i = 1; i <= n; ++i){
14 int x = i;
15 while(x != 0){
16 a[x % 10]++;
17 x /= 10;
18 }
19 }
20 for(int i = 0; i < 10; ++i){
21 cout << a[i] << " ";
22 }
23 cout << endl;
24 }
25 return 0;
26}

Giải thích một số chỗ trong code

  1. 2 dòng đầu tiên của hàm main (ios::sync_with_stdio(false); cin.tie(0);): hiểu cơ bản là nó tăng tốc luồng nhập xuất của C++ (tắt đồng bộ với C và bỏ ràng buộc giữa cin và cout), giúp cin/cout chạy nhanh ngang cú pháp nhập xuất của C. Một kinh nghiệm trong các bài cần nhập xuất nhiều, ví dụ như bài này ta thấy có tối đa $T = 1000$ bộ test, nghĩa là ta sẽ cần xuất cả $10 \cdot 1000 = 10{,}000$ số (vì mỗi giá trị $N$ sẽ xuất 10 số), thì ta hãy nên để 2 dòng này vào để chạy nhanh hơn. Tất nhiên nếu các bạn code theo cú pháp của C thì chẳng cần mấy dòng đó làm gì, tuy nhiên thường khi code các bài lập trình thi đấu thế này ta dùng C++ vì cú pháp ngắn gọn cũng như nhiều hỗ trợ của nó, vậy nên bạn nào code C++ thì nhớ thêm 2 dòng này vào đầu hàm main nhé.
  2. while(t--) là viết gọn của "kiểm tra t != 0 rồi giảm t", tức là vòng lặp chạy đủ $t$ lần. Trong lập trình khi gọi 1 biến ra trong câu lệnh điều kiện mà không có so sánh gì, các bạn hiểu là đang so sánh khác 0 nhé, tức là: if(x) ⇔ if(x != 0). Còn if(!x) ⇔ if(x == 0).
  3. int a[10] = {0}; nghĩa là khởi tạo mảng a 10 phần tử đều mang giá trị 0 hết. Nói chính xác hơn thì bản chất khi ta khởi tạo mảng = {0}; thì nghĩa là gán giá trị 0 cho ô đầu tiên của mảng (a[0]), và nguyên tắc là chỉ cần có danh sách khởi tạo thì toàn bộ các ô còn lại không được ghi rõ sẽ mặc định giá trị 0. Nên nếu ô đầu tiên cũng được khởi tạo là 0 thì toàn bộ mảng đều có giá trị 0 hết nha các bạn.
  4. Gán x = i: với mỗi giá trị $i$ trong vòng lặp for thì ta phải gán nó qua biến $x$ rồi từ đó đi phân tích biến $x$, chứ ta không làm trực tiếp trên biến $i$, bởi vì thao tác chia nguyên cho 10 sẽ làm thay đổi giá trị biến $i$. Khi đó vòng for sẽ bị sai lệch: $i$ bị chia về 0, rồi ++i thành 1, cứ thế $i$ luôn quanh quẩn ở giá trị 1 và vòng lặp for sẽ bị lặp vô tận.
  5. Khởi tạo mảng bên trong vòng lặp T: bước int a[10] = {0}; ta phải để bên trong vòng lặp $T$ để với mỗi bộ test đọc vào sẽ reset lại mảng toàn giá trị 0, chứ nếu ta mà để nó bên trên vòng lặp $T$ thì các giá trị chữ số cứ thế sẽ bị cộng dồn liên tục sẽ bị sai nha các bạn.

Đánh giá độ phức tạp

Không gian$O(1)$ – mảng 10 phần tử
Thời gian$O(T \cdot N \cdot \log N)$ → TLE

Vòng lặp while chạy $T$ lần để đọc $T$ bộ test, với mỗi bộ test đọc vào giá trị $N$ thì sẽ khởi tạo mảng 10 phần tử 0, sau đó có vòng lặp for $i$ chạy từ 1 tới $N$, với mỗi giá trị $i$ lại chạy đủ bấy nhiêu lần số lượng chữ số của $i$ để phân tích ra. Sau cùng thì chạy vòng lặp 10 lần để in ra kết quả của test case đó. Vậy nên độ phức tạp sẽ là $O(T \cdot (10 + N \cdot \log X + 10))$, với $X$ là các số từ 1 đến $N$. Theo quy tắc cộng Big O thì giữa các thành phần cộng với nhau sẽ lấy theo thành phần lớn nhất nên ta có thể kết luận độ phức tạp là $O(T \cdot N \cdot \log X)$.

Trường hợp xấu nhất: TLE
$T = 10^3$, mọi giá trị $N$ đọc vào đều là lớn nhất $N = 10^6$, $\log X$ cứ tính luôn là 6 cho tất cả. Số lần lặp sẽ là $10^3 \cdot 10^6 \cdot 6 = 6 \cdot 10^9$, tức là 6 tỷ lần lặp, vượt ngưỡng tối đa được phép (khoảng $10^8$) cả 60 lần, tức là phải cần cả 60 giây (1 phút) thì nó mới chạy in ra kết quả, trong khi đó đề yêu cầu chỉ được 1 giây thôi. Nên cách làm ở trên là toang rồi nha các bạn.

Tuy nhiên nếu bạn nào chưa nghĩ ra cách làm tối ưu thì cứ cài đặt theo cách bất kỳ, các bạn có thể giống như cách này chẳng hạn. Thậm chí theo mình là nên cài đặt theo cách này, để tuy nó chạy chậm nhưng nó luôn đảm bảo ra cho mình kết quả chính xác, từ đó ta có thể lấy nó làm công cụ kiểm tra tính chính xác cho một cách làm tối ưu khác tiếp theo (vét cạn + sinh test ngẫu nhiên rồi so sánh 2 cách).

Ý tưởngĐừng tính lại: a[N] = a[N − 1] + các chữ số của N

Lưu lại kết quả đã tính?

Giờ ta cần tìm ra một cách làm tốt hơn cho bài này. Ta cần nhìn thấy cách làm ở trên có cái dở là các giá trị $N$ nếu bị trùng nhau thì ta vẫn cứ phải đi tính lại. Tức là nếu giả sử trước đó ta đã xử lý tính ra kết quả của trường hợp $N = 10^6$ rồi, nếu bộ test tiếp theo đọc vào lại gặp $N = 10^6$ thì cách làm ở trên vẫn cứ đi xử lý lại, trong khi đó rõ ràng ta không cần phải xử lý lại vì ta đã tính ra trước đó kết quả của trường hợp này rồi mà. Vậy nên rõ ràng một cách để tốt hơn đó là ta có thể lưu trữ lại các kết quả $N$ đã tính, để mỗi giá trị $N$ mới ta có thể kiểm tra ngay trong $O(1)$ (nhờ vào kỹ thuật mảng đánh dấu) xem kết quả của trường hợp $N$ đó mình đã từng tính chưa. Nếu tính rồi thì ta lấy ra kết quả in ra thôi, khỏi phải đi xử lý.

Tuy nhiên việc tối ưu ở trên cũng chưa thực sự giải quyết được vấn đề. Giả sử $T$ bộ test toàn bộ đều là các giá trị phân biệt khác nhau thì với mỗi giá trị ta vẫn đều phải đi tính như bình thường. Ví dụ khi $T = 1000$ và các giá trị $N$ lần lượt từ $10^6$ trở về $999{,}001$ ($N = 1{,}000{,}000$, $N = 999{,}999$, $N = 999{,}998$, …) thì mọi giá trị $N$ ta đều phải đi tính chứ không có kết quả sẵn. Lúc này cứ căn cứ $N$ theo giá trị nhỏ nhất thì số lần lặp vẫn cỡ $1000 \cdot 999{,}000 \cdot 6 \approx 6$ tỷ lần lặp. Như thế cũng vẫn sẽ bị TLE thôi vì đã quá 100 triệu.

Mấu chốt: kết quả của N nằm sẵn trong kết quả của N − 1

Ta phải nhìn ra mấu chốt như sau: có phải nếu ta đã có kết quả thống kê số lượng các chữ số từ 0 đến 9 của trường hợp $N = 100$ rồi, thì qua trường hợp $N = 101$ tính ra nó vẫn sẽ là kết quả của trường hợp $N = 100$, chỉ đi phân tích thêm số 101 được thêm 1 chữ số 0 và 2 chữ số 1, lấy nó cộng thêm vào kết quả đã có của $N = 100$ là ta có kết quả của trường hợp $N = 101$ đúng không các bạn?

0123456789
a[100]11212020202020202020
+ số 101+1+2
a[101]12232020202020202020

Tức là nếu trường hợp $N = 100$ ta đã có dãy kết quả là 11 21 20 20 20 20 20 20 20 20 thì với trường hợp $N = 101$ ta sẽ cộng thêm vào dãy đó ở chữ số 0 thêm 1 đơn vị và ở chữ số 1 thêm 2 đơn vị, ta sẽ có dãy kết quả mới của trường hợp $N = 101$ là 12 23 20 20 20 20 20 20 20 20. Tại sao lại như vậy? Vì với $N = 100$ nghĩa là thống kê các chữ số từ 1 đến 100, với $N = 101$ nghĩa là thống kê các chữ số từ 1 đến 101, thì tính ra trường hợp $N = 101$ chỉ là có thêm số 101, còn lại từ số 1 đến 100 thì y như trường hợp $N = 100$.

Công thức chuyển
$a[N] = a[N - 1] + (\text{các chữ số của } N)$, với $a[0]$ là dãy 10 số 0. Tính lần lượt $a[1], a[2], \ldots, a[10^6]$ một lần duy nhất, sau đó mỗi bộ test chỉ việc in ra $a[N]$.

Ở trên là điểm mấu chốt, khi ta đã nhìn ra được rồi thì sẽ dễ dàng nhìn ra tư tưởng giải quyết cho bài này. Để tính ra trường hợp $N$ thì ta cần có kết quả từ trường hợp $N - 1$. Ta tiếp tục thấy $N$ lớn nhất là $10^6$ và nhỏ nhất là 1 theo như giới hạn đề bài. Vậy có phải ta đi tính ra các kết quả của các trường hợp từ 1 đến $10^6$ thì ta chỉ chạy vòng lặp đủ $10^6$ lần là có đủ kết quả các trường hợp $N$ từ 1 đến $10^6$ đúng không các bạn? Bước xử lý phân tách các chữ số của $N$ cứ cho rằng lớn nhất là có 6 chữ số thì bước này lặp 6 lần, và mỗi kết quả là dãy 10 số sẽ được gán qua cho trường hợp mới, nên cả thảy ta chỉ tính $10^6 \cdot (10 + 6) = 1.6 \cdot 10^7$ lần lặp và lưu trữ hết vào trong mảng 2 chiều a có $10^6$ dòng, 10 cột. Rồi sau đó chạy vòng lặp $T$ bộ test, với mỗi bộ test đọc vào số $N$ thì ta gọi a[N] chính là dãy 10 kết quả của trường hợp $N$.

Như thế thì số lần lặp cả thảy của bài chỉ là $1.6 \cdot 10^7 + T \cdot 10$, với trường hợp xấu nhất $T = 1000$ thì số lần lặp sẽ là $16{,}010{,}000$ (16 triệu + 10 ngàn) lần lặp. Như thế này sẽ yên tâm không bị TLE vì nó nhỏ hơn ngưỡng giới hạn 100 triệu, cụ thể ít hơn cả 6 lần. Việc tạo mảng 2 chiều a có $10^6$ dòng, 10 cột là cách dùng kỹ thuật mảng đánh dấu để sau đó mỗi khi gặp giá trị $N$ ta chỉ cần gọi a[N] trong $O(1)$ là ra kết quả ngay lập tức.

Mô phỏng ý tưởng

Em nhập $N$ (tối đa 60) rồi bấm Tiếp để xem từng dòng $a[i]$ được chép từ $a[i - 1]$ rồi cộng thêm các chữ số của $i$:

Ví dụ nhanh:

Cách 2Tính trước a[1..10⁶] – AC

Source code để các bạn tham khảo với hướng làm mình trình bày ở trên:

#include <bits/stdc++.h>
using namespace std;

void PrintArray(const vector<int> &x){
    for(int i = 0; i < (int)x.size(); ++i){
        cout << x[i] << " ";
    }
    cout << endl;
}

int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    vector<vector<int>> a(1000001, vector<int>(10, 0));
    for(int i = 1; i <= 1000000; ++i){
        a[i] = a[i - 1];
        int x = i;
        while(x != 0){
            a[i][x % 10]++;
            x /= 10;
        }
    }
    int t;
    cin >> t;
    while(t--){
        int n;
        cin >> n;
        PrintArray(a[n]);
    }
    return 0;
}
import sys
input = sys.stdin.readline

def PrintArray(x):
    print(*x)

# precompute
a = [[0] * 10 for _ in range(1_000_001)]
for i in range(1, 1_000_001):
    # copy previous row
    a[i] = a[i - 1].copy()
    x = i
    while x != 0:
        a[i][x % 10] += 1
        x //= 10

t = int(input())
for _ in range(t):
    n = int(input())
    PrintArray(a[n])
C++▶ Chạy tay từng dòng codeCách 2: tính trước a[i] = a[i − 1] + các chữ số của i cho mọi i, mỗi bộ test chỉ in a[N].
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3
4void PrintArray(const vector<int> &x){
5 for(int i = 0; i < (int)x.size(); ++i){
6 cout << x[i] << " ";
7 }
8 cout << endl;
9}
10
11int main(){
12 ios::sync_with_stdio(0);
13 cin.tie(0);
14 vector<vector<int>> a(1000001, vector<int>(10, 0));
15 for(int i = 1; i <= 1000000; ++i){
16 a[i] = a[i - 1];
17 int x = i;
18 while(x != 0){
19 a[i][x % 10]++;
20 x /= 10;
21 }
22 }
23 int t;
24 cin >> t;
25 while(t--){
26 int n;
27 cin >> n;
28 PrintArray(a[n]);
29 }
30 return 0;
31}
Lưu ý về bản Python
Code C++ thì AC nhưng code Python ở trên nộp sẽ bị TLE/MLE. Lý do Python khả năng chạy không được nhanh như C++ nên việc lặp qua cả chục triệu lần là không khả thi trong 1 giây; thêm nữa 1 triệu list Python, mỗi list 10 số tốn rất nhiều bộ nhớ (mình đo được khoảng 1.8 giây và gần 300 MB).

Giải thích một số chỗ trong code

  1. vector<vector<int>> a(1000001, vector<int>(10, 0)); nghĩa là khởi tạo vector 2 chiều a có 1,000,001 dòng và mỗi dòng có 10 cột mang giá trị mặc định là 0. Lý do phải khởi tạo 1,000,001 dòng là để index chạy từ 0 đến 1,000,000 (do index tính từ 0), coi như ta bỏ qua không dùng dòng index 0, chỉ dùng từ 1 đến 1 triệu. Ta làm thế vì đề cho $N$ tính từ 1, nên để thống nhất tránh cứ phải rắc rối tính $N - 1$, tốn thêm 10 ô nhớ mà không bị rối rắm có thể dẫn đến sai. (Dòng 0 toàn số 0 cũng chính là điểm xuất phát $a[0]$ cho công thức a[i] = a[i - 1].)
  2. Vì sao dùng vector thay vì mảng thường? Có nhiều lý do. Đầu tiên mảng cần có đến cả 10 triệu phần tử thì không thể khai báo int a[1000001][10] bên trong hàm main được vì sẽ bị tràn bộ nhớ stack, phải dùng con trỏ cấp 2 để cấp phát mảng động, hoặc nếu không thì đưa mảng ra ngoài hàm main làm biến toàn cục thì mới được. Còn nếu trong hàm main mà để int a[1000001][10] các bạn thử chạy trên máy đi sẽ thấy bị tắt chương trình luôn nhé. Tiếp theo đó do ta cần liên tục copy dãy 10 số kết quả của trường hợp $N - 1$ sang trường hợp $N$, nếu dùng vector sẽ rất tiện, chỉ cần làm phép gán a[N] = a[N - 1] thì nó sẽ tự lấy từng phần tử của mảng này gán qua từng phần tử của mảng kia. Nếu dùng mảng thường thì ta phải làm vòng lặp chạy 10 lần để gán qua. Nói chung là gọn hơn được chút xíu ở khâu cài đặt ^^
  3. Hàm PrintArray(const vector<int> &x): để cho gọn thì việc in ra dãy 10 số kết quả ta viết hàm PrintArray nhận tham số đầu vào là vector 1 chiều rồi in ra giá trị các phần tử. Bình thường lẽ ra chỉ cần là vector<int> x thôi đúng không các bạn? Nhưng tại sao lại có thêm const &? Khi truyền tham chiếu (có &) thì sẽ không phải tạo bản sao dữ liệu; ở đây dữ liệu là mảng 10 phần tử thì có & sẽ không phải sao chép 10 phần tử mỗi lần gọi, mà ta có $T$ tối đa 1000 lần nên sẽ tiết kiệm cả 10,000 câu lệnh gán. Còn có thêm const là vì nếu để tham chiếu &, lỡ ta vô tình thay đổi giá trị của nó bên trong hàm (trong khi hàm chỉ cần xuất ra giá trị đâu cần thay đổi gì) thì nó bị thay đổi lại luôn dưới hàm main. Từ khoá const ràng buộc để nếu ta vô tình thay đổi giá trị thì trình biên dịch sẽ báo lỗi ngay, vì const xem đó là hằng không thể thay đổi giá trị được.

Đánh giá độ phức tạp

Không gian$O(10 \cdot 10^6)$
Thời gian$O(10^6 \cdot (10 + 6) + 10 \cdot T)$

Không gian: bộ nhớ cố định 10 triệu phần tử kiểu int → 40 triệu byte (vì int là 4 byte) $/ 1024 \approx 39{,}062.5$ KB $/ 1024 \approx 38.14$ MB. Trong khi đó đề cho giới hạn 128 MB nên ngưỡng ta xài ít hơn cả 3 lần, yên tâm ta sẽ không bị lỗi MLE (Memory Limit Exceeded: quá giới hạn bộ nhớ). (Thực tế mỗi dòng vector còn tốn thêm phần quản lý nên mình đo được khoảng 71 MB, vẫn dư sức.)

Thời gian: như đã phân tích ở cuối phần ý tưởng, khoảng $1.6 \cdot 10^7 + 10 \cdot T$ phép tính.

Ghi chú kiểm tra của mình: cả 2 code C++ và bản Python cách 1 được so với code vét cạn (ghép chuỗi "123…N" rồi đếm từng chữ số) trên 200 bộ test ngẫu nhiên – khớp hết; bản Python cách 2 được so trên 1 bộ 300 test có $N$ tới $10^6$ – khớp. Với test lớn nhất ($T = 1000$, $N$ từ $999{,}001$ đến $10^6$): cách 1 (C++) chạy khoảng 9.6 giây → TLE đúng như phân tích; cách 2 (C++) chạy khoảng 0.1 giây, dùng khoảng 71 MB; cách 2 (Python) khoảng 1.8 giây và gần 300 MB.
Rất vui khi anh em nào đã đọc đến những dòng cuối cùng này. Anh em đọc chắc cũng mỏi mắt luôn chứ 😅 vậy anh em nghĩ người viết ra hết đống ở trên còn thế nào nữa 🤪. Vậy nên đừng tiếc gì cho mình một lời bình luận hoặc chia sẻ bài viết để động viên mình nè. Chúc anh em học tốt và gặt hái nhiều thành công 😇

Cộng một kiểu Đức

Cộng một kiểu Đức
Nguồn: Codeforces Group Độ khó: Trung bình Mảng đếm Quy hoạch động Tiền xử lý

Đề bài

Đức vừa phát minh ra một phép toán cộng độc đáo và đặt tên cho nó là “cộng một kiểu Đức”, cách cộng như sau: Mỗi lần cộng một số nguyên cho một thay vì tăng giá trị của số đó lên một đơn vị, Đức lại tăng tất cả các chữ số của số đó lên một đơn vị.

Ví dụ 1: $15$ “cộng một kiểu Đức” sẽ cho kết quả là $26$.

Ví dụ 2: $91$ “cộng một kiểu Đức” sẽ cho kết quả là $102$.

Yêu cầu: Với $T$ số nguyên dương $N$ cho trước hãy tính toán và cho biết sau khi thực hiện $K$ lần phép cộng một kiểu Đức thì kết quả sinh ra có bao nhiêu chữ số. Vì câu trả lời có thể rất lớn do đó hãy in ra kết quả sau khi chia lấy dư cho $10^9 + 7$.

Dữ liệu vào:

  • Dòng thứ nhất: chứa số nguyên $T$ $(1 \le T \le 2 \cdot 10^5)$.
  • $T$ dòng tiếp theo mỗi dòng chứa 2 số nguyên lần lượt là $N$ và $K$ $(1 \le N \le 10^9;\ 1 \le K \le 2 \cdot 10^5)$.

Kết quả: Với mỗi trường hợp in ra số lượng chữ số của kết quả chia lấy dư cho $10^9 + 7$.

Scoring

  • Có $30\%$ test với $T \le 10,\ K \le 20$;
  • Có $20\%$ test với $T \le 100,\ K \le 50$;
  • Có $50\%$ test không có giới hạn gì thêm.
Ví dụ – input
2
9812 2
9899 3
Ví dụ – output
6
8

Note:

  • $9812 \to 10923 \to 211034$ (6 chữ số)
  • $9899 \to 1091010 \to 21102121 \to 32213232$ (8 chữ số)

Lời mở đầu

Bài này mình được một bạn học viên hỏi nên sẵn tiện mình viết bài giải chi tiết để nhiều bạn được tham khảo luôn. Hy vọng mang đến được giá trị cho các bạn 😇.

Dành cho học viên · Ưu tiênNộp bài trên hệ thống CLUE OJ

Link này chỉ học viên đã tham gia tổ chức mới vào được. Học viên nên ưu tiên nộp ở đây vì sẽ được chấm nhanh hơn.

Nộp bài trên CLUE OJ
Dành cho tất cả mọi ngườiNộp bài trên Codeforces

Bạn nào muốn thử nộp bài để hệ thống chấm xem có đúng hay không thì dùng link này nhé. Chưa có tài khoản thì tạo tài khoản rồi đăng nhập nộp bài thử hen.

Nộp bài trên Codeforces

Lộ trình lời giải, bấm vào từng phần để nhảy tới:

Cách 1Vét cạn bằng chuỗi – 50% test

Theo quy luật như đề bài đưa ra thì khá dễ dàng để làm được code vét cạn em ha, và ít nhất nó sẽ lấy được 2 subtask đầu tiên tức là $50\%$ điểm. Cũng như ta có thể lợi dụng nó làm tool test cho cách tối ưu (sinh test ngẫu nhiên nhỏ, so kết quả 2 code với nhau).

Ý tưởng code vét cạn thì đơn giản là làm theo đúng như mô tả yêu cầu đề thôi: cho vòng lặp chạy $T$ lần, mỗi lần đọc dữ liệu số $n$ và $k$ thì đối với số $n$ ta nên đọc vào ở dạng chuỗi để từ đó việc xét từng ký tự số sẽ dễ dàng, chứ nếu dùng kiểu số thì sẽ phải phức tạp vấn đề lên. Khi xét từng ký tự số:

  • Nếu nó không phải ký tự '9' thì ta cứ ++ lên thôi.
  • Còn nếu là ký tự '9' thì sẽ biến nó lại thành ký tự '1' và theo lẽ thì phải có ký tự '0' được thêm vào ngay sau ký tự '1' đó.

Nhưng bản chất nếu làm vậy thì sẽ tăng độ phức tạp lên, vì chi phí thêm 1 phần tử vào giữa chuỗi (bản chất là mảng) sẽ là $O(\text{Length})$. Trong khi đó ta thấy rõ ràng dù 9 biến thành 10 rồi thì khi cộng tiếp các lần sau thì nó cũng cộng riêng lẻ cho ký tự 1 rồi riêng cho ký tự 0, chúng không bị phụ thuộc vào nhau. Vì thế ký tự 0 không nhất thiết phải nằm ngay sau ký tự 1, mà ta nên thêm nó vào sau cùng của chuỗi, như thế chỉ mất chi phí $O(1)$ sẽ tốt hơn. Rồi thì lặp lại quy trình cộng liên tục $k$ lần và sau đó thì in kết quả độ dài của chuỗi, nhớ chia dư cho $10^9 + 7$ như đề bài yêu cầu.

kSố thậtChuỗi trong codeDài
09899
9899
4
11091010
1911000
7
221102121
21221110
8
332213232
32332221
8

Ô cam = chữ số 0 mới được push_back vào cuối. Thứ tự chữ số khác số thật nhưng số lượng chữ số luôn bằng nhau – mà đề chỉ hỏi số lượng.

Source code sẽ như sau em tham khảo nhé:

#include <bits/stdc++.h>
using namespace std;
const int m = 1e9 + 7;
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    int t;
    cin >> t;
    while(t--){
        string s;
        int k;
        cin >> s >> k;
        while(k--){
            int length = s.length();
            for(int i = 0; i < length; ++i){
                if(s[i] < '9'){
                    s[i]++;
                }
                else{
                    s[i] = '1';
                    s.push_back('0');
                }
            }
        }
        cout << s.length() % m << '\n';
    }
    return 0;
}
import sys

m = 10**9 + 7
data = sys.stdin.read().split()
t = int(data[0])
out = []
for q in range(t):
    s = list(data[1 + 2 * q])
    k = int(data[2 + 2 * q])
    for _ in range(k):
        length = len(s)
        for i in range(length):
            if s[i] < '9':
                s[i] = chr(ord(s[i]) + 1)
            else:
                s[i] = '1'
                s.append('0')  # thêm '0' vào cuối cho nhanh, thứ tự không ảnh hưởng số chữ số
    out.append(len(s) % m)
print('\n'.join(map(str, out)))
C++▶ Chạy tay từng dòng codeCách 1 (vét cạn): đọc N dạng chuỗi, mỗi lần cộng duyệt từng ký tự, chữ số 9 thành 1 và đẩy 0 vào cuối.
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3const int m = 1e9 + 7;
4int main(){
5 ios::sync_with_stdio(0);
6 cin.tie(0);
7 int t;
8 cin >> t;
9 while(t--){
10 string s;
11 int k;
12 cin >> s >> k;
13 while(k--){
14 int length = s.length();
15 for(int i = 0; i < length; ++i){
16 if(s[i] < '9'){
17 s[i]++;
18 }
19 else{
20 s[i] = '1';
21 s.push_back('0');
22 }
23 }
24 }
25 cout << s.length() % m << '\n';
26 }
27 return 0;
28}

Để ý biến length được lấy trước vòng for: các chữ số 0 vừa được push_back trong lần cộng này sẽ không bị cộng thêm ngay trong cùng lần đó.

Độ phức tạp thời gian sẽ là $O(T \cdot K \cdot \text{Length})$, và rõ ràng nội việc chỉ cần $T \cdot K$ thì bài đó $T, K$ tối đa là $2 \cdot 10^5$, như vậy $T \cdot K = 4 \cdot 10^{10}$ tức là 40 tỷ, đã vượt quá ngưỡng giới hạn thời gian chạy cho 1 giây (tối đa chỉ khoảng 100 triệu phép tính). Tuy nhiên cái đáng sợ hơn là cái $\text{Length}$, vì nó sẽ là 1 con số rất khủng khiếp nếu rơi vào trường hợp xấu. Để anh sẽ cho em thấy:

Chuỗi dài khủng khiếp cỡ nào?

Nếu như dùng chuỗi là dở vì chuỗi sinh ra có thể sẽ rất dài, thì sợ rằng memory lưu không nổi luôn đó chứ. Vì em ngẫm thử 1 trường hợp xấu mà anh giả sử này nhé:

  • Giả sử $n$ ban đầu $= 999999999$ (9 chữ số 9), vậy thì sau 1 lần cộng sẽ tạo ra 18 chữ số: 101010101010101010.
  • Tiếp tục sau 9 lần cộng nữa thì 9 con số 1 sẽ tạo ra thêm 9 con số mới.
  • Rồi thêm 1 lần cộng nữa thì 9 con số 0 sẽ tạo ra thêm 9 con số mới.
  • Nên mình có thể hiểu là: cứ sau 1 (lần ban đầu) + 9 + 1 = 11 lần thì sẽ gấp đôi thành 36 chữ số; rồi sau 11 lần nữa thì lại gấp đôi thành 72 chữ số; rồi 144 chữ số…

Nếu không tin thì em thử chạy đoạn code vét cạn ở trên với các test 1 / 999999999 11 sẽ thấy kết quả là 36, 999999999 22 ra 72, 999999999 33 ra 144, 999999999 44 ra 288. Nó cứ thế nhân 2 lên liên tục em hình dung được không?

k = 118 = 18 · 20
k = 1136 = 18 · 21
k = 2272 = 18 · 22
k = 33144 = 18 · 23
k = 44288 = 18 · 24
Số chữ số của 999999999 sau k lần "cộng một kiểu Đức": cứ mỗi 11 lần là gấp đôi.

Mà em nhìn xem: $K$ đề cho tối đa là $2 \cdot 10^5$, vậy sẽ có $K / 11 \approx 18181$ lần gấp đôi (18 ngàn 181 lần). Nghĩa là kết quả sau cùng chuỗi sẽ có độ dài là khoảng $18 \cdot 2^{18181}$, quá lớn rồi đó em. Lũy thừa $2^{30}$ đã là hơn 1 tỷ rồi; kiểu số nguyên lớn nhất mà C++ hỗ trợ sẵn là kiểu 8 byte (64 bit) thì tối đa nó lưu được cỡ $2^{64}$, trong khi ở đây là cần đến số mũ 18 ngàn 181. Chưa kể là chuỗi sinh ra lại phải duyệt vòng lặp vào từng ký tự để ++ lên. Đó, nếu dùng chuỗi là thua: memory lưu không nổi, time chạy cũng không nổi.

Cách 2Mảng đếm 10 chữ số – O(T·K·10)

Mấu chốt em có độ phức tạp cao như vậy là bởi vì em sử dụng chuỗi. Em phải nghĩ theo hướng không dùng chuỗi thì mới tối ưu nhé: thực tế chỉ cần dùng mảng đánh dấu có 10 phần tử tương ứng từ 0 đến 9. Ví dụ mảng a thì a[i] sẽ lưu số lượng chữ số $i$, ví dụ a[3] sẽ lưu số lượng chữ số 3. Nếu số tạo ra có 100 chữ số 1, thay vì em tạo ra chuỗi "1111...111" (100 số 1) thì em chỉ cần thống kê tại index 1 mang giá trị 100, tức là a[1] = 100. Làm như thế thì em sẽ không phải lưu trữ quá nhiều và phải duyệt quá nhiều.

Sau 1 lần cộng, số lượng chữ số $d$ ($d < 9$) chuyển hết sang thành chữ số $d + 1$; còn số lượng chữ số 9 thì mỗi chữ số 9 biến thành 10, tức là góp vào cả số lượng chữ số 1 lẫn chữ số 0:

0123456789b[d + 1] = a[d] (d = 0..8)b[0] += a[9], b[1] += a[9]

Rồi sau cùng thì khi duyệt số lượng các chữ số để cộng tổng lại thì đó là đáp án. Thay vì cộng tổng lại hết rồi ta mới chia lấy dư thì ta sẽ áp dụng đồng dư thức $(a + b) \bmod c$ $= ((a \bmod c) + (b \bmod c)) \bmod c$, từ đó cứ mỗi lần cộng vào là ta chia lấy dư luôn để tránh trường hợp bị tràn số do quá lớn.

Em nhập $N$, $K$ rồi bấm Tiếp để xem cùng lúc số thật và mảng đếm biến đổi qua từng lần cộng:

Ví dụ nhanh:

Source code sẽ như sau:

#include <bits/stdc++.h>
using namespace std;
const int m = 1e9 + 7;
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    int t;
    cin >> t;
    while(t--){
        int n, k;
        cin >> n >> k;
        vector<int> a(10, 0);
        while(n != 0){
            a[n % 10]++;
            n /= 10;
        }
        while(k--){
            vector<int> b(10, 0);
            for(int i = 0; i < 9; ++i){
                b[i + 1] = a[i];
            }
            b[1] += a[9];
            b[0] += a[9];
            a = b;
        }
        int res = 0;
        for(int i = 0; i < 10; ++i){
            res += a[i] % m;
            res %= m;
        }
        cout << res << '\n';
    }
    return 0;
}
import sys

m = 10**9 + 7
data = sys.stdin.read().split()
t = int(data[0])
out = []
for q in range(t):
    n = int(data[1 + 2 * q])
    k = int(data[2 + 2 * q])
    a = [0] * 10          # a[d] = số lượng chữ số d
    while n != 0:
        a[n % 10] += 1
        n //= 10
    for _ in range(k):
        b = [0] * 10
        for i in range(9):
            b[i + 1] = a[i]
        b[1] += a[9]
        b[0] += a[9]
        a = b
    out.append(sum(a) % m)
print('\n'.join(map(str, out)))
C++▶ Chạy tay từng dòng codeCách 2: chỉ giữ mảng đếm 10 ô, mỗi lần cộng dịch số lượng sang phải, số lượng chữ số 9 cộng vào ô 0 và 1.
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3const int m = 1e9 + 7;
4int main(){
5 ios::sync_with_stdio(0);
6 cin.tie(0);
7 int t;
8 cin >> t;
9 while(t--){
10 int n, k;
11 cin >> n >> k;
12 vector<int> a(10, 0);
13 while(n != 0){
14 a[n % 10]++;
15 n /= 10;
16 }
17 while(k--){
18 vector<int> b(10, 0);
19 for(int i = 0; i < 9; ++i){
20 b[i + 1] = a[i];
21 }
22 b[1] += a[9];
23 b[0] += a[9];
24 a = b;
25 }
26 int res = 0;
27 for(int i = 0; i < 10; ++i){
28 res += a[i] % m;
29 res %= m;
30 }
31 cout << res << '\n';
32 }
33 return 0;
34}

Độ phức tạp thời gian giờ đây chỉ là $O(T \cdot K \cdot 10)$. Rõ ràng nó tốt hơn rất nhiều so với cách làm trước đó khi độ phức tạp là $O(T \cdot K \cdot \text{Length})$ với $\text{Length}$ là 1 con số rất khủng khiếp, thì giờ đây $\text{Length}$ đã thay bằng số 10. Tuy nhiên thì như đã nói $T, K$ tối đa là $2 \cdot 10^5$ nên $T \cdot K = 4 \cdot 10^{10}$, giờ tiếp tục nhân thêm 10 nữa thì thành 400 tỷ. Trong khi đó ngưỡng giới hạn 1 giây là 100 triệu, vì thế chắc chắn code này tuy tối ưu hơn nhưng vẫn sẽ bị TLE. Dù sao thì cũng rất đáng khen ngợi khi từ code trước đó ta tối ưu ra code này, vì trường hợp xấu nhất anh đề ra đó, nếu như code trước đó chắc là đến khi trái đất tận thế cũng chưa ra kết quả, còn với code này thì chỉ mất $400$ tỷ $/ 100$ triệu $= 4000$ giây, hơn 1 giờ là ra kết quả rồi 😄.

Lưu ý tràn số (mình ghi thêm)
Trong code trên các ô a[i], b[i] là int và không được chia dư trong lúc lặp, nên khi $K$ lớn (số chữ số tăng theo cấp số nhân, cỡ $2^{K/11}$) thì sẽ bị tràn số. Ví dụ test 1 / 999999999 200000 code này in ra số âm -212862980 thay vì 831614621. Ở 2 subtask đầu ($K \le 50$) thì không sao. Muốn an toàn em chỉ cần sửa 2 dòng thành b[1] = (b[1] + a[9]) % m; và b[0] = (b[0] + a[9]) % m; (giống như cách mod trong Cách 3).

Nhưng ta sẽ cần tìm ra cách tối ưu để đảm bảo ra được kết quả trong tối đa 1 giây. Hãy tiếp tục xem nhé:

Cách 3Bảng thống kê theo từng chữ số và từng K – AC

Mấu chốt của bài này là nhận ra như sau: ví dụ số $n$ có dạng abcde và cần đếm xem sau $k$ lần cộng sẽ có bao nhiêu chữ số theo quy tắc đề bài đưa ra. Thì sẽ hiểu là:

  • Đếm xem chữ số a sau $k$ lần cộng sẽ tạo ra bao nhiêu chữ số;
  • Đếm xem chữ số b sau $k$ lần cộng sẽ tạo ra bao nhiêu chữ số;
  • … tương tự với c, d, e.

Cộng tổng lại thì đó là đáp án (vì các chữ số biến đổi độc lập với nhau, như đã thấy ở Cách 1).

Vấn đề là có $T$ test case, mỗi test case thì $N$ trong đoạn $[1, 10^9]$, $K$ trong đoạn $[1, 2 \cdot 10^5]$, thì nhận thấy mấu chốt chỉ là xét các chữ số của số $N$ thôi, và ta biết chỉ có 10 chữ số từ 0 đến 9, $K$ thì cũng chỉ trong đoạn từ 1 đến $2 \cdot 10^5$. Nên là ta tạo bảng thống kê:

  • Ứng với từng chữ số (từ 0 đến 9),
  • Ứng với từng giá trị $k$ (từ 1 đến $2 \cdot 10^5$),
  • Thì xem coi chữ số đó, sau $k$ lần cộng, sẽ tạo ra bao nhiêu chữ số?

Khi có được bảng thống kê rồi thì với từng test case, với số $n$ ta lấy ra từng chữ số của nó, gọi vào chữ số tương ứng trong bảng thống kê và tại vị trí $k$ tương ứng, ta biết được chữ số đó sau $k$ lần cộng sẽ tạo ra bao nhiêu chữ số mới, rồi cộng tổng lại với toàn bộ các chữ số của $n$ là ra đáp án sau cùng.

Chi phí tạo bảng: ta có 10 chữ số, với 1 chữ số thì xét qua $2 \cdot 10^5$ cái $K$, và áp dụng tư tưởng dùng mảng đánh dấu thay vì chuỗi như ở trên thì mỗi cái $K$ sẽ xét qua mảng đánh dấu có 10 phần tử. Vậy $10 \cdot 2 \cdot 10^5 \cdot 10 = 2 \cdot 10^7$, nhỏ hơn ngưỡng $10^8$ nên vẫn ổn nè.

Ví dụ chữ số 9 đi ha: $9 \to 10 \to 21 \to 32 \to \ldots \to 98 \to 109 \to 2110 \to 3221 \to 4332$. Đó thì em cứ tạo bảng thống kê cho đến $+\,2 \cdot 10^5$. Và không nên lưu theo dạng chuỗi vì như anh đã nói nó sẽ có độ phức tạp cực kỳ lớn cả về thời gian lẫn bộ nhớ, hãy dùng mảng đánh dấu có 10 phần tử để lưu số lượng các chữ số thôi:

k9 sau k lần0123456789Tổng
0900000000011
11011000000002
22101100000002
33200110000002
44300011000002
55400001100002
66500000110002
77600000011002
88700000001102
99800000000112
1010911000000013
11211012100000004
12322101210000004
13433200121000004

Ví dụ số 9: trước khi cộng thì mảng đánh dấu của nó là 0, 0, 0, 0, 0, 0, 0, 0, 0, 1. Sau khi cộng 1 lần thành 10 thì số lượng chữ số 9 mất đi, thay vào đó số lượng chữ số 9 đó được cộng vào vị trí chữ số 0 và 1, ta có mảng đánh dấu là 1, 1, 0, 0, 0, 0, 0, 0, 0, 0. Nhìn dòng cuối cùng: $9$ cộng 13 lần $= 4332$ nên a[9][13] sẽ là mảng đánh dấu có 1 số 2, 2 số 3, 1 số 4, nghĩa là a[9][13] = {0, 0, 1, 2, 1, 0, 0, 0, 0, 0}.

Mảng 3 chiều int a[10][200001][10]
a[i][k][j] = số lượng chữ số $j$ sinh ra từ chữ số $i$ sau $k$ lần cộng.
• Số 10 đầu tiên: các chữ số $i$ từ 0 đến 9.
• Số 200001: $k$ từ 0 đến $2 \cdot 10^5$.
• Số 10 cuối: mảng đánh dấu 10 phần tử (chữ số $j$).
Công thức chuyển y như Cách 2: a[i][k][j] = a[i][k-1][j-1] với $j = 1..9$, rồi a[i][k][0] += a[i][k-1][9], a[i][k][1] += a[i][k-1][9].

Tạo ra được cái bảng này chỉ mất $10 \cdot 2 \cdot 10^5 \cdot 10 = 2 \cdot 10^7$ thôi, vẫn ổn. Rồi khi có nó rồi là xong, khỏe: với mỗi số $n$, lấy từng chữ số, với từng chữ số sẽ biết ngay chữ số ấy sau khi trải qua $K$ lần sẽ tạo ra bao nhiêu chữ số khác, cứ cộng tổng lại là xong.

Tuy nhiên lưu ý em rằng nhớ cái vụ mod $10^9 + 7$ thì nên mod ngay từ lúc tạo mảng 3 chiều và mod thêm lúc cộng tổng các chữ số của $N$. Như vậy bước sau đó chỉ cần xét qua toàn bộ chữ số của số $N$ là ra kết quả. Mà ta biết $N$ tối đa có 10 chữ số ($N \le 10^9$), $T$ tối đa $2 \cdot 10^5$ nên khúc sau cũng chỉ cần lặp $2 \cdot 10^5 \cdot 10 \cdot 10 = 2 \cdot 10^7$ lần là xong.

Source code làm theo đúng tư tưởng anh đã nói đó và đã nộp AC full 100% trên web em xem nhé:

#include <bits/stdc++.h>
using namespace std;
int a[10][200001][10];
const int m = 1e9 + 7;
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    for(int i = 0; i < 10; ++i){
        a[i][0][i] = 1;
    }
    for(int i = 0; i < 10; ++i){
        for(int k = 1; k <= 200000; ++k){
            for(int j = 1; j < 10; ++j){
                a[i][k][j] = a[i][k - 1][j - 1];
            }
            a[i][k][0] += a[i][k - 1][9] % m;
            a[i][k][0] %= m;
            a[i][k][1] += a[i][k - 1][9] % m;
            a[i][k][1] %= m;
        }
    }
    int t;
    cin >> t;
    while(t--){
        int n, k;
        cin >> n >> k;
        int res = 0;
        while(n != 0){
            int x = n % 10;
            n /= 10;
            for(int j = 0; j < 10; ++j){
                res += a[x][k][j];
                res %= m;
            }
        }
        cout << res << '\n';
    }
    return 0;
}
import sys

m = 10**9 + 7
MAXK = 200000
# Bảng 3 chiều 10 x 200001 x 10 quá nặng với Python, nên mình chỉ lưu
# tổng số chữ số: total[i][k] = tổng a[i][k][0..9] (cách tính y như C++)
total = [[0] * (MAXK + 1) for _ in range(10)]
for i in range(10):
    cur = [0] * 10
    cur[i] = 1
    total[i][0] = 1
    for k in range(1, MAXK + 1):
        # b[j] = a[j - 1] (j = 1..9), rồi b[0] += a[9], b[1] += a[9]
        c9 = cur[9]
        cur = [c9, (cur[0] + c9) % m] + cur[1:9]
        total[i][k] = sum(cur) % m

data = sys.stdin.buffer.read().split()
t = int(data[0])
out = []
for q in range(t):
    n = int(data[1 + 2 * q])
    k = int(data[2 + 2 * q])
    res = 0
    while n != 0:
        res = (res + total[n % 10][k]) % m
        n //= 10
    out.append(res)
print('\n'.join(map(str, out)))
C++▶ Chạy tay từng dòng codeCách 3: dựng sẵn bảng a[i][k][j] cho mọi chữ số i và mọi k, mỗi test chỉ cần cộng các ô tương ứng.
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3int a[10][200001][10];
4const int m = 1e9 + 7;
5int main(){
6 ios::sync_with_stdio(0);
7 cin.tie(0);
8 for(int i = 0; i < 10; ++i){
9 a[i][0][i] = 1;
10 }
11 for(int i = 0; i < 10; ++i){
12 for(int k = 1; k <= 200000; ++k){
13 for(int j = 1; j < 10; ++j){
14 a[i][k][j] = a[i][k - 1][j - 1];
15 }
16 a[i][k][0] += a[i][k - 1][9] % m;
17 a[i][k][0] %= m;
18 a[i][k][1] += a[i][k - 1][9] % m;
19 a[i][k][1] %= m;
20 }
21 }
22 int t;
23 cin >> t;
24 while(t--){
25 int n, k;
26 cin >> n >> k;
27 int res = 0;
28 while(n != 0){
29 int x = n % 10;
30 n /= 10;
31 for(int j = 0; j < 10; ++j){
32 res += a[x][k][j];
33 res %= m;
34 }
35 }
36 cout << res << '\n';
37 }
38 return 0;
39}
Không gian$O(10 \cdot K_{max} \cdot 10)$ ≈ 80 MB
Thời gian$O(K_{max} \cdot 100 + T \cdot 100)$

Bảng int a[10][200001][10] chiếm khoảng $2 \cdot 10^7 \cdot 4$ byte $\approx 80$ MB; khai báo toàn cục nên các ô ban đầu đều bằng 0 và không bị tràn stack. Bản Python mình viết thêm không lưu nguyên bảng 3 chiều (quá nặng với Python) mà chỉ lưu tổng số chữ số total[i][k], cách tính y hệt.

Cách 4Thu gọn bảng về 1 chiều (mình viết thêm)

Phần này mình viết thêm cho bạn nào muốn tối ưu tiếp nha 😉. Để ý 2 điều:

  1. Chữ số $d$ chính là chữ số $0$ sau $d$ lần cộng ($0 \to 1 \to \ldots \to d$). Vậy chữ số $d$ sau $k$ lần cộng = chữ số $0$ sau $k + d$ lần cộng. Chỉ cần 1 hàng của bảng thay vì 10 hàng!
  2. Gọi f[k] = số chữ số sinh ra từ chữ số 0 sau $k$ lần cộng. Với $k < 10$ thì chữ số 0 chưa vượt qua 9 nên f[k] = 1. Với $k \ge 10$: sau 10 lần, chữ số 0 thành 10, tức là 1 chữ số 1 và 1 chữ số 0, mỗi đứa còn $k - 10$ lần cộng. Chữ số 1 còn $k - 10$ lần = chữ số 0 còn $k - 9$ lần, nên:
f[k] = f[k − 10] + f[k − 9] (k ≥ 10), f[0..9] = 1

Đáp án của 1 test = tổng f[K + d] với mỗi chữ số $d$ của $N$. Kiểm tra với chữ số 9 cộng 13 lần: f[22] = f[12] + f[13] = 2 + 2 = 4, đúng bằng số chữ số của 4332.

#include <bits/stdc++.h>
using namespace std;
const int m = 1e9 + 7;
const int MAXK = 200000 + 10;
int f[MAXK];
int main(){
    ios::sync_with_stdio(0);
    cin.tie(0);
    // f[k] = số chữ số sinh ra từ chữ số 0 sau k lần cộng
    for(int k = 0; k < MAXK; ++k){
        if(k < 10){
            f[k] = 1;
        }
        else{
            f[k] = (f[k - 10] + f[k - 9]) % m;
        }
    }
    int t;
    cin >> t;
    while(t--){
        int n, k;
        cin >> n >> k;
        int res = 0;
        while(n != 0){
            int x = n % 10;
            n /= 10;
            res = (res + f[k + x]) % m;
        }
        cout << res << '\n';
    }
    return 0;
}
import sys

m = 10**9 + 7
MAXK = 200000 + 10
# f[k] = số chữ số sinh ra từ chữ số 0 sau k lần cộng
f = [1] * MAXK
for k in range(10, MAXK):
    f[k] = (f[k - 10] + f[k - 9]) % m

data = sys.stdin.buffer.read().split()
t = int(data[0])
out = []
for q in range(t):
    n = int(data[1 + 2 * q])
    k = int(data[2 + 2 * q])
    res = 0
    while n != 0:
        res = (res + f[k + n % 10]) % m
        n //= 10
    out.append(res)
print('\n'.join(map(str, out)))
C++▶ Chạy tay từng dòng codeCách 4: f[k] = số chữ số sinh ra từ chữ số 0 sau k lần, f[k] = f[k−10] + f[k−9]; đáp án = tổng f[K + d].
Dữ liệu mẫu:
1#include <bits/stdc++.h>
2using namespace std;
3const int m = 1e9 + 7;
4const int MAXK = 200000 + 10;
5int f[MAXK];
6int main(){
7 ios::sync_with_stdio(0);
8 cin.tie(0);
9 // f[k] = số chữ số sinh ra từ chữ số 0 sau k lần cộng
10 for(int k = 0; k < MAXK; ++k){
11 if(k < 10){
12 f[k] = 1;
13 }
14 else{
15 f[k] = (f[k - 10] + f[k - 9]) % m;
16 }
17 }
18 int t;
19 cin >> t;
20 while(t--){
21 int n, k;
22 cin >> n >> k;
23 int res = 0;
24 while(n != 0){
25 int x = n % 10;
26 n /= 10;
27 res = (res + f[k + x]) % m;
28 }
29 cout << res << '\n';
30 }
31 return 0;
32}
Không gian$O(K_{max})$
Thời gian$O(K_{max} + T \cdot 10)$
Ghi chú của mình: Code C++ của Cách 1, 2, 3 là của anh trong bài gốc; Cách 4 và các bản Python là mình viết thêm. Tất cả 8 bản đã được kiểm tra với 200 bộ test ngẫu nhiên ($T \le 10$, $K \le 50$, có cả $N = 999999999$, $N = 10^9$ chứa nhiều chữ số 0) so với cách cộng trực tiếp trên chuỗi đúng thứ tự chữ số: khớp 100%; thêm 10 bộ $K$ ngẫu nhiên tới $2 \cdot 10^5$ so Cách 3 với Cách 4: khớp. Với test lớn nhất ($T = K = 2 \cdot 10^5$), C++ Cách 3 chạy khoảng 0.2 giây, Cách 4 khoảng 0.05 giây; Python Cách 3 khoảng 2.7 giây (phần lớn là dựng bảng), Python Cách 4 khoảng 0.8 giây. Riêng code C++ Cách 2 bị tràn số khi $K$ lớn như đã lưu ý ở trên.
Rất vui khi anh em nào đã đọc đến những dòng cuối cùng này. Anh em đọc chắc cũng mỏi mắt luôn chứ 😅 vậy anh em nghĩ người viết ra hết đống ở trên còn thế nào nữa 🤪. Vậy nên đừng tiếc gì cho mình một lời bình luận hoặc chia sẻ bài viết để động viên mình nè. Chúc anh em học tốt và gặt hái nhiều thành công 😇