Эффективное нахождение всех делителей числа

поэтому я просто хочу найти все делители данного числа (за исключением самого числа). В настоящее время у меня есть это:

public static List<int> proper_divisors(int x)
{
    List<int> toreturn = new List<int>();
    toreturn.Add(1);
    int i = 0;
    int j=1;
    int z = 0;
    while (primes.ElementAt(i) < Math.Sqrt(x))
    {
        if (x % primes.ElementAt(i) == 0)
        {
            toreturn.Add(primes.ElementAt(i));
            toreturn.Add(x / primes.ElementAt(i));
            j = 2;
            z = (int)Math.Pow(primes.ElementAt(i), 2);
            while (z < x)
            {
                if (x % z == 0)
                {
                    toreturn.Add(z);
                    toreturn.Add(x / z);
                    j++;
                    z = (int)Math.Pow(primes.ElementAt(i), j);
                }
                else
                {
                    z = x;
                }
            }
        }
        i++;
    }
    toreturn = toreturn.Distinct().ToList<int>();
    return toreturn;
}

где простые числа-это список простых чисел (предположим, что он правильный и достаточно большой). Алгоритм работает в том смысле, что он находит все простые факторы, но не все факторы (т. е. данный 34534, он возвращает {1,2,17267,31,1114}, но пропускает {62, 557}, поскольку 62 является комбинацией, и, следовательно, пропускает 557.

Я также пробовал просто получая простые множители числа, но я не знаю как преобразовать это в список всех правильных комбинаций.

код для этого алгоритма выглядит следующим образом:

public static List<int> prime_factors(int x)
{
    List<int> toreturn = new List<int>();
    int i = 0;
    while (primes.ElementAt(i) <= x)
    {
        if (x % primes.ElementAt(i) == 0)
        {
            toreturn.Add(primes.ElementAt(i));
            x = x / primes.ElementAt(i);
        }
        else
        {
            i++;
        }
    }
    return toreturn;
}

любые идеи о том, как исправить первый, или как создать список комбинаций из второго (я бы предпочел, чтобы это было быстрее)?

9
c# factorization prime-factoring

2 ответов

так как у вас уже есть список простых факторов, то, что вы хотите сделать, это вычислить powerset этого списка.

теперь одна проблема заключается в том, что у вас могут быть дубликаты в списке (например, простые коэффициенты 20 = 2 * 2 * 5), но наборы не допускают дубликатов. Таким образом, мы можем сделать каждый элемент списка уникальным, проецируя его на структуру вида {x, y}, где x-простое число, а y-индекс простого числа в списке.

var all_primes = primes.Select((x, y) => new { x, y }).ToList();

теперь all_primes - это список формы {x, y}, где x-простое число, а y-индекс в списке.

затем мы вычисляем набор мощности (определение GetPowerSet ниже):

var power_set_primes = GetPowerSet(all_primes);

отсюда power_set_primes это IEnumerable<IEnumerable<T>> здесь T является анонимным типом {x, y} где x и y имеют тип int.

Далее мы вычисляем произведение каждого элемента питания set

foreach (var p in power_set_primes)
{
    var factor = p.Select(x => x.x).Aggregate(1, (x, y) => x * y);
    factors.Add(factor);
}

собираем все вместе:

var all_primes = primes.Select((x, y) => new { x, y }).ToList(); //assuming that primes contains duplicates.
var power_set_primes = GetPowerSet(all_primes);
var factors = new HashSet<int>();

foreach (var p in power_set_primes)
{
    var factor = p.Select(x => x.x).Aggregate(1, (x, y) => x * y);
    factors.Add(factor);
}

от http://rosettacode.org/wiki/Power_Set для реализаций powerset.

public IEnumerable<IEnumerable<T>> GetPowerSet<T>(List<T> list)
{
    return from m in Enumerable.Range(0, 1 << list.Count)
           select
               from i in Enumerable.Range(0, list.Count)
               where (m & (1 << i)) != 0
               select list[i];
}
7

был аналогичный вопрос до, который имеет интересное решение с использованием IEnumerable. Если вам нужны все делители, а не факторы, и если вы используете хотя бы C# 3.0, вы можете использовать что-то вроде этого:

static IEnumerable<int> GetDivisors(int n)
{
    return from a in Enumerable.Range(2, n / 2)
           where n % a == 0
           select a;                      
}

а затем используйте его так:

foreach(var divisor in GetDivisors(10))
    Console.WriteLine(divisor);

или, если вам нужен список, просто:

List<int> divisors = GetDivisors(10).ToList();
5