Produkty, które są równe sumie i odwrotnie

22

Zabawna para ekwiwalentów wynosi 1 + 5 = 2,3 i 1,5 = 2 + 3 . Istnieje wiele takich, inna to 1 + 1 + 8 = 1 · 2 · 5 i 1 · 1 · 8 = 1 + 2 + 5 . Zasadniczo iloczyn n dodatnich liczb całkowitych jest równy sumie n dodatnich liczb całkowitych i odwrotnie.

W tym wyzwaniu musisz wygenerować wszystkie takie kombinacje dodatnich liczb całkowitych dla wejścia n> 1 , z wyłączeniem permutacji. Możesz je wydrukować w dowolnym rozsądnym formacie. Na przykład wszystkie możliwe rozwiązania dla n = 3 to:

(2, 2, 2) (1, 1, 6)
(1, 2, 3) (1, 2, 3)
(1, 3, 3) (1, 1, 7)
(1, 2, 5) (1, 1, 8)

Program, który może wygenerować najwięcej kombinacji dla najwyższej n w ciągu jednej minuty na moim 2 GB pamięci RAM , wygrywa 64-bitowy laptop Intel Ubuntu. Jeśli twoja odpowiedź zużywa więcej niż 2 GB pamięci RAM lub jest napisana w języku, którego nie mogę przetestować za pomocą bezpłatnego oprogramowania, nie uzyskam odpowiedzi. Przetestuję odpowiedzi za dwa tygodnie i wybiorę zwycięzcę. Oczywiście później można opublikować niekonkurencyjne odpowiedzi.

Ponieważ nie wiadomo, jakie są pełne zestawy rozwiązań dla wszystkich n , możesz publikować odpowiedzi, które generują niekompletne rozwiązania. Jeśli jednak inna odpowiedź wygeneruje (bardziej) kompletne rozwiązanie, nawet jeśli ich maksymalna liczba n jest mniejsza , odpowiedź wygrywa.


Aby to wyjaśnić, oto proces punktacji w celu ustalenia zwycięzcy:

  1. Przetestuję twój program przy n = 2, n = 3 itd. Przechowuję wszystkie twoje wyjścia i zatrzymuję się, gdy twój program zajmie więcej niż minutę lub więcej niż 2 GB pamięci RAM. Za każdym razem, gdy program jest uruchamiany dla danego wejścia n, zostanie zakończony, jeśli zajmie to więcej niż 1 minutę.

  2. Patrzę na wszystkie wyniki dla wszystkich programów dla n = 2. Jeśli program wytworzył mniej poprawne rozwiązania niż inne, program ten jest eliminowany.

  3. Powtórz krok 2 dla n = 3, n = 4 itd. ... Ostatni stojący program wygrywa.

orlp
źródło
1
Więc nie ma odpowiedzi w językach wyłącznych dla systemu Windows?
Conor O'Brien
3
Osobiście nie lubię kryteriów punktacji. Nie możemy wiedzieć, czy nasze rozwiązania będą działać i gdzie ustawić progi, dopóki nie uzyskamy wyników testów na twoim komputerze. Wydaje mi się, że prosty golf w golfa byłby lepszym pytaniem.
musicman523
2
Zakładam, że kodowanie na stałe jest niedozwolone. Ale potem to ograniczenie jest prawie niemożliwe do zaobserwowania
Luis Mendo
1
@ user202729 Nie, muszę wypróbować każdy program dla każdego n, aby zobaczyć, który program generuje więcej rozwiązań.
orlp
2
„Za dwa tygodnie” są 3 dni temu.
GB

Odpowiedzi:

4

C (gcc) , n = 50000000 przy 6499 wynikach w 59 s

Aby uniknąć wytwarzania ponad terabajta wyjściowego składającego się prawie w całości z 1s, sekwencja (powiedzmy) 49999995 1s jest skracana jako 1x49999995.

#include <stdio.h>
#include <stdlib.h>

static int n, *a1, k1 = 0, *a2, k2 = 0, s1, p1, *factor;

static void out() {
  if (s1 == p1) {
    for (int i = 0; i < k1 && i < k2; i++) {
      if (a1[i] < a2[i])
        return;
      else if (a1[i] > a2[i])
        break;
    }
  }

  for (int i = 0; i < k1; i++)
    printf("%d ", a1[i]);
  printf("1x%d | ", n - k1);
  for (int i = 0; i < k2; i++)
    printf("%d ", a2[i]);
  printf("1x%d\n", n - k2);
}

static void gen2(int p, int s, int m);

static void gen3(int p, int s, int m, int x, int q) {
  int r = s - n + k2 + 2;
  int d = factor[q];
  do {
    if (x * d <= m)
      x *= d;
    q /= d;
  } while (q % d == 0);
  do {
    if (q == 1) {
      a2[k2++] = x;
      gen2(p / x, s - x, x);
      k2--;
    } else {
      gen3(p, s, m, x, q);
    }
    if (x % d != 0)
      break;
    x /= d;
  } while (p / (x * q) >= r - x * q);
}

static void gen2(int p, int s, int m) {
  int n2 = n - k2;
  if (p == 1) {
    if (s == n2)
      out();
  } else if (n2 >= 1 && m > 1) {
    int r = s - n2 + 1;
    if (r < 2 || p < r)
      return;
    if (m > r)
      m = r;
    if (factor[p] <= m)
      gen3(p, s, m, 1, p);
  }
}

static void gen1(int p, int s, int m) {
  int n1 = n - k1;
  p1 = p;
  s1 = s + n1;
  gen2(s1, p1, s + n1 + 1 - n);
  if (n1 != 0) {
    int *p1 = &a1[k1++];
    for (int x = 2; x <= m && p * x <= s + x + n1 - 1; x++) {
      *p1 = x;
      gen1(p * x, s + x, x);
    }
    k1--;
  }
}

int main(int argc, char **argv) {
  if (argc < 2)
    return 1;
  n = atoi(argv[1]);
  if (n < 2)
    return 1;
  a1 = malloc(n * sizeof(int));
  a2 = malloc(n * sizeof(int));
  factor = calloc(4 * n - 1, sizeof(int));
  for (int p = 2; p < 4 * n - 1; p++)
    if (factor[p] == 0) {
      factor[p] = p;
      for (int i = p; i <= (4 * n - 2) / p; i++)
        factor[p * i] = p;
    } else if (factor[p] < factor[p / factor[p]]) {
      factor[p] = factor[p / factor[p]];
    }
  gen1(1, 0, 3 * n - 1);
  return 0;
}

Wypróbuj online!

Anders Kaseorg
źródło
3

Mathematica, n = 293 z 12 rozwiązaniami

OP zmieniło wyzwanie i prosi o dane wejściowe
Oto nowy kod, który przyjmuje dowolne n jako dane wejściowe
Dla n = 293 otrzymujesz 12 rozwiązań

If[#<5,Union[Sort/@Select[Tuples[{1,2,3,4,5,6,7,8,9},{#}],Tr@#==Times@@#&]],For[a=1,a<3,a++,For[b=a,b<3,b++,For[c=b,c<5,c++,For[d=c,d<10,d++,For[e=d,e<300,e++,If[Tr[s=Join[Table[1,#-5],{a,b,c,d,e}]]==Times@@s,Print[s]]]]]]]]&


wkład

[n]

Możesz przetestować ten algorytm na Wolfram Sandbox, który jest darmowym oprogramowaniem online.
Wystarczy kliknąć link, wkleić kod (ctrl + v), wkleić dane wejściowe na końcu kodu i nacisnąć shift + enter, aby uruchomić.
Wszystkie moje rozwiązania otrzymasz w kilka sekund

Oto także Wypróbuj online!w C ++ (gcc)
(wielkie dzięki @ThePirateBay za wsparcie i tłumaczenie mojego kodu na wolny język)

ten program generuje tylko rozwiązania o postaci {a, b, c} {a, b, c}
co oznacza a + b + c = a * b * c

Obliczenie zajmuje 1 sekundę

dwanaście rozwiązań to:

{1,1 ..., 1,1,1,2,293} {1,1 ..., 1,1,1,2,293}
{1,1 ..., 1,1,1,3,147} {1 , 1 ..., 1,1,1,3,147}
{1,1 ..., 1,1,1,5,74} {1,1 ..., 1,1,1,5,74}
{1,1 ..., 1,1,2,2,98} {1,1 ..., 1,1,2,2,98}
{1,1 ..., 1,1,2, 3,59} {1,1 ..., 1,1,2,3,59}
{1,1 ..., 1,1,2,5,33} {1,1 ..., 1, 1,2,5,33}
{1,1 ..., 1,1,2,7,23} {1,1 ..., 1,1,2,7,23}
{1,1 .. ., 1,1,2,8,20} {1,1 ..., 1,1,2,8,20}
{1,1 ..., 1,1,3,3,37} {1 , 1 ..., 1,1,3,3,37}
{1,1 ..., 1,1,3,4,27} {1,1 ..., 1,1,3,4, 27}
{1,1 ..., 1,1,3,7,15} {1,1 ..., 1,1,3,7,15}
{1,1 ..., 1,2, 2,6,13} {1,1 ..., 1,2,2,6,13}

J42161217
źródło
1
„Jeśli twoja odpowiedź [...] jest napisana w języku, którego nie mogę przetestować za pomocą swobodnie dostępnego oprogramowania, nie uzyskam odpowiedzi.”
Leaky Nun
4
@GB „możesz publikować odpowiedzi, które generują niekompletne rozwiązania”
202729,
1
mój program „.. generuje najwięcej kombinacji dla najwyższego n w ciągu jednej minuty”. Nie jest zakodowany. Po prostu znajduje pierwsze 12 „najłatwiejszych” rozwiązań w mniej niż minutę
J42161217
1
Mogło być bardziej zrozumiałe, że n miało być wejściem. Wyjaśniłem to teraz. Nie wygląda na to, że twój program pobiera dane wejściowe n .
orlp
2
@orlp Naprawiono! Mój program przyjmuje dowolne n jako dane wejściowe. Dla n = 293 otrzymujesz 12 rozwiązań. anuluj głosowanie, jeśli wszystko działa!
J42161217,
2

Python 2 , n = 175, 28 wyników w 59s

Spowolniłem go, stosując współczynnik redukcji 2, ale dostaje więcej rozwiązań, zaczynając od n = 83

Otrzymuję wyniki dla n do 92 na TIO w jednym przebiegu.

def submats(n, r):
    if n == r:
        return [[]]
    elif r > 6:
        base = 1
    else:
        base = 2
    mx = max(base, int(n*2**(1-r)))

    mats = []
    subs = submats(n, r+1)
    for m in subs:
        if m:
            mn = m[-1]
        else:
            mn = 1
        for i in range(mn, mx + 1):
            if i * mn < 3*n:
                mats += [m + [i]]
    return mats

def mats(n):
    subs = []
    for sub in submats(n, 0):
        sum = 0
        prod = 1
        for m in sub:
            sum += m
            prod *= m
        if prod > n and prod < n*3:
            subs += [[sub, sum, prod]]
    return subs

def sols(n):
    mat = mats(n)
    sol = [
        [[1]*(n-1)+[3*n-1],[1]*(n-2)+[2,2*n-1]],
    ]
    if n > 2:
        sol += [[[1]*(n-1)+[2*n+1],[1]*(n-2)+[3,n]]]
    for first in mat:
        for second in mat:
            if first[2] == second[1] and first[1] == second[2] and [second[0], first[0]] not in sol:
                sol += [[first[0], second[0]]];
    return sol

Wypróbuj online!

GB
źródło
1
„zachowaj 5 elementów [1..2] i ogranicz 3n ...” Cieszę się, że polubiłeś mój algorytm ;-)
J42161217
Zrobiłem już coś podobnego w wersji Ruby, a teraz próbuję usunąć to ograniczenie.
GB
Dla danego n, ile rozwiązań jest zakodowanych na stałe w twoim algorytmie?
J42161217,
Niezbyt zakodowane: 2 standardowe rozwiązania można wygenerować za pomocą skrótu (z wyjątkiem n = 2, gdzie są one tą samą kombinacją), a pomijając je, mogę ograniczyć zakres do 2n zamiast 3n. Jeśli zostanie to uznane za zakodowane, zmienię to.
GB
1
Dla 61 mój wynik to 28, pamiętam, że to 27 ... Być może popełniłem błąd
RosLuP
1

Rubin , n = 12 otrzymuje 6 rozwiązań

Przynajmniej w TIO, zwykłe wyniki dla 1 do 11

->n{
  arr=[*1..n*3].product(*(0..n-2).map{|x|
    [*1..[n/3**x,2].max]|[1]
  }).select{|a|
    a.count(1) >= n-4
  }.map(&:sort).uniq
  arr.product(arr).map(&:sort).uniq.select{|r|
    r[0].reduce(&:+) == r[1].reduce(&:*) &&
    r[0].reduce(&:*) == r[1].reduce(&:+)
  }
}

Wypróbuj online!

Pobiera 10 wyników poniżej minuty dla n = 13 na moim laptopie.

GB
źródło
1

Mathematica, n = 19 z 11 rozwiązaniami

to moja nowa odpowiedź zgodnie z nowymi kryteriami PO

(SOL = {};
For[a = 1, a < 3, a++, 
For[b = a, b < 3, b++, 
For[c = b, c < 5, c++, 
 For[d = c, d < 6, d++, 
  For[e = d, e < 3#, e++, 
   For[k = 1, k < 3, k++, 
    For[l = k, l < 3, l++, 
     For[m = l, m < 5, m++, 
      For[n = m, n < 6, n++, For[o = n, o < 3#, o++,
        s = Join[Table[1, # - 5], {a, b, c, d, e}];
        t = Join[Table[1, # - 5], {k, l, m, n, o}];            
        If[Tr[s[[-# ;;]]] == Times @@ t[[-# ;;]] && 
          Tr[t[[-# ;;]]] == Times @@ s[[-# ;;]], 
         AppendTo[SOL,{s[[-#;;]],t[[-#;;]]}]]]]]]]]]]]];
Union[SortBy[#,Last]&/@SOL])&

jeśli wpiszesz [n] na końcu, program wyświetli rozwiązania

oto moje wyniki (na moim starym laptopie 64-bitowym 2,4 GHz)

n-> rozwiązania
2 -> 2
3 -> 4
4 -> 3
5 -> 5
6 -> 4
7 -> 6
8 -> 5
9 -> 7
10 -> 7
11 -> 8
12 -> 6 (w 17 sekund)
13 -> 10 (w 20 sekund)
14 -> 7 (w 25 sekund)
15 -> 7 (w 29 sekund)
16 -> 9 (w 34 sekund)
17 -> 10 (w 39 sekund)
18 - > 9 (w 45 sekund)
19 -> 11 (w 51 sekund)
20 -> 7 (w 58 sekund)

J42161217
źródło
1

Haskell, wiele rozwiązań szybko

import System.Environment

pr n v = prh n v v

prh 1 v l = [ [v] | v<=l ]
prh n 1 _ = [ take n $ repeat 1 ]
prh _ _ 1 = []
prh n v l = [ d:r | d <-[2..l], v `mod` d == 0, r <- prh (n-1) (v`div`d) d ]

wo n v = [ (c,k) | c <- pr n v, let s = sum c, s>=v,
                   k <- pr n s, sum k == v, s>v || c>=k ]

f n = concatMap (wo n) [n+1..3*n]

main = do [ inp ] <- getArgs
          let results = zip [1..] $ f (read inp)
          mapM_ (\(n,s) -> putStrLn $ (show n) ++ ": " ++ (show s)) results

foblicza rozwiązania, mainfunkcja dodaje pobieranie danych z wiersza poleceń oraz pewne formatowanie i liczenie.

Christian Sievers
źródło
Skompiluj w ten sposób: ghc -O3 -o prodsum prodsum.hsi uruchom z argumentem wiersza poleceń:./prodsum 6
Christian Sievers
0

Haskell , n = 10 z 2 roztworami


import           Data.List

removeDups = foldl' (\seen x -> if x `elem` seen then seen else x : seen) []
removeDups' = foldl' (\seen x -> if x `elem` seen then seen else x : seen) []

f n= removeDups $ map sort filterSums
  where maxNumber = 4
        func x y = if (((fst x) == (fst.snd$y)) && ((fst y) == (fst.snd$x)))
                     then [(snd.snd$x),(snd.snd$y)]
                     else [[],[]]
        pOf = removeDups' $ (map sort (mapM (const [1..maxNumber]) [1..n]))
        sumOf = map (\x->((sum x),((product x), x))) pOf
        filterSums = filter (\x-> not$(x == [[],[]])) (funcsumOfsumOf)

Działa to jak bzdury, ale przynajmniej to naprawiłem, więc teraz podejmuję wyzwanie!

Wypróbuj online!

wałek klonowy
źródło
dla n = 2 otrzymujesz [„[3,3] [2,3]”, „[2,2] [2,2]”, „[1,3] [2,2]”, „[1, 2] [1,3] "," [1,1] [1,2] "] co jest złe
J42161217,
Wygląda na to, że wszystkie rozwiązania są błędne
GB
@Jenny_mathy Jak to jest źle? 3 + 3 to 6, a 2 * 3 to 6. Czy źle rozumiem pytanie.
maple_shaft
brakuje Ci „vice versa”
J42161217,
@Jenny_mathy Głupi błąd z mojej strony! Naprawiłem to, powinno działać teraz.
maple_shaft
0

Aksjomat, n = 83 w 59 sekund tutaj

-- copy the below text in the file name "thisfile.input" 
-- and give something as the command below in the Axiom window:
-- )read C:\Users\thisuser\thisdirectory\thisfile

)cl all
)time on

-- controlla che l'array a e' formato da elementi  a.i<=a.(i+1)
tv(a:List PI):Boolean==(for i in 1..#a-1 repeat if a.i> a.(i+1) then return false;true)

-- funzione incremento: incrementa a, con #a=n=b/3,sotto la regola di "reduce(+,a)+#a-1>=reduce(*,a)"
-- e che n<reduce(*,a)<3*n ed reduce(+,a)<3*n 
inc3(a:List PI):INT==
   i:=1; n:=#a; b:=3*n
   repeat
      if i>n  then return 0
      x:=reduce(*,a)
      if x>=b then a.i:=1
      else
          y:=reduce(+,a)
          if y>b then a.i=1
          else if y+n-1>=x then
                      x:=x quo a.i
                      a.i:=a.i+1
                      x:=x*a.i
                      if tv(a) then break
                      else a.i:=1
          else a.i:=1
      i:=i+1
   if x<=n then return inc3(a) -- x<=n non va
   x

-- ritorna una lista di liste di 4 divisori di n
-- tali che il loro prodotto e' n
g4(n:PI):List List PI==
  a:=divisors(n)
  r:List List PI:=[]
  for i in 1..#a repeat
     for j in i..#a repeat
        x:=a.i*a.j
        if x*a.j>n then break
        for k in j..#a repeat
            y:=x*a.k
            if y*a.k>n then break
            for h in k..#a repeat
                z:=y*a.h
                if z=n  then r:=cons([a.h,a.k,a.j,a.i],r)
                if z>=n then break 
  r

-- ritorna una lista di liste di 3 divisori di n
-- tali che il loro prodotto e' n
g(n:PI):List List PI==
  a:=divisors(n)
  r:List List PI:=[]
  for i in 1..#a repeat
     for j in i..#a repeat
        x:=a.i*a.j
        if x*a.j>n then break
        for k in j..#a repeat
            y:=x*a.k
            if y=n  then r:=cons([a.k,a.j,a.i],r)
            if y>=n then break
  r

-- cerca che [a,b] nn si trovi gia' in r
searchr(r:List List List PI,a:List PI,b:List PI):Boolean==
  aa:=sort(a); bb:=sort(b)
  for i in 1..#r repeat
      x:=sort(r.i.1);y:=sort(r.i.2)
      if x=aa and y=bb then return false
      if x=bb and y=aa then return false
  true

-- input n:PI
-- ritorna r, tale che se [a,b] in r
-- allora #a=#b=n
--        ed reduce(+,a)=reduce(*,b) ed reduce(+,b)=reduce(*,a)
f(n:PI):List List List PI==
  n>100000 or n<=1 =>[]
  a:List PI:=[]; b:List PI:=[]; r:List List List PI:=[]
  for i in 1..n repeat(a:=cons(1,a);b:=cons(1,b))
  if n~=72 and n<86 then  m:=min(3,n)
  else                    m:=min(4,n) 
  q:=reduce(*,a) 
  repeat
    w:=reduce(+,a)
    if n~=72 and n<86 then x:= g(w)
    else                   x:=g4(w)
    if q=w then r:=cons([copy a, copy a],r)
    for i in 1..#x repeat
           for j in 1..m repeat
                  b.j:=(x.i).j
           -- per costruzione abbiamo che reduce(+,a)= prodotto dei b.i=reduce(*,b)
           -- manca solo di controllare che reduce(+,b)=reduce(*,a)=q
           if reduce(+,b)=q and searchr(r,a,b) then r:=cons([copy a, copy b],r)
    q:=inc3(a)
    if q=0 then break
  r

wyniki:

 for i in 2..83 repeat output [i, # f(i)]
   [2,2][3,4][4,3][5,5][6,4][7,6][8,5][9,7][10,7][11,8][12,6][13,10][14,7][15,7]
   [16,10][17,10][18,9][19,12][20,7][21,13][22,9][23,14][24,7][25,13][26,11]
   [27,10][28,11][29,15][30,9][31,16][32,11][33,17][34,9][35,9][36,13][37,19]
   [38,11][39,14][40,12][41,17][42,11][43,20][44,12][45,16][46,14][47,14][48,13]
   [49,16][50,14][51,17][52,11][53,20][54,15][55,17]
   [56,14][57,20][58,17][59,16][60,15][61,28][62,15][63,16][64,17][65,18]
   [66,14][67,23][68,20][69,19][70,13][71,18][72,15][73,30][74,15][75,17][76,18]
   [77,25][78,16][79,27][80,9][81,23][82,17][83,26]


 f 3
    [[[1,2,5],[8,1,1]],[[1,3,3],[7,1,1]],[[1,2,3],[1,2,3]],[[2,2,2],[6,1,1]]]
                                     Type: List List List PositiveInteger
                                   Time: 0.07 (IN) + 0.05 (OT) = 0.12 sec

Sposób uruchamiania nad tekstem w Axiomie polega na skopiowaniu całego tekstu do pliku, zapisaniu pliku o nazwie: Name.input, w oknie Axiom użyj „) read absolutepath / Name”.
wyniki: (# f (i) znajduje długość tablicy f (i), czyli liczbę rozwiązań)

RosLuP
źródło