题目描述
给定一个长度为 nnn 的数组 A1,A2,⋅⋅⋅,AnA_1,A_2,⋅⋅⋅,A_nA1,A2,⋅⋅⋅,An。
你可以从中选出两个数 AiA_iAi 和 AjA_jAj(iii 不等于 jjj),然后将 AiA_iAi 和 AjA_jAj 一前一后拼成一个新的整数。
例如 12 和 345 可以拼成 12345 或 34512。
注意交换 AiA_iAi 和 AjA_jAj 的顺序总是被视为 2 种拼法,即便是 Ai=AjA_i=A_jAi=Aj 时。
请你计算有多少种拼法满足拼出的整数是 KKK 的倍数。
输入格式
第一行包含 2 个整数 nnn 和 KKK。
第二行包含 n 个整数 A1,A2,⋅⋅⋅,AnA_1,A_2,⋅⋅⋅,A_nA1,A2,⋅⋅⋅,An。
输出格式
一个整数代表答案。
数据范围
1≤n≤1051≤n≤10^51≤n≤105,
1≤K≤1051≤K≤10^51≤K≤105,
1≤Ai≤1091≤A_i≤10^91≤Ai≤109
输入样例
4 2
1 2 3 4
输出样例
6
算法思想(哈希)
根据题目描述,求AiA_iAi 和 AjA_jAj组成的整数AiAj‾\overline {A_iA_j}AiAj能够整除KKK的方案总数。进一步分析AiAj‾=Ai×10t+Aj\overline {A_iA_j}=A_i\times10^{t}+A_jAiAj=Ai×10t+Aj,其中ttt表示AjA_jAj的位数。那么对于每个AjA_jAj来说,要求Ai×10t+Aj≡0(modK)A_i\times10^{t}+A_j\equiv0\pmod{K}Ai×10t+Aj≡0(modK)的方案数,即求Ai×10t≡−Aj(modK)A_i\times10^{t}\equiv -A_j\pmod{K}Ai×10t≡−Aj(modK)的方案数。
通过上述分析,可以使用哈希表存储所有Ai×10t%KA_i\times10^{t} \% KAi×10t%K的方案数。不妨设v=Ai×10t%Kv = A_i\times10^{t} \% Kv=Ai×10t%K,那么h[t][v]h[t][v]h[t][v]表示Ai×10t%KA_i\times10^{t} \% KAi×10t%K值为vvv的方案数。由于1≤Ai≤1091≤A_i≤10^91≤Ai≤109,所以0≤t≤100\le t \le 100≤t≤10;而1≤K≤1051≤K≤10^51≤K≤105,所以0≤v≤100\le v \le 100≤v≤10。
那么对于每个AiA_iAi来说,只需要累加哈希表中h[t][v]h[t][v]h[t][v]的方案数即可,其中t=Ait = A_it=Ai的位数,v=(K−Ai%K)%Kv = (K-A_i \% K)\%Kv=(K−Ai%K)%K,即−A[i](modK)-A[i] \pmod{K}−A[i](modK)的值。
注意,如果Ai×10t%KA_i\times10^{t} \% KAi×10t%K的值也等于vvv,为了保证选出的两个数 AiA_iAi 和 AjA_jAj(iii 不等于 jjj),所以方案总数要减去1。
时间复杂度
需要遍历每个AiA_iAi,查询哈希表中h[t][v]h[t][v]h[t][v]的值,所以时间复杂度为O(n×t)=106O(n \times t) = 10^6O(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;
}