Problem 659. How long is the longest prime diagonal?

Created by Ned Gulley

Solve Later

{"relationships":[{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/document","targetMode":"","relationshipId":"rId1","target":"/matlab/document.xml"},{"relationshipType":"http://schemas.mathworks.com/matlab/code/2013/relationships/output","targetMode":"","relationshipId":"rId2","target":"/matlab/output.xml"}],"parts":[{"partUri":"/matlab/document.xml","relationship":[],"contentType":"application/vnd.mathworks.matlab.code.document+xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n<w:document xmlns:w=\"http://schemas.openxmlformats.org/wordprocessingml/2006/main\"><w:body><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Stanislaw Ulam once observed that if the counting numbers are</w:t></w:r><w:r><w:t> </w:t></w:r><w:hyperlink w:docLocation=\"http://en.wikipedia.org/wiki/Ulam_spiral\"><w:r><w:t>arranged in a spiral</w:t></w:r></w:hyperlink><w:r><w:t>, the prime numbers contained in it form a surprising pattern. They appear to cluster along diagonals of the spiral matrix.</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Given n, return the length of the longest diagonal sequence of primes in spiral(n).</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Example:</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[ Input  n = 7\n Output d = 4]]></w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"text\"/></w:pPr><w:r><w:t>Since isprime(spiral(n)) is</w:t></w:r></w:p><w:p><w:pPr><w:pStyle w:val=\"code\"/></w:pPr><w:r><w:t><![CDATA[     1     0     0     0     1     0     0\n     0     0     0     1     0     0     0\n     1     0     1     0     0     0     0\n     0     1     0     0     1     1     0\n     0     0     1     0     1     0     1\n     0     1     0     0     0     1     0\n     1     0     0     0     0     0     1]]></w:t></w:r></w:p></w:body></w:document>"},{"partUri":"/matlab/output.xml","contentType":"text/xml","content":"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?><embeddedOutputs><metaData><evaluationState>manual</evaluationState><layoutState>code</layoutState><outputStatus>ready</outputStatus></metaData><outputArray type=\"array\"/><regionArray type=\"array\"/></embeddedOutputs>"}]}

Solve

Solution Stats

1215 Solutions
361 Solvers

Last Solution submitted on Jan 12, 2022

Last 200 Solutions

Problem Comments

9 Comments

9 Comments

Show 6 older comments

Jacob on 11 Jun 2012

I agree that the solution is mis-coded. isprime(spiral(4)) gives:
1 0 0 0
0 0 1 1
1 0 1 0
0 0 0 1
which means that d=2 for this case.

Ned Gulley on 11 Jun 2012

Rescoring based on corrected test suite. Sorry about that.

Jonathan Campelli on 11 May 2015

"find the longest diagonal arrangement of primes in spiral(n)" leaves room for misinterpretation.

I would have quickly understood:
"return the length of the longest diagonal sequence of primes in spiral(n)"

David Verrelli on 17 Nov 2017

Like Jonathan Campelli, I still feel the problem statement could be clearer. From the example I eventually figure that reference is being made to the diagonally oriented elements running between elements (1,5) and (4,2) of isprime(spiral(n)). But arguably these are not on a/the "diagonal" of that matrix, per http://mathworld.wolfram.com/Diagonal.html (cf. https://en.wikipedia.org/wiki/Main_diagonal ).

David Verrelli on 17 Nov 2017

Furthermore, there is no mention in the current online MATLAB documentation that "spiral" is actually a function that is already provided with MATLAB — i.e. that the user does not have to write their own function to make the spiral matrix. Apparently, "spiral can be found in the demo folder" [ref: https://stackoverflow.com/questions/29466058/spiral-loop-on-a-matrix-from-a-point].

Tung Nguyen on 29 Oct 2019

hard

Highphi on 24 Jul 2020

eh, I think the problem statement made sense. especially if you just run 'isprime(spiral(n))' in your command window. I liked this problem. Ended up using the technique I almost used on another problem awhile back but it made more sense to use here, so I had the basic framework already started :D

mahesh kumar on 8 Nov 2020

for i=1:50
spy(isprime(spiral(i)))
pause(0.1)
end

Dara Murray on 19 Jul 2021

Nice Problem

Solution Comments

Solution 6037528

Solution 5954230

1 Comment

1 Comment

J.R.! Menzinger on 6 Aug 2021

Dear Bora,
You asked for help resolving this but I was "far away from Cody" for a long time... ;-)

I celebrate you were at least successful with it.

Solution 2045288

1 Comment

1 Comment

Miles Walker on 4 Dec 2019

Interesting! This error stems from using "Ix" in both function and sub-function. As I test my code in scripts not contained within a function I was missing it.

Why are parameters made available to sub-functions by default?

Solution 1854983

1 Comment

1 Comment

jianrui fan on 21 Jun 2019

this is really awesome and impressive.I do learn something

Solution 1699697

Solution 593819

1 Comment

1 Comment

jianrui fan on 21 Jun 2019

thank you very much

Solution 443569

1 Comment

1 Comment

Ramoflaple on 11 Aug 2015

a good precise algorithm

Solution 117430

2 Comments

2 Comments

Jon on 21 Mar 2014

Best solution that doesn't cheat with a lookup or regexp

empad on 7 Dec 2020

Think I might have found a problem with merging the original and rotated matrices. In the case of n=7, if you change (1,6) and (2,7) to ones, the output becomes 6 when it should still be 4. If you evaluate the matrices separately then it should work fine, but apart from that this is a really nice solution

Solution 102911

1 Comment

1 Comment

Dirk Engel on 26 Jun 2012

max(diff(find(~zeroOneVector)))-1 returns the longest run of ones in a vector of zeros and ones, but it FAILS IN SEVERAL CASES, e.g.:
- all ones: [1 1 1 1 1] returns empty
- just a single zero: [1 1 0 1 1] returns empty
- all zeros at one end: [1 1 1 0 0] returns 0
- longest run at one end: [1 1 0 1 0] returns the longest run between zeros, here 1

As a workaround for all cases, add zeros to both ends of the vector:
max(diff(find(~[0 zeroOneVector 0])))-1

However, this is not necessary for this solution because the matrix is very sparse.

Problem Recent Solvers361

Problem Tags

asee diagonal pattern poor_problem_statement primes spiral

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Problem 659. How long is the longest prime diagonal?

Solution Stats

Last 200 Solutions

Problem Comments

Solution Comments

Problem Recent Solvers361

Suggested Problems

More from this Author50

Problem Tags

Community Treasure Hunt

Problem 659. How long is the longest prime diagonal?

Solution Stats

Last 200 Solutions

Problem Comments

Solution Comments

Problem Recent Solvers361

Suggested Problems

More from this Author50

Problem Tags

Community Treasure Hunt

Players