มีฟังก์ชั่นการสร้างที่สำคัญจำนวนมาก ค่อนข้างทั้งหมดของพวกเขาถูกสร้างขึ้นและอยู่บนพื้นฐานของ Eratosthenes, ฟังก์ชั่นMöbiusหรือทฤษฎีบทของวิลสันและโดยทั่วไปจะไม่สามารถคำนวณได้ในทางปฏิบัติ แต่ยังมีเครื่องกำเนิดไฟฟ้าที่มีโครงสร้างที่ง่ายมากและถูกค้นพบโดยบังเอิญ

ในปี 2003 Stephen Wolfram ได้สำรวจสมการการเกิดซ้ำแบบซ้อนในการทดสอบคอมพิวเตอร์สดที่ NKS Summer School กลุ่มคนรอบ ๆ Matthew Frank ติดตามการทดลองเพิ่มเติมและค้นพบคุณสมบัติที่น่าสนใจของการกลับเป็นซ้ำ

a(n) = a(n-1) + gcd(n,a(n-1))

a(1) = 7ที่มีค่าเริ่มต้นของ ความแตกต่างa(n) - a(n-1) = gcd(n,a(n-1))มักจะเป็น 1 หรือนายก ความแตกต่างแรก ๆ คือ ( OEIS A132199 ):

1, 1, 1, 5, 3, 1, 1, 1, 1, 11, 3, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 23, 3, 1, 1, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 47, 3, 1, 5, 3, ...

หากเราไม่ใช้ 1s เราจะได้รับลำดับต่อไปนี้ ( OEIS A137613 ):

5, 3, 11, 3, 23, 3, 47, 3, 5, 3, 101, 3, 7, 11, 3, 13, 233, 3, 467, 3, 5, 3, 
941, 3, 7, 1889, 3, 3779, 3, 7559, 3, 13, 15131, 3, 53, 3, 7, 30323, 3, ...

Eric S. Rowlandพิสูจน์ความเป็นอันดับหนึ่งของแต่ละองค์ประกอบในรายการนี้ไม่กี่ปีต่อมา อย่างที่คุณเห็นช่วงเวลาต่าง ๆ และบางส่วนก็ปรากฏขึ้นหลายครั้ง มันยังได้รับการพิสูจน์แล้วว่าลำดับนั้นรวมถึงช่วงเวลาที่แตกต่างกันมากมาย ยิ่งไปกว่านั้นมันคือการคาดเดาว่าทุกช่วงเวลาแปลก ๆ ปรากฏขึ้น

เนื่องจากตัวสร้างหลักนี้ไม่ได้ถูกสร้างขึ้น แต่ค้นพบโดยบังเอิญโดยที่ตัวสร้างหลักเรียกว่า "เกิดขึ้นตามธรรมชาติ" แต่ขอให้สังเกตว่าในทางปฏิบัติเครื่องกำเนิดไฟฟ้านี้ยังไม่สามารถคำนวณได้ เมื่อปรากฎออกมาไพรม์ p จะปรากฏขึ้นหลังจาก(p–3)/21s ติดต่อกันเท่านั้น อย่างไรก็ตามการนำเครื่องมือสร้างสำคัญนี้ไปใช้จะเป็นงานของคุณ

ท้าทาย:

เขียนฟังก์ชั่นหรือโปรแกรมที่พิมพ์nองค์ประกอบแรกของลำดับA137613(ลำดับที่ไม่มี 1s) คุณสามารถอ่านหมายเลขอินพุตn >= 0ผ่านทาง STDIN, อาร์กิวเมนต์บรรทัดคำสั่ง, พรอมต์หรือฟังก์ชันอาร์กิวเมนต์ เอาท์พุทnองค์ประกอบแรกในรูปแบบที่อ่านได้ไปยัง STDOUT หรือส่งกลับอาร์เรย์หรือรายการที่มีค่าเหล่านี้

นี่คือรหัสกอล์ฟ ดังนั้นรหัสที่สั้นที่สุดชนะ

ลีดเดอร์บอร์ด:

นี่คือ Stack Snippet เพื่อสร้างทั้งกระดานผู้นำปกติและภาพรวมของผู้ชนะตามภาษา เพื่อให้แน่ใจว่าคำตอบของคุณปรากฏขึ้นโปรดเริ่มคำตอบด้วยหัวข้อโดยใช้เทมเพลต Markdown ต่อไปนี้:

# Language Name, N bytes

โดยที่ N คือขนาดของการส่งของคุณ หากคุณปรับปรุงคะแนนของคุณคุณสามารถเก็บคะแนนเก่าไว้ในบรรทัดแรกโดยการตีพวกเขาผ่าน ตัวอย่างเช่น

# Ruby, <s>104</s> <s>101</s> 96 bytes

var QUESTION_ID=55272;function answersUrl(e){return"http://api.stackexchange.com/2.2/questions/"+QUESTION_ID+"/answers?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+ANSWER_FILTER}function getAnswers(){jQuery.ajax({url:answersUrl(page++),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){answers.push.apply(answers,e.items),e.has_more?getAnswers():process()}})}function shouldHaveHeading(e){var a=!1,r=e.body_markdown.split("\n");try{a|=/^#/.test(e.body_markdown),a|=["-","="].indexOf(r[1][0])>-1,a&=LANGUAGE_REG.test(e.body_markdown)}catch(n){}return a}function shouldHaveScore(e){var a=!1;try{a|=SIZE_REG.test(e.body_markdown.split("\n")[0])}catch(r){}return a}function getAuthorName(e){return e.owner.display_name}function process(){answers=answers.filter(shouldHaveScore).filter(shouldHaveHeading),answers.sort(function(e,a){var r=+(e.body_markdown.split("\n")[0].match(SIZE_REG)||[1/0])[0],n=+(a.body_markdown.split("\n")[0].match(SIZE_REG)||[1/0])[0];return r-n});var e={},a=1,r=null,n=1;answers.forEach(function(s){var t=s.body_markdown.split("\n")[0],o=jQuery("#answer-template").html(),l=(t.match(NUMBER_REG)[0],(t.match(SIZE_REG)||[0])[0]),c=t.match(LANGUAGE_REG)[1],i=getAuthorName(s);l!=r&&(n=a),r=l,++a,o=o.replace("{{PLACE}}",n+".").replace("{{NAME}}",i).replace("{{LANGUAGE}}",c).replace("{{SIZE}}",l).replace("{{LINK}}",s.share_link),o=jQuery(o),jQuery("#answers").append(o),e[c]=e[c]||{lang:c,user:i,size:l,link:s.share_link}});var s=[];for(var t in e)e.hasOwnProperty(t)&&s.push(e[t]);s.sort(function(e,a){return e.lang>a.lang?1:e.lang<a.lang?-1:0});for(var o=0;o<s.length;++o){var l=jQuery("#language-template").html(),t=s[o];l=l.replace("{{LANGUAGE}}",t.lang).replace("{{NAME}}",t.user).replace("{{SIZE}}",t.size).replace("{{LINK}}",t.link),l=jQuery(l),jQuery("#languages").append(l)}}var ANSWER_FILTER="!t)IWYnsLAZle2tQ3KqrVveCRJfxcRLe",answers=[],page=1;getAnswers();var SIZE_REG=/\d+(?=[^\d&]*(?:&lt;(?:s&gt;[^&]*&lt;\/s&gt;|[^&]+&gt;)[^\d&]*)*$)/,NUMBER_REG=/\d+/,LANGUAGE_REG=/^#*\s*([^,]+)/;

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>

Pyth, 14 13 ไบต์

meaYhP+sY-5dQ

ใช้a(n) = Lpf(6 - n + sum(a(i) for i in range(1, n))ที่Lpfหมายถึงปัจจัยสำคัญอย่างน้อยที่สุด

ลองที่นี่ออนไลน์

Python 3.5.0b1 +, 95 93 ไบต์

ลิงก์ไปยังPython เวอร์ชัน 3.5.0b1 +

import math
def f(k,n=2,a=7,L=[]):x=math.gcd(n,a);return k and f(k-1%x,n+1,a+x,L+1%x*[x])or L

การนำไปปฏิบัติโดยตรงของการเกิดซ้ำซึ่งมี:

• เพื่อนที่ดีของเรา1%xและ
• math.gcdfractions.gcdเมื่อเทียบกับ

Julia, 110 ไบต์

n->(a(n)=(n≥1&&(n==1?7:a(n-1)+gcd(n,a(n-1))));i=2;j=0;while j<n x=a(i)-a(i-1);x>1&&(j+=1;println(x));i+=1end)

Ungolfed:

function a(n::Int)
    n ≥ 1 && (n == 1 ? 7 : a(n-1) + gcd(n, a(n-1)))
end

function f(n::Int)
    i = 2;
    j = 0;
    while j < n
        x = a(i) - a(i-1)
        if x > 1
            j += 1
            println(x)
        end
        i += 1
    end
end

PHP, 101 96 99 98 77 72 ไบต์

<?for(;2>$t=gmp_strval(gmp_gcd(~++$i,7+$e+=$t))or$argv[1]-=print"$t ";);

การใช้งาน:
โทรหาสคริปต์พร้อมอาร์กิวเมนต์: php -d error_reporting=0 script.php 30
หากคุณต้องการทดสอบคุณต้องยกเลิกการคอมเม้นต์;extension=php_gmp.dllใน php.ini ของคุณ
-> extension=php_gmp.dll
ฉันควรเพิ่มส่วนขยายไปยังจำนวนไบต์ของฉันหรือไม่ ความคิดใด ๆ

บันทึก:
บันทึกแล้ว 3 ไบต์ขอบคุณ Ismael Miguel
บันทึก 26 ไบต์ขอบคุณ primo

Haskell, 51 ไบต์

d=zipWith gcd[2..]$scanl(+)7d
f=($filter(>1)d).take

โปรดทราบว่าfเป็นฟังก์ชั่นที่จะคืนองค์ประกอบnแรก

มากกว่าการคำนวณa(n)แล้วการทำงานออกความแตกต่างที่เราคำนวณความแตกต่างและรวมพวกเขาร่วมกันที่จะได้รับd(n) a(n)(ผู้ที่ไม่คุ้นเคยกับ Haskell อาจประท้วงว่าเราต้องการa(n)ก่อนเพื่อให้ได้d(n)แต่แน่นอนว่าการประเมินผลที่ขี้เกียจทำให้เราได้รับปัญหานี้!)

Ungolfed:

a = scanl (+) 7 d        -- yielding a(n) = 7 + d(1) + d(2) + ... + d(n-1)
d = zipWith gcd [2..] a  -- yielding d(n) = gcd(n+1, a(n))

f n = take n $ filter (> 1) d -- get rid of 1s and take the first n

Pyth, 30 ไบต์

กอล์ฟที่แย่มากสามารถลดได้อย่างมาก กำหนดฟังก์ชั่นวนซ้ำที่ด้านหน้าตัวกรอง.first-n จากนั้นจับคู่ความแตกต่าง

L?tb+KytbibK7m-yhdyd.ft-yhZyZQ

ลองมันนี่เกมออนไลน์

Julia, 69 67 ไบต์

n->(i=1;a=7;while n>0 x=gcd(i+=1,a);a+=x;x>1&&(n-=1;println(x))end)

นี่เป็นวิธีแก้ปัญหาแบบวนซ้ำที่ง่ายสำหรับปัญหา xคือความแตกต่าง (ซึ่งเป็นgcd) และการปรับปรุงแล้วฉันโดยการเพิ่มax

JavaScript (ES6), 91

gcd แบบเรียกซ้ำ, ฟังก์ชั่นหลักซ้ำ ไม่เร็วนัก

หมายเหตุปกติ: ทดสอบการเรียกใช้ตัวอย่างข้อมูลบนเบราว์เซอร์ที่สอดคล้องกับ EcmaScript 6 (ไม่ใช่ Chrome ไม่ใช่ MSIE โดยเฉพาะอย่างยิ่งฉันทดสอบบน Firefox, Safari 9 แล้ว)

F=m=>{
  for(G=(a,b)=>b?G(b,a%b):a,o=[],p=7,n=1;m;d>1&&(o.push(d),--m))
    p+=d=G(++n,p);
  return o
}

O.innerHTML=F(+I.value)

<input id=I value=10><button onclick='O.innerHTML=F(+I.value)'>-></button>
<pre id=O></pre>

Haskell, 74 71 66 ไบต์

f=($filter(>1)$tail>>=zipWith(-)$scanl(\x->(x+).gcd x)7[2..]).take

ใช้เคล็ดลับที่นี่: https://codegolf.stackexchange.com/a/39730/43318และทำให้ไม่มีจุด

(ก่อนหน้า: 71 ไบต์)

a=scanl(\x->(x+).gcd x)7[2..]
f m=take m$filter(>1)$zipWith(-)(tail a)a

ก่อนอื่นให้สร้างลำดับของ a แล้วจึงหาความแตกต่าง

(ก่อนหน้า: 74 ไบต์)

f m=take m$filter(>1)$map snd$scanl(\(x,d)->(\y->(x+y,y)).gcd x)(7,1)[2..]

ฟังก์ชั่นรายการมาตรฐานรวมถึงการใช้ฟังก์ชั่นแลมบ์ดาอย่างชาญฉลาด โปรดทราบว่านี่คือ 1 ไบต์สั้นกว่าที่ชัดเจนยิ่งขึ้น

g m=take m$filter(>1)$map snd$scanl(\(x,d)n->(x+gcd x n,gcd x n))(7,1)[2..]

หากเราไม่นับการนำเข้าฉันสามารถทำให้มันลดลงเหลือ 66

import Data.List
h m=take m$filter(>1)$snd$mapAccumL(\x->(\y->(x+y,y)).gcd x)7[2..]

PARI / GP, 60 ไบต์

a(n)=a=7;i=1;while(n,if(1<d=gcd(i+=1,a),n-=1;print(d));a+=d)

ถ่ายภาพโดยตรงจากคำนิยามa (n) - a (n-1) = gcd (n, a (n-1))

เอาท์พุทสำหรับa(20):

5
3
11
3
23
3
47
3
5
3
101
3
7
11
3
13
233
3
467
3

C ++, 193 182 180 172 ไบต์

ขอบคุณ @Jakube - บันทึก 8 ไบต์กับเอาต์พุต

int g(int a,int b){return a==b?a:a>b?g(b,a-b):g(a,b-a);}void f(int *r,int n){int m=1,i=0,a=7,b;while(i<n){b=g(a,++m);if(b>1){r[i]=b;++i;}a+=b;}}int main(){int r[6];f(r,6);}

Mathematica, 59 ไบต์

For[i=j=1;a=7,i<=#,,d=GCD[++j,a];If[d>1,Print@d;i++];a+=d]&