查找总和除以给定数字的数组的子数组的数量

Find numbers of subarray of an array whose sum is divided by given number

本文关键字:数组 数字 查找      更新时间:2023-10-16

我陷入了一个算法问题。请为我建议一些针对以下问题的有效算法。

问题是

查找其总和可被给定数字整除的子数组的数量。

我的工作

我做了一个算法,其复杂度是O(N^2),这里,N =数组的大小。

我的代码

#include<stdio.h>
using namespace std;
 main() {
    int N;
    int P;
    int T;
    int val;
    long long int count = 0;
    long long int answer = 0;
    scanf("%d", &T);
    //T = 20;
    for(int k = 1; k <= T; k++) {
        scanf("%d", &N);
        scanf("%d", &P);
        count = 0;
        answer = 0;
        for(int i = 0; i < N; i++) {
            scanf("%d", &val);
            count += val;
            workingArray[i] = count;
        }
        for(int length = 1; length <= N; length++) {
            for(int start = 0; start <= (N-length); start++) {
                if( start == 0 ) {
                    if(workingArray[start+length-1]%P == 0) answer++;
                }
                else if( (workingArray[start+length-1] - workingArray[start-1])%P == 0) answer++;
            }
        }
        printf("Case #%dn%lldn", k, answer);
    }
    return 0;
 }

对于给定的数量X...

基本思想:(带有非正式的正确性证明)

如果范围[a, b]中的数字之和能被 X 整除,则:

(∑i=1 to a-1input[i]) % X = (∑i=1 to binput[i]) % X

用不太数学的术语来说:

the sum from the first element to b = the sum from the first element to a
                                    + the sum of the elements between the two

所以:

the sum of the elements between the two = the sum from the first element to b
                                        - the sum from the first element to a

然后,如果右边的和在除以X时都有相同的余数,余数将抵消,两者之间的元素之和将被X整除。阐述:

C = the sum of the elements between the two
B = the sum from the first element to b
A = the sum from the first element to a

现在我们可以将B转换为形式PX + QA转换为形式RX + S,对于一些整数PQRS,带有0 <= Q, S < X。在这里,根据定义,QS将是BA各自的余数除以X

然后我们有:

C = (PX + Q) - (RX + S)
C = PX + Q - RX - S
C = PX - RX + Q - S
C = (P-R)X + Q - S

显然,(P-R)X可以被X整除(结果很简单(P-R))。现在我们只需要Q - S能被X整除,但是,既然0 <= Q, S < X,它们就需要相等。

例:

B = 13A = 7X = 3

这里B % X = 1A % X = 1.

我们可以将B重写为4*3 + 1A重写为2*3 + 1.

然后C = 4*3 + 1 - 2*3 - 1 = 2*3,它能被3整除。

高级方法:

构造一个 key -> value 的哈希映射,其中每个值表示您可以从数组的开头开始到在某个给定位置结束的几种方式,加起来为 sum mod X = key(请参阅"Mod 3"行和下面的示例中的映射值)。

现在,基于上面的逻辑,我们知道如果两个子数组分别从位置开始和结束于位置 ab,两者都具有相同的sum mod X,子数组[a, b]将被 X 整除。

因此,哈希映射中的每个值都表示一组可能的起点和终点的大小,这将为我们提供一个可被X整除的子数组(任何点可以是起点或终点)。

选择这些起点和终点的可能方法的数量很简单
value choose 2 = value!/(2*(value-2)!)(如果值为 1,则为 0)。

因此,我们计算哈希映射中的每个值,并将它们全部相加以得到可被X整除的子数组的数量。

算法:

构造一个哈希映射,它将存储到目前为止所有数字的累积总和,mod X映射到该余数值出现的频率(以预期O(n)构造)。

0 的值增加 1 - 这对应于数组的开头。

将计数初始化为 0。

对于哈希映射中的每个值,将value!/(2*(value-2)!)添加到计数中。

计数是所需的值。

运行时间:

预期O(n) .

例:

Input:    0  5  3  8  2  1
X = 3
Sum:   0  0  5  8 16 18 19
Mod 3: 0  0  2  2  1  0  1
Map:
  0 -> 3
  1 -> 2
  2 -> 2
Count = 3! / 2*(3-2)! = 3  +
        2! / 2*(2-2)! = 1  +
        2! / 2*(2-2)! = 1
      = 5

子数组将是:

0  5  3  8  2  1
-                     0                 =  0 % 3 = 0
-------------         0 + 5 + 3 + 8 + 2 = 18 % 3 = 0
   ----------         5 + 3 + 8 + 2     = 18 % 3 = 0
      -               3                 =  3 % 3 = 0
            ----      2 + 1             =  3 % 3 = 0

我可能有一个更简单的解决方案。 在 O(n) 时间和 O(n+k) 空间中。 其中n是数组的大小,k是我们正在检查可除性的数字。

将数组视为 A[n],数字为 K

  1. 创建另一个数组 SUM_TILL_NOW[n]。
  2. 对于每个 A[i] 填充SUM_TILL_NOW [i]= SUM_TILL_NOW[i-1]+A[i] %K;(SUM_TILL_NOW[0]= A[0])
  3. 在这个新数组中找到两个相等的数字。

为此,创建一个大小为 K 的新数组 CHECK[]。

遍历SUM_TILL_NOW数组并检查是否设置了 CHECK[SUM_TILL_NOW[i]]。

如果未设置为 i。

否则 CHECK[SUM_TILL_NOW[i]],i 是总和可被 K 整除的子集。

下面是一个相同的 c++ 函数。

#include <iostream>
#include <string.h>
using namespace std;
void printrange(int* A, int N, int K)
{
    int STN[N], C[K];
    memset(C, -1, K*sizeof(int));
    int i;
    int sum=A[0];
    STN[0]= (A[0]%K);
    for (i= 1; i< N; i++)
    {
        sum+= A[i];
        STN[i]= sum%K;
    }
    for(i=0; i< N; i++)
    {
        if(C[STN[i]] == -1)
            C[STN[i]] =i;
        else
        {
            cout<< C[STN[i]]+1 <<" "<< i;
            break;
        }
    }
}
int main()
{
    int A[]= {6, 9, 2, 1, 8, 6, 2, 5};
    printrange(A, sizeof(A)/sizeof(A[0]), 7);
    return 0;
}