只考加法的面试题

我们知道:

1+2 = 3；

4+5 = 9；

2+3+4 = 9。

等式的左边都是两个或两个以上连续的自然数相加，是不是所有的整数都可以写成这样的形式呢？

问题1：  对于一个64位正整数，输出它所有可能的连续自然数（两个以上）之和的算式。

问题2：  大家在测试上面程序的过程中，肯定会注意到有一些数字不能表达为一系列连续的自然数

之和，例如32好像就找不到。那么，这样的数字有什么规律呢？能否证明你的结论？

问题3： 在64位正整数范围内，子序列数目最多的数是哪一个？

解答1. 首先分析，

对于正整数N，如果表示成2个连续自然数相加，N = m + (m + 1) = 2m + 1，则N为奇数；

如果表示成3个连续的自然数相加，N = (m - 1) + m + (m + 1) = 3m，则N为3的倍数；

如果表示成4个连续的自然数相加，N = (m - 1) + m + (m + 1) + (m + 2) = 2 * (2m + 1)，则N为某奇数的2倍；

如果表示成5个连续的自然数相加，……，则N为5的倍数；

如果表示成5个连续的自然数相加，……，则N为某奇数的3倍；

……

 

已经找到规律了：对于N，

如果表示成偶数个(2 * k) 连续的自然数相加，      则N为k的倍数；

如果表示成奇数个(2 * k + 1)连续的自然数相加， 则N为(2 * k + 1)的倍数。

 

这样程序就好写啦：

#include <math.h>   
#include <stdio.h>   
#include <stdlib.h>   
  
bool AddSubN(__int64 N)  
{  
    if(N < 3)  
    {  
        printf("No sequences fit N./n");  
        return false;  
    }  
    bool bFind = false;  
    int num = 0;  
    printf("/n %I64d ./n", N);  
  
    __int64 maxLoop = sqrt(2 * N);  
    //当N = 1 + 2 + ... + m = m * (m + 1) / 2时的序列长度最大，   
    //所以maxLoop比这个小   
    //注意：这里可能得用大整数开方，我懒得写了   
  
    __int64 i, j, testN;  
            // 看N是否能被表达成 i 个连续自然数之和   
    for(i = 2; i <= maxLoop; i ++)  
    {  
        if(!(i & 0x1))  //如果 i 是偶数，则N为 i/2 的倍数，且 N /(i/2) 为奇数；   
        {  
            if((!(N % (i >> 1))) && ((N / (i >> 1))) & 1)  
            {  
                __int64 sub0 = (N / (i >> 1) - 1) / 2;    // i 个数中，中间偏左的一个   
                sub0 -= i / 2 - 1;                      // i 个数中最左边的一个   
                testN = 0;  
                for(j = 0; j < i; j ++)  
                    testN += sub0 + j;  
                    //打印出来，测试是否跟输入的N一样   
                printf("/n %I64d = ", N, testN);  
  
                    //打印连续自然数序列   
                for(j = 0; j < i - 1; j ++)  
                    printf("%I64d + ", sub0 + j);  
                printf("%I64d ./n", sub0 + j);  
  
                if(!bFind)  
                    bFind = true;  
                num ++;  
  
            }  
        }  
        else        // 如果 i 是奇数，则N为 i 的倍数   
        {  
            if(!(N % i))  
            {  
                __int64 sub0 = N / i;  
                sub0 -= i / 2;          //找到i个数中最左边的一个   
  
                testN = 0;  
                for(j = 0; j < i; j ++)  
                    testN += sub0 + j;  
                //打印出来，测试是否跟输入的N一样   
                printf("/n %I64d = ", N, testN);  
                  
                //打印连续自然数序列   
                for(j = 0; j < i - 1; j ++)  
                    printf("%I64d + ", sub0 + j);  
                printf("%I64d ./n", sub0 + j);  
                  
                if(!bFind)  
                    bFind = true;  
  
                num ++;  
            }  
        }  
    }  
    if(!bFind)  
    {  
        printf("No sequences fit N./n");  
        return false;  
    }  
  
    printf("/n---------------/n %d sequences fit N found ./n", num);  
  
    return true;  
}  
  
void main()  
{  
    __int64 N = 3;  
    while(N >= 1)  
    {  
        printf("Please input a int64 number. (input 0 or -1 to escape)/n");  
        scanf("%I64d", &N);  
        system("cls");  
        AddSubN(N);  
        printf("/n/n/n");  
        system("pause");  
        system("cls");  
    }  
}

 

解答2. 从分析中可以看出，N的因式分解必然要有奇数。证明如下：

 a. 首先证明，只要N的因式分解中有奇数，N就能表示为自然数连和。

     如果N的因式分解中有奇数，假设为s，且N= k * s，

     如果k > s/2，

     则 N可以表示为这个序列的和：k - (s / 2), k - (s / 2) + 1 ... k + (s / 2);

     例如 54 = 6 x 9, 可表示为 2 + 3 + 4 + 5 + 6 + 7 + 8 + 9 + 10；（6 x 9，偶数为6，则6在这9个自然数的序列的中间）

     如果 k < s/2，

     则N可以表示为这个序列的和：(s + 1)/2 - k, (s+1)/2 - k + 1,..., (s + 1)/2 + k - 1;

     如54 = 2 x 27，则可表示为 12 + 13 + 14 + 15；（中间两个数 13 + 14 = 27）

 b. 再证明：如果N的因式分解中没有奇数，则N不能表示成连和的形式。 

     反证：如果能，假设N能表示成k个自然数的连和。

     如果k为偶数，设这k个数的中数（中间偏左的一个）为m，将这k个数收尾相加得到k/2个自然数的序列，

     易证这k/2个自然数都等于2m+1，

     那么N = (2m + 1) * k / 2，与N的因式分解中没有奇数矛盾；

 

     如果k为奇数，设中数为n，那么N = k * n，与N的因式分解中没有奇数矛盾。

     综上，得证。 

 ------------------------------------------------------------------------------------------------------------

假设N的质因子分解中有k种奇数（不是k个）记为odd[1,..,k]，odd[i]有n[i]个，那么N能表示的序列数为：

    Ns = (n[0] + 1) * (n[2] + 1) * (n[3] + 1) * ... * (n[k] + 1) - 1

 

验证：3*3*3*3*5*5*7*7*7 = 694575，我用程序打印，输出了59个。

这里有4个3，2个5，3个7，(4+1)*(2+1)*(3+1) - 1 = 60 - 1 = 59，验证成功。而694575乘以2的任意幂次方，所得结果都是59.

显然，拆分个数，只与奇质因子的数目有关。

2 ^ 64 = 1.8e19

3 * 5 *7 *11 *13 *17 * 19 *23 *29 *31 * 37 *41 * 43 *47 *53 = 1.6e19

 

假设N是有最多因子个数的最小64位奇数，设 N = 3^a3 * 5^a5 * 7^a7 …

则一定有 a3 >= a5 >= a7 … 否则只要交换不满足条件的那两个数，得到相同因子个数但比N更小的数，这与假设矛盾。

  S = 2 ^ 64 = 1.8e19

M=3*5*7*11*13*17*19*23*29*31*37*41*43*47*53=1.6e19（因子个数2^15）

因而，N的最大质因子一定小等于53

 

由S / (M / 53) = 60  可将60拆分成3^3（因子数5*2^13）  3^2 * 5（因子数3*2^14）

可得局部最优解：R1 = 3^3 *5^2 *7*11*13*17*19*23*29*31*37*41*43*47

如果N不等于R1，则a47 = 0（要将S / (M / 53/47)) = 2820 拿出来拆分）

  若N包含k个质因子a， t为满足a^t > 47（显然t >= 2）的最小整数，则 k < 2*t-1

（否则若将t个a拆分成47，由 (k+1)*1 – (k-t+1) * 2 = 2*t-k-1 <=0，

可知拆分后得到的数更优，与N最优矛盾）。

因此a3 <=2*4-2=6，

a5 <= 2*3 – 2 = 4, 

a7 <= 2*2-2 = 2

a11 <= 2*2-2 = 2

 

若a7 <=1， 则a3<=4，否则可以将2个3拆成1个7，得到更优解。由

S/(3^4*5^4)/ (7*11*13*17*19*23*29*31*37*41)  = 35

（能得到的最多因子个数为25*2^10 < 3*2^14不是最优解）

因而 a7 = 2

 

（ 若az = 2, ax = a, ay =b 且 z > x * y，若不能将 z拆分成 x * y，则有

   (a+1)*(b+1)*3 > (a+2)*(b+2)*2，即 (a-1)*(b-1) >= 7 ）

 

若a23=2则可将1个23拆成3和7，由 (1+a3)*3*3 – (1+a3+1)*4*2 = a3-7<0

可知得到的数更优，与假设矛盾，因而 a23<=1，

由于 S/(3^6*5^4)/(7*11*13*17*19)^2 = 387 > 23因而 一定含有因子23，a23 = 1

 

若a31=0，则 a5 = 2（否则，5*7合并成31，得到更优解）

由 2^64 / (3^6*(5*7*11*13*17*19)^2 * 23 * 29) = 14

可知，该情况下得到的最大数不是最优， 因而 a31 = 1

 

（若a17 =2则 a3>=5, a5=3 或 a3>=4 a5=4，否则可以将17拆分成3*5）

 

利用前面的结论，

  a3 >= a5 >= a7 …

  a3 <= 6  a5 <= 4  a7 = 2  a23 = 1  a31 = 1  a47 = 0

可在较短时间内搜索出满足上述条件的因子个数最多的奇数，再与局部最优解R1进行比较，就可以确定最优解。