Це номер Сміта?


28

Опис виклику

Число Сміт є композитом числа, сума цифр дорівнює сумі сум цифр його простих множників. Давши ціле число N, визначте, чи це номер Сміта чи ні.

Перші кілька номерів Сміта є 4, 22, 27, 58, 85, 94, 121, 166, 202, 265, 274, 319, 346, 355, 378, 382, 391, 438(послідовність A006753 в OEIS).

Зразок вводу / виводу

18: False (sum of digits: 1 + 8 = 9; factors: 2, 3, 3; sum of digits of factors: 2 + 3 + 3 = 8)
22: True
13: False (meets the digit requirement, but is prime)
666: True (sum of digits: 6 + 6 + 6 = 18; factors: 2, 3, 3, 37; sum of digits of factors: 2 + 3 + 3 + 3 + 7 = 18)
-265: False (negative numbers can't be composite)
0: False (not composite)
1: False (not composite)
4937775: True

Примітки

  • Ваш код може бути функцією (методом) або повноцінною робочою програмою,
  • Замість таких слів, як Trueі False, ви можете використовувати будь-які неправдиві та хибні значення, якщо зрозуміло, що вони є,
  • Це завдання з , тому зробіть свій код якомога коротшим!

6
Мені довелося прочитати це: "сума цифр дорівнює сумі цифр її основних факторів" кілька разів: P
Стюі Гріффін

@StewieGriffin: Так, це досить складне речення, але я відчував, що мені потрібно дати правильне визначення, а не покладатися лише на приклади :)
shooqie

2
Це одне з тих питань, де я думав, що "Java + це = ні", але я прихилився до цієї ідеї: P
Shaun Wild

3
Я іноді помічаю шаблони в цифрах, сумі цифр тощо, але чи справді люди помічають подібні речі: "Альберт Віланський придумав термін Сміт номер, коли помітив визначальну властивість у номері телефону свого зятя" ?
Стюі Гріффін

1
@StewieGriffin: Так, це як Рамануджан і 1729 рік, завжди також здивував мене.
shooqie

Відповіді:


9

Желе , 12 11 байт

Æfḟȯ.DFżDSE

Повертає 1 для номерів Сміта і 0 в іншому випадку. Спробуйте в Інтернеті! або перевірити всі тестові випадки .

Фон

Æf(проста факторизація) та D(ціле число до десяткових) реалізуються так, що P(добуток) та (десятковий до цілого числа) становлять ліві оберти.

Для цілих чисел -4 до 4 , Æfповертає наступне.

-4 -> [-1, 2, 2]
-3 -> [-1, 3]
-2 -> [-1, 2]
-1 -> [-1]
 0 -> [0]
 1 -> []
 2 -> [2]
 3 -> [3]
 4 -> [2, 2]

Для чисел -10, -1, -0.5, 0, 0,5, 1, 10 , Dповертає наступне.

-11   -> [-1, -1]
-10   -> [-1, 0]
 -1   -> [-1]
 -0.5 -> [-0.5]
  0   -> [0]
  0.5 -> [0.5]
  1   -> [1]
 10   -> [1, 0]
 11   -> [1, 1]

Як це працює

Æfḟȯ.DFżDSE  Main link. Argument: n (integer)

Æf           Yield the prime factorization of n.
  ḟ          Filter; remove n from its prime factorization.
             This yields an empty array if n is -1, 0, 1, or prime.
   ȯ.        If the previous result is an empty array, replace it with 0.5.
     D       Convert all prime factors to decimal.
      F      Flatten the result.
        D    Yield n in decimal.
       ż     Zip the results to both sides, creating a two-column array.
         S   Compute the sum of each column.
             If n is -1, 0, 1, or prime, the sum of the prime factorization's
             digits will be 0.5, and thus unequal to the sum of the decimal array.
             If n < -1, the sum of the prime factorization's digits will be
             positive, while the sum of the decimal array will be negative.
          E  Test both sums for equality.

2
Це серйозно класне рішення, про яке я повинен сказати!
Емінья

@Emigna - Це те, що я зробив, але реалізував набагато перевершив: D
Джонатан Аллан

@JonathanAllan На жаль, я не розмовляю з Jelly, тому я не маю поняття, що робить ваш код :)
Emigna

1
@Emigna - так, я планував розробити, як пограти в гольф, перш ніж додати розділ, як це працює.
Джонатан Аллан

9

Python 2, 122 115 110 106 байт

n=m=input()
s=0
for d in range(2,n):
 while n%d<1:n/=d;s+=sum(map(int,`d`))
print n<m>s==sum(map(int,`m`))

Збережено 4 байти завдяки Деннісу

Спробуйте це на ideone.com

Пояснення

Читає число на stdin і виводить, Trueякщо число є номером Сміта або Falseякщо його немає.

n=m=input()                  # stores the number to be checked in n and in m
s=0                          # initializes s, the sum of the sums of digits of prime factors, to 0
for d in range(2,n):         # checks all numbers from 2 to n for prime factors
 while n%d<1:                # while n is divisible by d
                             #   (to include the same prime factor more than once)
  n/=d                       # divide n by d
  s+=sum(map(int,`d`))       # add the sum of the digits of d to s
print                        # print the result: "True" if and only if
      n<m                    #   n was divided at least once, i.e. n is not prime
      >                      #   and m>s (always true) and
      s==sum(map(int,`m`))   #   s is equal to the sum of digits of m (the input)

1
Вниз виборця - можливо, буде корисно додати коментар, щоб пояснити, чому
Джонатан Аллан

6
@JonathanAllan Режим голосування був автоматично переданий користувачем спільноти, коли відповідь була відредагована. Я вважаю це помилкою .
Денніс

1
Останній рядок можна переписати як print n<m>s==sum(map(int,`m`)).
Денніс

@Dennis Це чудове використання ланцюгового порівняння!
LevitatingLion

8

Брахілог , 19 байт

@e+S,?$pPl>1,P@ec+S

Спробуйте в Інтернеті!

Пояснення

@e+S,                 S is the sum of the digits of the input.
     ?$pP             P is the list of prime factors of the input.
        Pl>1,         There are more than 1 prime factors.
             P@e      Split each prime factor into a list of digits.
                c     Flatten the list.
                 +S   The sum of this list of digits must be S.

2
@JonathanAllan Це так . У Брахілогу від’ємний знак для чисел - це _(так званий низький мінус ).
Фаталізувати


5

PowerShell v3 +, 183 байти

param($n)$b=@();for($a=$n;$a-gt1){2..$a|?{'1'*$_-match'^(?!(..+)\1+$)..'-and!($a%$_)}|%{$b+=$_;$a/=$_}}$n-notin$b-and(([char[]]"$n")-join'+'|iex)-eq(($b|%{[char[]]"$_"})-join'+'|iex)

Немає вбудованої простої перевірки. Немає вбудованого факторингу. Немає вбудованої цифри-суми. Все зроблено рукою. : D

Приймає введення $nяк ціле число, встановлює $bрівний порожній масив. Ось $bнаша колекція основних факторів.

Далі - forпетля. Спочатку встановлюємо $aрівне нашому вхідному числу, а умовне значення поки $aне дорівнює або дорівнює 1. Цей цикл збирається знайти основні фактори.

Ми обертаємо цикл 2до $a, використовуємо Where-Object( |?{...}), щоб витягнути праймери , які також є чинниками !($a%$_). Вони подаються у внутрішню петлю, |%{...}яка розміщує фактор на $bта ділиться $a(таким чином ми врешті-решт дістанемося 1).

Отже, зараз у нас є всі наші основні фактори $b. Час сформулювати наш булевий вихід. Нам необхідно переконатися в тому , що $nє -notin $b, тому що , якщо це те , що означає , що $nпростий, і тому не є число Сміт. Крім того, ( -and) нам потрібно переконатися, що наші два набори розрядних сум є -eqнереальними. Булевий залишок залишається на конвеєрі і вихід неявний.

Примітка. Потрібна версія v3 або новіша для -notinоператора. Я все ще виконую введення для 4937775(це повільно обчислювати), тому я оновлю це, коли це закінчиться. Через 3+ годин у мене з’явилася помилка stackoverflow. Отже, десь є верхня межа. Ну добре.

Це буде спрацьовувати при негативному введенні, нулі чи одиниці, тому що права рука -andзафіксує помилку, намагаючись обчислити цифри суми (показано нижче), що призведе до того, що половина піде на $falseоцінку. Оскільки STDERR за замовчуванням ігнорується , а правильний вихід все ще відображається, це добре.


Тестові справи

PS C:\Tools\Scripts\golfing> 4,22,27,58,85,94,18,13,666,-265,0,1|%{"$_ -> "+(.\is-this-a-smith-number.ps1 $_)}
4 -> True
22 -> True
27 -> True
58 -> True
85 -> True
94 -> True
18 -> False
13 -> False
666 -> True
Invoke-Expression : Cannot bind argument to parameter 'Command' because it is an empty string.
At C:\Tools\Scripts\golfing\is-this-a-smith-number.ps1:1 char:179
+ ... "$_"})-join'+'|iex)
+                    ~~~
    + CategoryInfo          : InvalidData: (:String) [Invoke-Expression], ParameterBindingValidationException
    + FullyQualifiedErrorId : ParameterArgumentValidationErrorEmptyStringNotAllowed,Microsoft.PowerShell.Commands.InvokeExpressionCommand

-265 -> False
Invoke-Expression : Cannot bind argument to parameter 'Command' because it is an empty string.
At C:\Tools\Scripts\golfing\is-this-a-smith-number.ps1:1 char:179
+ ... "$_"})-join'+'|iex)
+                    ~~~
    + CategoryInfo          : InvalidData: (:String) [Invoke-Expression], ParameterBindingValidationException
    + FullyQualifiedErrorId : ParameterArgumentValidationErrorEmptyStringNotAllowed,Microsoft.PowerShell.Commands.InvokeExpressionCommand

0 -> False
Invoke-Expression : Cannot bind argument to parameter 'Command' because it is an empty string.
At C:\Tools\Scripts\golfing\is-this-a-smith-number.ps1:1 char:179
+ ... "$_"})-join'+'|iex)
+                    ~~~
    + CategoryInfo          : InvalidData: (:String) [Invoke-Expression], ParameterBindingValidationException
    + FullyQualifiedErrorId : ParameterArgumentValidationErrorEmptyStringNotAllowed,Microsoft.PowerShell.Commands.InvokeExpressionCommand

1 -> False


3

Желе , 27 25 23 байт

(подальший гольф, ймовірно, напевно можливий)

ḢDS×
ÆFÇ€SḢ
DS=Ça<2oÆP¬

Повертається 0за хибну чи 1істинну

Усі тестові приклади на TryItOnline

Як?

DS=Ça<2oÆP¬ - main link takes an argument, n
DS          - transform n to a decimal list and sum up
   Ç        - call the previous link (ÆFÇ€SḢ)
  =         - test for equality
     <2     - less than 2?
    a       - logical and
        ÆP  - is prime?
       o    - logical or
          ¬ - not
            - all in all tests if the result of the previous link is equal to the digit
              sum if the number is composite otherwise returns 0.

ÆFÇ€SḢ - link takes an argument, n again
ÆF     - list of list of n's prime factors and their multiplicities
  Ç€   - apply the previous link (ḢDS×) for each
    S  - sum up
     Ḣ - pop head of list (there will only be one item)

ḢDS× - link takes an argument, a factor, multiplicity pair
Ḣ    - pop head, the prime factor - modifies list leaving the multiplicity
 DS  - transform n to a decimal list and sum up
   × - multiply the sum with the multiplicity

3

Власне, 18 байт

На жаль, насправді немає вбудованої факторизації, яка дає множину простих факторів для кратності, тому мені довелося зламати один разом. Пропозиції з гольфу вітаються. Спробуйте в Інтернеті!

;w`i$n`MΣ♂≈Σ@$♂≈Σ=

Ungolfing

         Implicit input n.
;w       Duplicate n and get the prime factorization of a copy of n.
`...`M   Map the following function over the [prime, exponent] lists of w.
  i        Flatten the list. Stack: prime, exponent.
  $n       Push str(prime) to the stack, exponent times.
            The purpose of this function is to get w's prime factors to multiplicity.
Σ        sum() the result of the map.
          On a list of strings, this has the same effect as "".join()
♂≈Σ      Convert every digit to an int and sum().
@        Swap the top two elements, bringing other copy of n to TOS.
$♂≈Σ     Push str(n), convert every digit to an int, and sum().
=        Check if the sum() of n's digits is equal 
          to the sum of the sum of the digits of n's prime factors to multiplicity.
         Implicit return.

3

Haskell, 120 105 байт

1%_=[];a%x|mod a x<1=x:div a x%x|0<1=a%(x+1)
p z=sum[read[c]|c<-show z]
s x|z<-x%2=z<[x]&&sum(p<$>z)==p x

2

Октава, 80 78 байт

t=num2str(factor(x=input('')))-48;disp(any(t<0)&~sum([num2str(x)-48 -t(t>0)]))

Пояснення:

factor(x=input(''))                 % Take input, store as x and factor it
num2str(factor(x=input('')))-48     % Convert it to an array (123 -> [1 2 3]) 
                                    % and store as t
any(t<0)                            % Check if there are several prime factors
                                    % [2 3] -> [2 -16 3]
sum([num2str(x)-48 -t(t>0)])        % Check if sum of prime factor
                                    % is equal the sum of digits

Спробуйте в Інтернеті .


1
Це any(t<0)для не первинності дуже розумно
Луїс Мендо

2

Pyth, 21 байт

&&>Q1!P_QqsjQTssmjdTP

Програма, яка приймає введення цілого числа і друкує Trueабо Falseяк доречно.

Спробуйте в Інтернеті

Як це працює

&&>Q1!P_QqsjQTssmjdTP  Program. Input: Q
           jQT         Yield digits of the base-10 representation of Q as a list
          s            Add the digits
                    P  Yield prime factors of Q (implicit input fill)
                mjdT   Map base-10 representation across the above, yielding digits of each
                       factor as a list of lists
               s       Flatten the above
              s        Add up the digits
         q             Those two sums are equal
&                      and
  >Q1                  Q>1
 &                     and
     !P_Q              Q is not prime
                       Implicitly print

2

Perl 6 , 92 88 87 байт

{sub f(\i){my \n=first i%%*,2..i-1;n??n~f i/n!!i}
!.is-prime&&$_>1&&.comb.sum==.&f.comb.sum}

{sub f(\i){my \n=first i%%*,2..^i;n??[n,|f i/n]!!|i}
$_>.&f>1&&.comb.sum==.&f.comb.sum}

Анонімна функція, яка повертає Bool.

  • Тепер проводиться 100% ручна факторизація та перевірка первинності.
  • Збережено кілька байтів, випробувавши як "вхід> 1", так і "кількість факторів> 1" з одним ланцюговим порівнянням, оскільки m> Ω (m) .

( спробуйте в Інтернеті )

EDIT: -1 байт завдяки b2gills


2..i-1краще написано як 2..^i.
Бред Гілберт b2gills

2

Java 7, 509 506 435 426 419 230 байт

boolean c(int n){return n<2|p(n)?0>1:d(n)==f(n);}int d(int n){return n>9?n%10+d(n/10):n;}int f(int n){int r=0,i;for(i=1;++i<=n;)for(;n%i<1;n/=i,r+=i>9?d(i):i);return r;}boolean p(int n){int i=2;while(i<n)n=n%i++<1?0:n;return n>1;}

Я повинен був слухати коментар @BasicallyAlanTuring ..

Це одне з тих запитань, де я думав, що "Java + це = ні", проте я виступав за цю ідею: P

Ну добре .. Деякі мови програмування використовують один байт для прайме-факторів або прайм-перевірки, але Java, безумовно, не з них.

EDIT: Вдвічі зменшив кількість байтів, коли я трохи задумався.

Безлічі (розбір ..) та тестові справи:

Спробуйте тут.

class M{
  static boolean c(int n){
    return n < 2 | p(n)
            ? 0 > 1 //false
            : d(n) == f(n);
  }

  // Sums digits of int
  static int d(int n) {
    return n > 9
            ? n%10 + d(n/10)
            : n;
  }

  // Convert int to sum of prime-factors
  static int f(int n) {
    int r = 0,
        i;
    for(i = 1; ++i <= n; ){
      for( ; n % i < 1; n /= i,
                        r += i > 9 ? d(i) : i);
    }
    return r;
  }

  // Checks if the int is a prime
  static boolean p(int n){
    int i = 2;
    while(i < n){
      n = n % i++ < 1
           ? 0
           : n;
    }
    return n > 1;
  }

  public static void main(String[] a){
    System.out.println(c(18));
    System.out.println(c(22));
    System.out.println(c(13));
    System.out.println(c(666));
    System.out.println(c(-256));
    System.out.println(c(0));
    System.out.println(c(1));
    System.out.println(c(4937775));
  }
}

Вихід:

false
true
false
true
false
false
false
true

2

Брахілог (новіший) , 11 байт

¬ṗ&ẹ+.&ḋcẹ+

Спробуйте в Інтернеті!

Присудок досягає успіху, якщо вхід є номером Сміта, а невдалий, якщо його немає.

               The input
¬ṗ             is not prime,
  &            and the input's 
   ẹ           digits
    +          sum to
     .         the output variable,
      &        and the input's 
       ḋ       prime factors' (getting prime factors of a number < 1 fails)
        c      concatenated
         ẹ     digits
          +    sum to
               the output variable.



1

Пайк, 16 байт

Pm[`mbs(sQ[qRlt*

Спробуйте тут!


1
Помилки без результату для введення менше2
Джонатан Аллан

@JonathanAllan жоден вихід до stdout неправдивий. Якщо попередження вимкнено, stderr також ігнорується
Blue

Я знав, що ми можемо ігнорувати stderr, але жоден вихід не здається трохи дивним ... але якщо це прийнятно, то він прийнятний.
Джонатан Аллан

Особисто я не впевнений, чи це прийнятно, але можу сказати, що це правильно?
Блакитний


1

APL (Dyalog Extended) , 36 29 байт SBCS

Ця відповідь завдячує своєю голісністю монаді Extended за повернення основних факторів числа, і це краще при базовому перетворенні, ніж у Dyalog Unicode.

Редагувати: -7 байт завдяки dzaima.

{2>⍵:0⋄(⊃=+/-⊃×2<≢)+⌿10⊤⍵,⍭⍵}

Спробуйте в Інтернеті!

Пояснення

{1⋄(3)2}  A dfn, a function in brackets.  is a statement separator.
          The numbers signify the sections in the order they are explained.

2>⍵:0  If we have a number less than 2,
       we immediately return 0 to avoid a DOMAIN ERROR.

+⌿10⊤⍵,⍭⍵
        ⍭⍵  We take the factors of ⍵, our input as our right argument,
      ⍵,    and append it to our input again.
   10      before converting the input and its factors into a matrix of their base-10 digits
            (each row is the places, units, tens, hundreds, etc.)
+⌿         And taking their sum across the columns of the resulting matrix,
            to give us the sum of their digits, their digit-sums.

(⊃=+/-⊃×2<≢)  We run this section over the list of sums of digits above.
 ⊃=+/-⊃       We check if the main digit-sum (of our input)
               Is equal to the sum of our digit-sums
               (minus our main digit-sum that is also still in the list)
        ×2<≢   The trick here is that we can sneak in our composite check
               (if our input is prime there will be only two numbers, 
               the digit-sum of the prime,
               and the digit-sum of its sole prime factor, itself)
               So if we have a prime, we zero our (minus our main sum)
               in the calculation above, so that primes will not succeed in the check.
               We return the result of the check.

29 байт -{2>⍵:0⋄(⊃=+/-⊃×2<≢)+⌿10⊤⍵,⍭⍵}
dzaima


1

C (gcc) , 139 136 байт

S(m,i,t,h,_){t=m=m<2?2:m;for(_=h=i=1;m>1;h=1){while(m%++h);for(m/=h;i+=h%10,h/=10;);}while(t%++h);for(m=t;_+=m%10,m/=10;);m=t-h?i==_:0;}

Спробуйте в Інтернеті!

-3 байти завдяки стельовій коті

Пояснення:

/* 
 * Variable mappings:
 *  is_smith      => S
 *  argument      => m
 *  factor_digits => i
 *  arg_copy      => t
 *  least_factor  => h
 *  digit_sum     => _    
 */
int is_smith(int argument){                     /* S(m,i,t,h,_){ */
    int factor_digits;
    int arg_copy;
    int least_factor;
    int digit_sum;

    /* 
     * The cases of 0 and 1 are degenerate. 
     * Mapping them to a non-degenerate case with the right result.
     */
    if (argument < 2) {                         /* t=m=m<2?2:m; */
        argument = 2;
    }
    arg_copy = argument;

    /* 
     * Initializing these to 1 instead of zero is done for golf reasons.
     * In the end we just compare them, so it doesn't really matter.
     */
    factor_digits = 1;                          /* for(_=h=i=1; */
    digit_sum = 1;

    /* Loop over each prime factor of argument */
    while (argument > 1) {                      /* m>1; */

        /*
         * Find the smallest factor 
         * Note that it is initialized to 1 in the golfed version since prefix
         * increment is used in the modulus operation.
         */
        least_factor = 2;                       /* h=1){ */
        while (argument % least_factor != 0)    /* while(m% */
            least_factor++;                     /* ++h); */
        argument /= least_factor;               /* for(m/=h; */

        /* Add its digit sum to factor_digits */
        while (least_factor > 0) {
            factor_digits += least_factor % 10; /* i+=h%10, */
            least_factor /= 10;                 /* h/=10;) */
        }                                       /* ; */

    }                                           /* } */

    /* In the golfed version we get this for free in the for loop. */
    least_factor = 2;
    while (arg_copy % least_factor != 0)        /* while(t% */
        least_factor++;                         /* ++h); */

    /* Restore the argument */
    argument = arg_copy;                        /* for(m=t; */

    /* Compute the arguments digit sum */
    while (argument > 0) {
        digit_sum += argument % 10;             /* _+=m%10, */
        argument /= 10;                         /* m/=10;) */
    }                                           /* ; */

    /* This return is done by assigning to first argument when golfed. */
                                                /* m= */
    if (arg_copy == least_factor) {             /* t==h? */
        return 0; /* prime input */             /* 0 */
    } else {                                    /* : */
        return digit_sum == factor_digits;      /* i == _ */
    }                                           /* ; */
}                                               /* } */

Це вніс кілька помилок (наприклад, 2 і 3), але я думаю, що це все-таки має бути досяжним.
LambdaBeta

Запропонувати t-h&&i==_замістьt-h?i==_:0
roofcat

0

Ракетка 176 байт

(define(sd x)(if(= x 0)0(+(modulo x 10)(sd(/(- x(modulo x 10))10)))))
(require math)(define(f N)
(if(=(for/sum((i(factorize N)))(*(sd(list-ref i 0))(list-ref i 1)))(sd N))1 0))

Повертає 1, якщо істинно, і 0, якщо невірно:

(f 27)
1
(f 28)
0
(f 85)
1
(f 86)
0

Детальна версія:

(define (sd x)   ; fn to find sum of digits
  (if (= x 0)
      0
      (+ (modulo x 10)
         (sd (/ (- x (modulo x 10)) 10)))))

(require math)
(define (f N)
  (if (= (for/sum ((i (factorize N)))
           (* (sd (list-ref i 0))
              (list-ref i 1)))
         (sd N)) 1 0))

0

Іржа - 143 байти

fn t(mut n:u32)->bool{let s=|k:u32| (2..=k).fold((0,k),|(a,m),_|(a+m%10,m/10));s(n).0==(2..n).fold(0,|mut a,d|{while n%d<1{n/=d;a+=s(d).0};a})}

запозичене рішення пітона @levitatinglion ... принаймні, це коротше, ніж у Java ...

degolfed на play.rust-lang.org


0

APL (NARS), 33 знака, 66 байт

{1≥≢k←π⍵:0⋄s←{+/⍎¨⍕⍵}⋄(s⍵)=+/s¨k}

"π⍵" коефіцієнти повернення списку ⍵, припустимо, що вхід є одним натуральним числом> = 1; тест:

  h←{1≥≢k←π⍵:0⋄s←{+/⍎¨⍕⍵}⋄(s⍵)=+/s¨k}
  (h¨1..100)/1..100
4 22 27 58 85 94 

0

C (gcc), 177 байт

Визначає функцію, Qяка повертає 0 для номерів Сміта, а ненульова - для номерів Сміта

#define r return
O(D,i){for(i=0;D>0;i+=D%10,D-=D%10,D/=10);r i;}D(O,o){for(o=1;o<O;)if(O%++o<1)r o;r O;}Q(p,q,i,j){if(p^(q=D(i=p))){for(j=0;p>1;q=D(p/=q))j+=O(q);r j^O(i);}r 1;}

Спробуйте в Інтернеті!

Пояснення:

// Return the sum of digits of D if D > 0, otherwise 0
O(D,i){
    // While D is greater than 0:
    // Add the last digit of D to i, and remove the last digit from D
    for(i=0;D>0;i+=D%10,D-=D%10,D/=10);
    return i;
}
// Return the smallest prime factor of O if O>1 else O
D(O,o){
    // Iterate over numbers less than O
    for(o=1;o<O;)
        // If O is divisible by o return o
        if(O%++o<1)
            return o;
    // Otherwise return O
    return O;
}
Q(p,q,i,j){
    // Set q to D(p) and i to p
    // If p != D(p) (i.e, p is composite and > 0)
    if(p^(q=D(i=p))){
        // Iterate over the prime factors of p and store their digit sum in j
        for(j=0;p>1;q=D(p/=q))
            j+=O(q);
        // i is the original value of p. If O(i)^j == 0, O(i) == j
        return j^O(i);
    }
    // If p was composite or < 0, return 1
    return 1;
}


Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.