题目描述
给定一个长度为 n n n 的数组 A 1 , A 2 , ⋅ ⋅ ⋅ , A n A_1,A_2,⋅⋅⋅,A_n A1,A2,⋅⋅⋅,An。
你可以从中选出两个数 A i A_i Ai 和 A j A_j Aj( i i i 不等于 j j j),然后将 A i A_i Ai 和 A j A_j Aj 一前一后拼成一个新的整数。
例如 12 和 345 可以拼成 12345 或 34512。
注意交换 A i A_i Ai 和 A j A_j Aj 的顺序总是被视为 2 种拼法,即便是 A i = A j A_i=A_j Ai=Aj 时。
请你计算有多少种拼法满足拼出的整数是 K K K 的倍数。
输入格式
第一行包含 2 个整数 n n n 和 K K K。
第二行包含 n 个整数 A 1 , A 2 , ⋅ ⋅ ⋅ , A n A_1,A_2,⋅⋅⋅,A_n A1,A2,⋅⋅⋅,An。
输出格式
一个整数代表答案。
数据范围
1
≤
n
≤
1
0
5
1≤n≤10^5
1≤n≤105,
1
≤
K
≤
1
0
5
1≤K≤10^5
1≤K≤105,
1
≤
A
i
≤
1
0
9
1≤A_i≤10^9
1≤Ai≤109
输入样例
4 2
1 2 3 4
输出样例
6
算法思想(哈希)
根据题目描述,求 A i A_i Ai 和 A j A_j Aj组成的整数 A i A j ‾ \overline {A_iA_j} AiAj能够整除 K K K的方案总数。进一步分析 A i A j ‾ = A i × 1 0 t + A j \overline {A_iA_j}=A_i\times10^{t}+A_j AiAj=Ai×10t+Aj,其中 t t t表示 A j A_j Aj的位数。那么对于每个 A j A_j Aj来说,要求 A i × 1 0 t + A j ≡ 0 ( m o d K ) A_i\times10^{t}+A_j\equiv0\pmod{K} Ai×10t+Aj≡0(modK)的方案数,即求 A i × 1 0 t ≡ − A j ( m o d K ) A_i\times10^{t}\equiv -A_j\pmod{K} Ai×10t≡−Aj(modK)的方案数。
通过上述分析,可以使用哈希表存储所有 A i × 1 0 t % K A_i\times10^{t} \% K Ai×10t%K的方案数。不妨设 v = A i × 1 0 t % K v = A_i\times10^{t} \% K v=Ai×10t%K,那么 h [ t ] [ v ] h[t][v] h[t][v]表示 A i × 1 0 t % K A_i\times10^{t} \% K Ai×10t%K值为 v v v的方案数。由于 1 ≤ A i ≤ 1 0 9 1≤A_i≤10^9 1≤Ai≤109,所以 0 ≤ t ≤ 10 0\le t \le 10 0≤t≤10;而 1 ≤ K ≤ 1 0 5 1≤K≤10^5 1≤K≤105,所以 0 ≤ v ≤ 10 0\le v \le 10 0≤v≤10。
那么对于每个 A i A_i Ai来说,只需要累加哈希表中 h [ t ] [ v ] h[t][v] h[t][v]的方案数即可,其中 t = A i t = A_i t=Ai的位数, v = ( K − A i % K ) % K v = (K-A_i \% K)\%K v=(K−Ai%K)%K,即 − A [ i ] ( m o d K ) -A[i] \pmod{K} −A[i](modK)的值。
注意,如果 A i × 1 0 t % K A_i\times10^{t} \% K Ai×10t%K的值也等于 v v v,为了保证选出的两个数 A i A_i Ai 和 A j A_j Aj( i i i 不等于 j j j),所以方案总数要减去1。
时间复杂度
需要遍历每个 A i A_i Ai,查询哈希表中 h [ t ] [ v ] h[t][v] h[t][v]的值,所以时间复杂度为 O ( n × t ) = 1 0 6 O(n \times t) = 10^6 O(n×t)=106 。
代码实现
#include <iostream>
#include <cstdio>
using namespace std;
typedef long long LL;
const int N = 100010;
int a[N];
int h[11][N];
int main()
{
int n, k;
scanf("%d%d", &n, &k);
for(int i = 0; i < n; i ++) scanf("%d", &a[i]);
for(int i = 0; i < n; i ++)
{
LL v = a[i] % k; //求余数v
//处理哈希表中h[t][v]的值,即a[i]×10^t % k的方案数
for(int t = 0; t < 11; t ++)
{
h[t][v] ++;
v = v * 10 % k;
}
}
LL res = 0;
for(int i = 0; i < n; i ++)
{
//计算a[i]的长度
int t = to_string(a[i]).size();
//计算−a[i](mod K)的值
LL v = (k - a[i] % k) % k;
//累加哈希表中h[t][v]的方案数
res += h[t][v];
//特殊处理a[i]×10^t%k的值也等于v的情况
LL r = a[i] % k;
while(t --) r = r * 10 % k;
if(r == v) res --;
}
cout << res << endl;
return 0;
}