29

Dato un numero intero positivo n, semplifica la radice quadrata √nnella forma a√bestraendo tutti i fattori quadrati. Gli output a,bdovrebbero essere numeri interi positivi n = a^2 * bcon il bpiù piccolo possibile.

È possibile produrre ae bin qualsiasi ordine in qualsiasi formato ragionevole. Non è possibile omettere gli output 1come impliciti.

Le uscite per n=1..36come (a,b):

1 (1, 1)
2 (1, 2)
3 (1, 3)
4 (2, 1)
5 (1, 5)
6 (1, 6)
7 (1, 7)
8 (2, 2)
9 (3, 1)
10 (1, 10)
11 (1, 11)
12 (2, 3)
13 (1, 13)
14 (1, 14)
15 (1, 15)
16 (4, 1)
17 (1, 17)
18 (3, 2)
19 (1, 19)
20 (2, 5)
21 (1, 21)
22 (1, 22)
23 (1, 23)
24 (2, 6)
25 (5, 1)
26 (1, 26)
27 (3, 3)
28 (2, 7)
29 (1, 29)
30 (1, 30)
31 (1, 31)
32 (4, 2)
33 (1, 33)
34 (1, 34)
35 (1, 35)
36 (6, 1)

Questi sono OEIS A000188 e A007913 .

Correlati: una versione più complessa .

Mostra snippet di codice

var QUESTION_ID=83814,OVERRIDE_USER=20260;function answersUrl(e){return"https://api.stackexchange.com/2.2/questions/83814/answers?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+ANSWER_FILTER}function commentUrl(e,s){return"https://api.stackexchange.com/2.2/answers/"+s.join(";")+"/comments?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+COMMENT_FILTER}function getAnswers(){jQuery.ajax({url:answersUrl(answer_page++),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){answers.push.apply(answers,e.items),answers_hash=[],answer_ids=[],e.items.forEach(function(e){e.comments=[];var s=+e.share_link.match(/\d+/);answer_ids.push(s),answers_hash[s]=e}),e.has_more||(more_answers=!1),comment_page=1,getComments()}})}function getComments(){jQuery.ajax({url:commentUrl(comment_page++,answer_ids),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){e.items.forEach(function(e){e.owner.user_id===OVERRIDE_USER&&answers_hash[e.post_id].comments.push(e)}),e.has_more?getComments():more_answers?getAnswers():process()}})}function getAuthorName(e){return e.owner.display_name}function process(){var e=[];answers.forEach(function(s){var r=s.body;s.comments.forEach(function(e){OVERRIDE_REG.test(e.body)&&(r="<h1>"+e.body.replace(OVERRIDE_REG,"")+"</h1>")});var a=r.match(SCORE_REG);a&&e.push({user:getAuthorName(s),size:+a[2],language:a[1],link:s.share_link})}),e.sort(function(e,s){var r=e.size,a=s.size;return r-a});var s={},r=1,a=null,n=1;e.forEach(function(e){e.size!=a&&(n=r),a=e.size,++r;var t=jQuery("#answer-template").html();t=t.replace("{{PLACE}}",n+".").replace("{{NAME}}",e.user).replace("{{LANGUAGE}}",e.language).replace("{{SIZE}}",e.size).replace("{{LINK}}",e.link),t=jQuery(t),jQuery("#answers").append(t);var o=e.language;/<a/.test(o)&&(o=jQuery(o).text()),s[o]=s[o]||{lang:e.language,user:e.user,size:e.size,link:e.link}});var t=[];for(var o in s)s.hasOwnProperty(o)&&t.push(s[o]);t.sort(function(e,s){return e.lang>s.lang?1:e.lang<s.lang?-1:0});for(var c=0;c<t.length;++c){var i=jQuery("#language-template").html(),o=t[c];i=i.replace("{{LANGUAGE}}",o.lang).replace("{{NAME}}",o.user).replace("{{SIZE}}",o.size).replace("{{LINK}}",o.link),i=jQuery(i),jQuery("#languages").append(i)}}var ANSWER_FILTER="!t)IWYnsLAZle2tQ3KqrVveCRJfxcRLe",COMMENT_FILTER="!)Q2B_A2kjfAiU78X(md6BoYk",answers=[],answers_hash,answer_ids,answer_page=1,more_answers=!0,comment_page;getAnswers();var SCORE_REG=/<h\d>\s*([^\n,]*[^\s,]),.*?(\d+)(?=[^\n\d<>]*(?:<(?:s>[^\n<>]*<\/s>|[^\n<>]+>)[^\n\d<>]*)*<\/h\d>)/,OVERRIDE_REG=/^Override\s*header:\s*/i;

body{text-align:left!important}#answer-list,#language-list{padding:10px;width:290px;float:left}table thead{font-weight:700}table td{padding:5px}

<script src="https://ajax.googleapis.com/ajax/libs/jquery/2.1.1/jquery.min.js"></script> <link rel="stylesheet" type="text/css" href="//cdn.sstatic.net/codegolf/all.css?v=83c949450c8b"> <div id="answer-list"> <h2>Leaderboard</h2> <table class="answer-list"> <thead> <tr><td></td><td>Author</td><td>Language</td><td>Size</td></tr></thead> <tbody id="answers"> </tbody> </table> </div><div id="language-list"> <h2>Winners by Language</h2> <table class="language-list"> <thead> <tr><td>Language</td><td>User</td><td>Score</td></tr></thead> <tbody id="languages"> </tbody> </table> </div><table style="display: none"> <tbody id="answer-template"> <tr><td>{{PLACE}}</td><td>{{NAME}}</td><td>{{LANGUAGE}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table> <table style="display: none"> <tbody id="language-template"> <tr><td>{{LANGUAGE}}</td><td>{{NAME}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table>

Espandi frammento

code-golf math arithmetic

— XNOR
fonte

Lo abbiamo già fatto in precedenza , e questo è stato chiuso come duplicato della sfida qui collegata.

— Flawr,

13

Gelatina , 9 byte

ÆE;0d2ZÆẸ

Provalo online! o verifica tutti i casi di test .

Come funziona

ÆE;0d2ZÆẸ  Main link. Argument: n

ÆE         Exponents; generate the exponents of n's prime factorization.
  ;0       Append 0 since 1ÆE returns [].
    d2     Divmod by 2.
      Z    Zip/transpose to group quotients and remainders.
       ÆẸ  Unexponent; turn the exponents of prime factorizations into integers.

— Dennis
fonte

3

In UTF-8 lo è, ma Jelly usa una tabella codici personalizzata. Il collegamento byte nell'intestazione punta ad esso.

— Dennis,

Si registra quel commento molto, quindi forse si dovrebbe fare i byte come più chiara (ad esempio: [bytes](link-to-byes) (not UTF-8).

— NoOneIsHere

12

PARI / GP, 12 byte

n->core(n,1)

corerestituisce la parte squarefree di ndefault, ma impostando il flag del secondo argomento su 1 fa sì che restituisca entrambe le parti. L'ordine di uscita è (b, a), ad es (n->core(n,1))(12) -> [3, 2].

— SP3000
fonte

11

Python 2, 43 byte

k=n=input()
while n%k**2:k-=1
print k,n/k/k

Provalo su Ideone .

— Dennis
fonte

6

MATL , 12 byte

t:U\~f0)GyU/

Provalo online!

Spiegazione

t     % Take input n implicitly. Duplicate
:U    % Push [1 4 9 ... n^2]
\~    % True for entries that divide the input
f0)   % Get (1-based) index of the last (i.e. largest) dividing number
G     % Push input again
y     % Duplicate index of largest dividing number
U     % Square to recover largest dividing number
/     % Divide input by that. Implicitly display stack

— Luis Mendo
fonte

4

Julia, 32 byte

\(n,k=n)=n%k^2>0?n\~-k:[k,n/k^2]

Provalo online!

— Dennis
fonte

2

Mathematica 34 byte

#/.{a_ b_^_:>{a, b},_[b_,_]:>{1,b}}&

Questo dice di sostituire tutto l'input ( #) secondo le seguenti regole: (1) un numero, a , volte la radice quadrata di b dovrebbe essere sostituita da {a, b} e una funzione b alla potenza di qualunque cosa dovrebbe essere sostituita da {1, b }. Si noti che la funzione presuppone che l'input sia nella forma Sqrt[n],. Non funzionerà con altri tipi di input.

Questa funzione senza nome è insolitamente criptica per Mathematica. Può essere chiarito in qualche modo mostrando la sua forma completa, seguita da sostituzioni delle forme più brevi originali.

Function[
   ReplaceAll[
      Slot[1],
      List[
         RuleDelayed[Times[Pattern[a,Blank[]],Power[Pattern[b,Blank[]],Blank[]]],List[a,b]],
         RuleDelayed[Blank[][Pattern[b,Blank[]],Blank[]],List[1,b]]]]]

che è lo stesso di

   ReplaceAll[
      #,
      List[
         RuleDelayed[Times[Pattern[a,Blank[]],Power[Pattern[b,Blank[]],Blank[]]],List[a,b]],
         RuleDelayed[Blank[][Pattern[b,Blank[]],Blank[]],List[1,b]]]]&

e

ReplaceAll[#, 
  List[RuleDelayed[
    Times[Pattern[a, Blank[]], 
     Power[Pattern[b, Blank[]], Blank[]]], {a, b}], 
   RuleDelayed[Blank[][Pattern[b, Blank[]], Blank[]], {1, b}]]] &

e

ReplaceAll[#, 
  List[RuleDelayed[Times[a_, Power[b_, _]], {a, b}], 
   RuleDelayed[Blank[][b_, _], {1, b}]]] &

e

ReplaceAll[#, {RuleDelayed[a_*b^_, {a, b}], RuleDelayed[_[b_, _], {1, b}]}]&

e

ReplaceAll[#, {a_*b^_ :> {a, b}, _[b_, _] :> {1, b}}] &

— DavidC
fonte

1

Pyth, 15 byte

/Q^
ef!%Q^T2SQ2

Suite di test.

— Perdente Suora
fonte

1

Matlab, 51 byte

x=input('');y=1:x;z=y(~rem(x,y.^2));a=z(end)
x/a^2

Spiegazione

x=input('')       -- takes input
y=1:x             -- numbers from 1 to x
z=y(~rem(x,y.^2)) -- numbers such that their squares divide x
a=z(end)          -- biggest such number (first part of output)
x/a^2             -- remaining part

— Pajonk
fonte

1

JavaScript (ECMAScript 2016), 40 byte

n=>{for(k=n;n%k**2;k--);return[k,n/k/k]}

Fondamentalmente una porta JavaScript della risposta Python 2 di Dennis .

Provalo su JSBin .

Nota: non funziona in modalità rigorosa, perché knon è inizializzato da nessuna parte. Per farlo funzionare in modalità rigorosa, k=nnel ciclo dovrebbe essere cambiato in let k=n.

— Michał Perłakowski
fonte

1

Haskell, 43> 42 byte

Soluzione di forza bruta.

Salvato 1 byte grazie a Xnor

f n=[(x,y)|y<-[1..],x<-[1..n],x*x*y==n]!!0

— Damien
fonte

Bella soluzione, mi piace come non lo usa mod o div. Penso che tu possa fare a y<-[1..]causa della pigrizia.

— xnor

Si hai ragione. Non è stato possibile con la mia prima soluzione, last[(x,y)|x<-[1..n],y<-[1..n],x*x*y==n]ma ora funzionerà. Grazie. Hai la tua soluzione in Haskell?

— Damien,

1

05AB1E, 14 byte

Lv¹ynÖi¹yn/y‚ï

spiegato

Lv              # for each x in range(1,N) inclusive
  ¹ynÖi         # if N % x^2 == 0
       ¹yn/y‚ï  # create [N/x^2,x] pairs, N=12 -> [12,1] [3,2]
                # implicitly output last found pair

Provalo online

— Emigna
fonte

1

Python, 74 byte

def e(k):a=filter(lambda x:k/x**2*x*x==k,range(k,0,-1))[0];return a,k/a**2

Abbastanza semplice.

— Backerupper
fonte

0

Python 2.7 (non golfizzato) - 181 byte

def e(n):   
 for x in range(1,n+1):
  r=(1,x)
  for i in range(1,x+1):
   l=i**.5
   for j in range(1,x+1): 
    if i*j==x and l%1==0 and j<r[1]:r=(int(l),j)                
  print x,r

Esegui come: e (numero) ad es. e (24)

Uscita campione:

>> e(24)
1 (1, 1)
2 (1, 2)
3 (1, 3)
4 (2, 1)
5 (1, 5)
6 (1, 6)
7 (1, 7)
8 (2, 2)
9 (3, 1)
10 (1, 10)
11 (1, 11)
12 (2, 3)
13 (1, 13)
14 (1, 14)
15 (1, 15)
16 (4, 1)
17 (1, 17)
18 (3, 2)
19 (1, 19)
20 (2, 5)
21 (1, 21)
22 (1, 22)
23 (1, 23)
24 (2, 6)

— Swadhikar C
fonte

1

Per favore,

— prova a golfare la

0

APL, 25 caratteri

 {(⊢,⍵÷×⍨)1+⍵-0⍳⍨⌽⍵|⍨×⍨⍳⍵}

In inglese:

0⍳⍨⌽⍵|⍨×⍨⍳⍵: indice del più grande dei quadrati fino ad n che divide completamente n;
1+⍵-: l'indice si trova nell'array invertito, quindi regola l'indice
(⊢,⍵÷×⍨): produce il risultato: l'indice stesso (a) e il quoziente b (cioè n ÷ a * a)

Test:

     ↑{(⊢,⍵÷×⍨)⊃z/⍨0=⍵|⍨×⍨z←⌽⍳⍵}¨⍳36
1  1
1  2
1  3
2  1
1  5
1  6
1  7
2  2
3  1
1 10
1 11
2  3
1 13
1 14
1 15
4  1
1 17
3  2
1 19
2  5
1 21
1 22
1 23
2  6
5  1
1 26
3  3
2  7
1 29
1 30
1 31
4  2
1 33
1 34
1 35
6  1

— lstefano
fonte

0

JavaScript (ECMAScript 6), 35 byte

f=(n,k=n)=>n/k%k?f(n,--k):[k,n/k/k]

JavaScript 1+, 37 B

for(k=n=prompt();n/k%k;--k);[k,n/k/k]

— l4m2
fonte

0

J , 19 byte

(0|:0 2#:])&.(_&q:)

Provalo online!

Come la soluzione Jelly .

— FrownyFrog
fonte

Semplifica una radice quadrata

Gelatina , 9 byte

Come funziona

PARI / GP, 12 byte

Python 2, 43 byte

MATL , 12 byte

Spiegazione

Julia, 32 byte

Mathematica 34 byte

Pyth, 15 byte

Matlab, 51 byte

JavaScript (ECMAScript 2016), 40 byte

Haskell, 43> 42 byte

05AB1E, 14 byte

Python, 74 byte

APL, 25 caratteri

JavaScript (ECMAScript 6), 35 byte

JavaScript 1+, 37 B

J , 19 byte