重数计算的最佳有效解决方案?

Optimal and efficient solution for the heavy number calculation?

我需要找出两个整数 AB 之间的重整数的数量,其中总是 A <= B

An integer is considered heavy whenever the average of it's digit is larger than 7.

For example: 9878 is considered heavy, because (9 + 8 + 7 + 8)/4 = 8 , while 1111 is not, since (1 + 1 + 1 + 1)/4 = 1.

我有下面的解决方案,但它绝对糟糕,当 运行 有大量输入时它会超时。我该怎么做才能提高效率?

int countHeavy(int A, int B) {
    int countHeavy = 0;

    while(A <= B){
        if(averageOfDigits(A) > 7){
            countHeavy++;
        }
        A++;
    }

    return countHeavy;
}
 
float averageOfDigits(int a) {
    float result = 0;
    int count = 0;

    while (a > 0) {
        result += (a % 10);
        count++;
        a = a / 10;
    }

    return result / count;
}

我对问题的看法与你不同。我的看法是问题是基于数字的 10 进制表示,因此您应该做的第一件事是将数字放入 10 进制表示中。可能有更好的方法,但是 Java 字符串表示以 10 为基数的整数,所以我使用了它们。将单个字符转换为整数实际上非常快,所以这并不会花费太多时间。

最重要的是,你在这件事上的计算永远不需要使用除法或浮点数。问题的核心是仅与整数有关。数字(整数)中的所有数字(整数)加起来是否大于或等于七(整数)乘以数字(整数)的数目?

警告 - 我并不是说这是最快的方法,但这可能比您原来的方法更快。

这是我的代码:

package heavyNum;

public class HeavyNum
{
    public static void main(String[] args)
    {
        HeavyNum hn = new HeavyNum();
        long startTime = System.currentTimeMillis();
        hn.countHeavy(100000000, 1);
        long endTime = System.currentTimeMillis();
        System.out.println("Time elapsed: "+(endTime- startTime));
    }

    private void countHeavy(int A, int B)
    {
        int heavyFound = 0;
        for(int i = B+1; i < A; i++)
        {
            if(isHeavy(i))
                heavyFound++;
        }
        System.out.println("Found "+heavyFound+" heavy numbers");
    }

    private boolean isHeavy(int i)
    {
        String asString = Integer.valueOf(i).toString();
        int length = asString.length();
        int dividingLine = length * 7, currTotal = 0, counter = 0;
        while(counter < length)
        {
            currTotal += Character.getNumericValue(asString.charAt(counter++));
        }
        return currTotal > dividingLine;
    }
}

感谢 this SO Question for how I get the number of digits in an integer and this SO Question 如何在 java

中快速将字符转换为整数

运行 在一台功能强大的计算机上,没有针对 1 到 100,000,000 之间的数字的调试器导致此输出:

Found 569484 heavy numbers

Time elapsed: 6985

编辑:我最初是在寻找位数大于或等于位数的 7 倍的数字。我之前在 7025 毫秒内得到了 843,453 个数字的结果。

您确实可以使用字符串来获取位数,然后将各个数字的值相加以查看它们是否 sum > 7 * length,正如 Jeutnarg 似乎所做的那样。我拿了他的代码并添加了我自己的代码,简单 isHeavyRV(int):

private boolean isHeavyRV(int i)
{
    int sum = 0, count = 0;
    while (i > 0)
    {
        sum += i % 10;
        count++;
        i = i / 10;
    }
    return sum >= count * 7;
}

现在,而不是

        if(isHeavy(i))

我试过了

        if(isHeavyRV(i))

我实际上首先使用字符串测试了他对 isHeavy() 的实现,并且 运行 在我的机器(一台旧 iMac)上用了 12388 毫秒,它找到了 843453 个重数。

使用我的实现,我发现了完全相同数量的重数,但只用了 5416 毫秒。

字符串可能很快,但它们无法击败一个简单的循环,基本上可以完成 Integer.toString(i, 10) 所做的事情,但没有绕行字符串。

当您将一个数字加 1 时,您将增加一位数字,并将所有较小的数字更改为零。如果递增从重数变为非重数,那是因为太多的低位数字被归零了。在这种情况下,无需检查中间的所有数字就可以很容易地找到下一个重数:

public class CountHeavy
{
    public static void main(String[] args)
    {
        long startTime = System.currentTimeMillis();
        int numHeavy = countHeavy(1, 100000000);
        long endTime = System.currentTimeMillis();
        System.out.printf("Found %d heavy numbers between 1 and 100000000\n", numHeavy);
        System.out.println("Time elapsed: "+(endTime- startTime)+" ms");
    }

    static int countHeavy(int from, int to)
    {
        int numdigits=1;
        int maxatdigits=9;
        int numFound = 0;
        if (from<1)
        {
            from=1;
        }
        for(int i = from; i < to;)
        {
            //keep track of number of digits in i
            while (i > maxatdigits)
            {
                long newmax = 10L*maxatdigits+9;
                maxatdigits = (int)Math.min(Integer.MAX_VALUE, newmax);
                ++numdigits;
            }
            //get sum of digits
            int digitsum=0;
            for(int digits=i;digits>0;digits/=10)
            {
                digitsum+=(digits%10);
            }

            //calculate a step size that increments the first non-zero digit
            int step=1;
            int stepzeros=0;
            while(step <= (Integer.MAX_VALUE/10) && to-i >= step*10 && i%(step*10) == 0)
            {
                step*=10;
                stepzeros+=1;
            }
            //step is a 1 followed stepzeros zeros

            //how much is our sum too small by?
            int need = numdigits*7+1 - digitsum;
            if (need <= 0)
            {
                //already have enough.  All the numbers between i and i+step are heavy
                numFound+=step;
            }
            else if (need <= stepzeros*9)
            {
                //increment to the smallest possible heavy number. This puts all the
                //needed sum in the lowest-order digits
                step = need%9;
                for(;need >= 9;need-=9)
                {
                    step = step*10+9;
                }
            }
            //else there are no heavy numbers between i and i+step
            i+=step;
        }
        return numFound;
    }
}

Found 569484 heavy numbers between 1 and 100000000

Time elapsed: 31 ms

请注意,答案与@JeutNarg 的答案不同,因为您要求平均值 > 7,而不是平均值 >= 7。

通过查找计算数字 table

您可以生成一个 table 来存储有多少个具有 d 位的整数,其数字之和大于数字 x 。然后,您可以快速查找在 10、100、1000 ... 整数的任意范围内有多少重数。这些table只包含9×d个值,所以占用的space很少,可以快速生成。

然后,要检查范围 A-B,其中 B 有 d 个数字,您为 1 到 d[=66= 构建 tables ]-1 位,然后将范围 A-B 分成 10、100、1000 ... 的块,并查找 table 中的值,例如对于范围 A = 782,B = 4321:

   RANGE      DIGITS  TARGET     LOOKUP      VALUE

 782 -  789     78x    > 6    table[1][ 6]     3    <- incomplete range: 2-9
 790 -  799     79x    > 5    table[1][ 5]     4
 800 -  899     8xx    >13    table[2][13]    15
 900 -  999     9xx    >12    table[2][12]    21
1000 - 1999    1xxx    >27    table[3][27]     0
2000 - 2999    2xxx    >26    table[3][26]     1
3000 - 3999    3xxx    >25    table[3][25]     4
4000 - 4099    40xx    >24    impossible       0
4100 - 4199    41xx    >23    impossible       0
4200 - 4299    42xx    >22    impossible       0
4300 - 4309    430x    >21    impossible       0
4310 - 4319    431x    >20    impossible       0
4320 - 4321    432x    >19    impossible       0    <- incomplete range: 0-1
                                              --
                                              48

如果第一个和最后一个范围不完整(不是*0 - *9),请对照目标检查起始值或结束值。 (在示例中,2 不大于 6,因此所有 3 个重数都包含在范围内。)

生成查询 table

对于一位十进制整数,大于值x的整数n的个数是:

x:  0  1  2  3  4  5  6  7  8  9
n:  9  8  7  6  5  4  3  2  1  0

如您所见,这很容易通过取 n = 9-x 来计算。

对于2位十进制整数,整数n的位数之和大于值x的个数是:

x:   0  1  2  3  4  5  6  7  8  9 10 11 12 13 14 15 16 17 18
n:  99 97 94 90 85 79 72 64 55 45 36 28 21 15 10  6  3  1  0

对于三位十进制整数,整数n的位数之和大于值x的个数是:

x:   0   1   2   3   4   5   6   7   8   9  10  11  12  13  14  15  16  17  18  19  20  21  22  23  24  25  26  27
n: 999 996 990 980 965 944 916 880 835 780 717 648 575 500 425 352 283 220 165 120  84  56  35  20  10   4   1   0

这些序列中的每一个都可以从前一个序列生成:从值 10d 开始,然后从该值中减去前一个序列(跳过第一个零)。例如。要从 2 位数字的序列生成 3 位数字的序列,从 103 = 1000 开始,然后:

 0. 1000 -   1      = 999
 1.  999 -   3      = 996
 2.  996 -   6      = 990
 3.  990 -  10      = 980
 4.  980 -  15      = 965
 5.  965 -  21      = 944
 6.  944 -  28      = 916
 7.  916 -  36      = 880
 8.  880 -  45      = 835
 9.  835 -  55      = 780
10.  780 -  64 +  1 = 717  <- after 10 steps, start adding the previous sequence again
11.  717 -  72 +  3 = 648
12.  648 -  79 +  6 = 575
13.  575 -  85 + 10 = 500
14.  500 -  90 + 15 = 425
15.  425 -  94 + 21 = 352
16.  352 -  97 + 28 = 283
17.  283 -  99 + 36 = 220
18.  220 - 100 + 45 = 165  <- at the end of the sequence, keep subtracting 10^(d-1)
19.  165 - 100 + 55 = 120
20.  120 - 100 + 64 =  84
21.   84 - 100 + 72 =  56
22.   56 - 100 + 79 =  35
23.   35 - 100 + 85 =  20
24.   20 - 100 + 90 =  10
25.   10 - 100 + 94 =   4
26.    4 - 100 + 97 =   1
27.    1 - 100 + 99 =   0

顺便说一下,如果 "heavy" 数字定义为 7 以外的值,您可以使用相同的 tables。


代码示例

下面是演示该方法的 Java 脚本代码片段(我不会说话 Java)。它非常未经优化,但它在不到 0.07 毫秒的时间内完成了 0→100,000,000 示例。它也适用于 7 以外的权重。转换为 Java,它应该很容易击败任何实际运行数字并检查其权重的算法。

function countHeavy(A, B, weight) {
    var a = decimalDigits(A), b = decimalDigits(B);        // create arrays
    while (a.length < b.length) a.push(0);                 // add leading zeros
    var digits = b.length, table = weightTable();          // create table
    var count = 0, diff = B - A + 1, d = 0;                // calculate range
    for (var i = digits - 1; i >= 0; i--) if (a[i]) d = i; // lowest non-0 digit
    while (diff) {                                         // increment a until a=b
        while (a[d] == 10) {                               // move to higher digit
            a[d++] = 0;
            ++a[d];                                        // carry 1
        }
        var step = Math.pow(10, d);                        // value of digit d
        if (step <= diff) {
            diff -= step;
            count += increment(d);                         // increment digit d
        }
        else --d;                                          // move to lower digit
    }
    return count;

    function weightTable() {                               // see above for details
        var t = [[],[9,8,7,6,5,4,3,2,1,0]];
        for (var i = 2; i < digits; i++) {
            var total = Math.pow(10, i), final = total / 10;
            t[i] = [];
            for (var j = 9 * i; total > 0; --j) {
                if (j > 9) total -= t[i - 1][j - 10]; else total -= final;
                if (j < 9 * (i - 1)) total += t[i - 1][j];
                t[i].push(total);
            }
        }
        return t;
    }
    function increment(d) {
        var sum = 0, size = digits;
        for (var i = digits - 1; i >= d; i--) {
            if (a[i] == 0 && i == size - 1) size = i;      // count used digits
            sum += a[i];                                   // sum of digits
        }
        ++a[d];
        var target = weight * size - sum;
        if (d == 0) return (target < 0) ? 1 : 0;           // if d is lowest digit
        if (target < 0) return table[d][0] + 1;            // whole range is heavy
        return (target > 9 * d) ? 0 : table[d][target];    // use look-up table
    }
    function decimalDigits(n) {
        var array = [];
        do {array.push(n % 10);
            n = Math.floor(n / 10);
        } while (n);
        return array;
    }
}
document.write("0 &rarr; 100,000,000 = " + countHeavy(0, 100000000, 7) + "<br>");
document.write("782 &rarr; 4321 = " + countHeavy(782, 4321, 7) + "<br>");
document.write("782 &rarr; 4321 = " + countHeavy(782, 4321, 5) + " (weight: 5)");

我真的很喜欢@m69 的 post 所以我写了受它启发的实现。 table 的创建不是那么优雅,但是有效。对于 n+1 位长整数,我从 n 位长整数中(最多)求和 10 个值,每个数字 0-9 一个。

我使用这种简化来避免任意范围计算:

countHeavy(A, B) = countHeavy(0, B) - countHeavy(0, A-1)

结果分两次循环计算。一个用于比给定数字短的数字,一个用于其余数字。我无法轻松合并它们。 getResult只是通过范围检查查找 table,其余代码应该很明显。

public class HeavyNumbers {
    private static int maxDigits = String.valueOf(Long.MAX_VALUE).length();
    private int[][] table = null;

    public HeavyNumbers(){
        table = new int[maxDigits + 1][];
        table[0] = new int[]{1};

        for (int s = 1; s < maxDigits + 1; ++s) {
            table[s] = new int[s * 9 + 1];
            for (int k = 0; k < table[s].length; ++k) {
                for (int d = 0; d < 10; ++d) {
                    if (table[s - 1].length > k - d) {
                        table[s][k] += table[s - 1][Math.max(0, k - d)];
                    }
                }
            }
        }
    }

    private int[] getNumberAsArray(long number) {
        int[] tmp = new int[maxDigits];
        int cnt = 0;

        while (number != 0) {
            int remainder = (int) (number % 10);
            tmp[cnt++] = remainder;
            number = number / 10;
        }

        int[] ret = new int[cnt];
        for (int i = 0; i < cnt; ++i) {
            ret[i] = tmp[i];
        }
        return ret;
    }

    private int getResult(int[] sum, int digits, int fixDigitSum, int heavyThreshold) {
        int target = heavyThreshold * digits - fixDigitSum + 1;
        if (target < sum.length) {
            return sum[Math.max(0, target)];
        }
        return 0;
    }

    public int getHeavyNumbersCount(long toNumberIncl, int heavyThreshold) {
        if (toNumberIncl <= 0) return 0;

        int[] numberAsArray = getNumberAsArray(toNumberIncl);

        int res = 0;

        for (int i = 0; i < numberAsArray.length - 1; ++i) {
            for (int d = 1; d < 10; ++d) {
                res += getResult(table[i], i + 1, d, heavyThreshold);
            }
        }

        int fixDigitSum = 0;
        int fromDigit = 1;
        for (int i = numberAsArray.length - 1; i >= 0; --i) {
            int toDigit = numberAsArray[i];
            if (i == 0) {
                toDigit++;
            }
            for (int d = fromDigit; d < toDigit; ++d) {
                res += getResult(table[i], numberAsArray.length, fixDigitSum + d, heavyThreshold);
            }

            fixDigitSum += numberAsArray[i];
            fromDigit = 0;
        }

        return res;
    }

    public int getHeavyNumbersCount(long fromIncl, long toIncl, int heavyThreshold) {
        return getHeavyNumbersCount(toIncl, heavyThreshold) -
                getHeavyNumbersCount(fromIncl - 1, heavyThreshold);
    }
}

是这样使用的:

HeavyNumbers h = new HeavyNumbers();
System.out.println( h.getHeavyNumbersCount(100000000,7));

打印出569484,table不初始化的重复计算时间小于1us

这是一个带有记忆的漂亮准系统递归,它逐一枚举固定数字的数字可能性。在计算相应的位数时,可以通过控制i的范围来设置AB

看起来很快(查看 20 位数字的结果)。

JavaScript代码:

var hash = {}
 
function f(k,soFar,count){
  if (k == 0){
    return 1;
  }

  var key = [k,soFar].join(",");
  
  if (hash[key]){
    return hash[key];
  }

  var res = 0;

  for (var i=Math.max(count==0?1:0,7*(k+count)+1-soFar-9*(k-1)); i<=9; i++){
    res += f(k-1,soFar+i,count+1);
  }

  return hash[key] = res;
}

// Output:

console.log(f(3,0,0)); // 56
hash = {};

console.log(f(6,0,0)); // 12313
hash = {};

console.log(f(20,0,0)); // 2224550892070475