Is this a consecutive-prime/constant-exponent number?

A while ago, I had a look at the prime factorization of 27000:

27000 = 23 × 33 × 53

There are two special things about that:

• consecutive-prime: The primes are consecutive: 2 is the 1st prime, 3 is the 2nd prime, 5 is the 3rd prime.
• constant-exponent: The exponent is the same for every prime (always 3)

Mathematically expressed:

An integer x is a consecutive-prime/constant-exponent number if there exist strictly positive integers n, k, m such that x = pnm × pn+1m × ... × pn+km, where pj is the j-th prime

Your task is to test if a positive integer fulfills these conditions.

Input:

A positive integer > 1, in any reasonable form.

Output:

One of two values, at least one of which has to be constant, indicating whether the input is a consecutive-prime/constant-exponent number.

Edge cases:

• primes are truthy, as the factorization for prime p is p1
• other numbers that can be written as pm where p is a prime are also truthy.

Rules:

• Standard loopholes apply.
• No worries about integer overflow, but numbers up to 255 must work.
• Shortest code in bytes wins.

Test cases:

Truthy:

2
3
4
5
6
7
8
9
11
13
15
27000
456533


Falsy:

10
12
14
72
10000000


Here is a python script generating some test cases.

The fact that I accepted an answer does not mean that the challenge is over; the winner can still change!

• You could probably come at this the other way by generating a list of all such numbers and checking if the input is in the list May 16, 2018 at 15:48
• @EngineerToast There are infinitely many truthy numbers though. May 16, 2018 at 21:58
• @AlexisOlson Sure, but a finite that can be handled as integers by many languages. May 16, 2018 at 22:11
• Your mathematical expression has Pj doesn't relate to the x = Pn^m part. I'm assuming you meant Pn is the n-th prime May 16, 2018 at 22:37
• @Veskah n has a specific value (index of the first prime dividing x), so saying Pn is the n-th prime is awkward if you also want to imply that Pn+1 is the n+1-th prime. May 17, 2018 at 3:53

05AB1E, 4 bytes

Ó0ÛË


Try it online!

Explanation

Ó     # get a list of prime exponents
Ë  # all remaining elements are equal

• ÎÓKË is all I can think of other than this, nice one... I was thinking ß but it does the opposite of what I thought. May 18, 2018 at 1:53

Regex (ECMAScript), 276205201193189177 174 bytes

Comparing the multiplicities (exponents) of different prime factors is an interesting problem for solving with ECMAScript regex – the lack of backreferences that persist through iterations of a loop makes it a challenge to count anything. Even if counting the numerical trait in question is possible, often a more indirect approach makes for better golf.

As with my other ECMA regex posts, I'll give a spoiler warning: I highly recommend learning how to solve unary mathematical problems in ECMAScript regex. It's been a fascinating journey for me, and I don't want to spoil it for anybody who might potentially want to try it themselves, especially those with an interest in number theory. See this earlier post for a list of consecutively spoiler-tagged recommended problems to solve one by one.

So do not read any further if you don't want some advanced unary regex magic spoiled for you. If you do want to take a shot at figuring out this magic yourself, I highly recommend starting by solving some problems in ECMAScript regex as outlined in that post linked above.

The main payload from a regex I previously developed turned out to be very applicable to this challenge. That is the regex that finds the prime(s) of highest multiplicity. My first solution for that was very long, and I later golfed it way down in stages, first rewriting it to use molecular lookahead, and then porting it back to plain ECMAScript using an advanced technique to work around the lack of molecular lookahead, and subsequently golfing it down to be much smaller than the original plain ECMAScript solution.

The part from that regex that applies to this problem is the first step, which finds Q, the smallest factor of N that shares all of its prime factors. Once we have this number, all we have to do to show that N is a "constant-exponent number" is divide N by Q until we can't any longer; if the result is 1, all primes are of equal multiplicity.

After submitting an answer using my previously developed algorithm for finding Q, I realized that it could be calculated in an entirely different way: Find the largest square-free factor of N (using the same algorithm as my Carmichael number regex). As it turns out, this poses no difficulty at all in terms of stepping around the lack of molecular lookahead and variable-length lookbehind (no need to pull in the advanced technique previously used), and is 64 bytes shorter! [Note, at the time, molecular lookahead still offered a slight benefit, but thanks to later golfs, it no longer does at all.] Additionally it eliminates the complexity of treating square-free N and prime N as different special cases, dropping another 7 bytes from this solution.

(There still remain other problems that require the advanced technique formerly used here to golf down the calculation of Q, but currently none of them are represented by my PPCG posts.)

I put the multiplicity test before the consecutive-primes test because the latter is much slower; putting tests that can fail more quickly first makes the regex faster for uniformly distributed input. It's also better golf to put it first, because it uses more backreferences (which would cost more if they were double-digit). [The multiplicity test is now split in two for golf, with the consecutive-primes test sandwiched inside.]

I was able to drop 4 bytes from this regex (193 → 189) using a trick found by Grimmy that can futher shorten the form of division in which the assertion that $$\dividend-quotient\$$ is divisible by $$\divisor-1\$$ can be skipped. (Skipping this assertion was already saving 4 bytes, so now it's saving 8 bytes.) This algorithm (in both its 4 byte longer form and the currently used one) is guaranteed to match the correct quotient when at least one of these constraints is met (along with satisfying the assertion that $$\dividend\$$ is divisible by $$\divisor\$$):

• $$\divisor^2+2(divisor)+1 < 4(dividend)\$$
• $$\quotient=1\$$  or equivalently  $$\divisor=dividend\$$
• $$\divisor\$$ is a prime power
• One of the cases in which this must be true is when $$\dividend\$$ is semi-prime and neither of the two other conditions are met. It also must be true if $$\dividend\$$ is a prime power.

We're actually not guaranteed in this context that any of those constraints will be met. Most notably, we're not even asserting that $$\quotient\$$ (which will become the next $$\dividend\$$) is still divisible by $$\divisor\$$ after each iteration of the loop. But an additional fact comes to the rescue: if the division algorithm matches an incorrect quotient (which it doesn't always do even when none of the above constraints are met), the incorrect answer is guaranteed to be less than $$\divisor\$$ but greater than $$\1\$$.

So if this happens when repeatedly dividing N by Q, the loop will terminate, and the final result will not be 1, thus guaranteeing no false positives. If N is actually a constant-exponent number, the correct quotient will always be greater than or equal to Q, so this doesn't cause any false negatives either.

Further golfs:

• -3 bytes by dividing by N just like any other largest square-free factor, instead of treating N being square-free as a special case (inspired by H.PWiz but goes farther then what he did, saving 1 extra byte – he made it so the regex attempted to divide by 0 in the case that Q=N, and then asserted that the end result of the repeated division was Q rather than 1)
• -6 bytes by changing (?=A)(?!B) to (?!B)A (some of this is from multiple uses of backrefs changing from double to single digit, with less of the opposite occurring) – thanks to H.PWiz
• -2 bytes by capturing the smallest prime as $$\\ge 2\$$ instead of augmenting the non-composite test to be a full prime test – thanks to H.PWiz
• -2 bytes by removing an extraneous pair of parentheses – thanks to H.PWiz
• -2 bytes by changing an outermost (|(x+)\2*(?=\2$)) to (|x*)\1*(?=\1\b), similarly to what I did in the alternate answer here; one of these bytes is from a backref changing from double to single digit And the regex: ^(?=(|x*)\1*(?=\1\b)(((?=(xx+?)\4*$)(?=(x+)(\5+$))\6(?!\4*$))*x$))(?!(((xx+)(?=\9+$)(x+))(?!\8+$)(x+))\7*(?=\7$)(?!(\11\10?)?(xx+)\13+$))((?=((x*)(?=\2\16*$)x)(\15*$))\17)*x$

Try it online!

^                           # tail = N = input number

# Assert that all of N's prime factors are of equal multiplicity. There are two stages
# to this, which we will call Step 1 and Step 2. For golf reasons, a negative assertion
# is inserted between them.

# Step 1: Find Q, the largest square-free factor of N (which will also be the smallest
# factor of N that has all the same prime factors as N).
(?=
(|x*)\1*(?=\1\b)        # cycle tail through all of the divisors of N, including N,
# from largest to smallest;
# \1 = the divisor, or zero if the divisor is N itself
# Assert that tail is square-free (its prime factors all have single multiplicity)
(                       # \2 = tail = the divisor
(
(?=(xx+?)\4*$) # \4 = smallest prime factor of tail (?=(x+)(\5+$))  # \5 = tail / \4 (implicitly); \6 = tool to make tail = \5
\6              # tail = \5
(?!\4*$) # Assert that tail is no longer divisible by \4, i.e. that # that prime factor was of exactly single multiplicity. )*x$
)
)

# Assert that there exists no trio of prime numbers such that N is divisible by the
# smallest and largest prime but not the middle prime.
(?!
(                       # \7 = a factor of N
(                   # \8 = a non-factor of N between \7 and \9
(xx+)(?=\9+$) # \9 = a factor of N smaller than \7 and greater than 1 (x+) # \10 = tool (with \11) to make tail = \9 ) (?!\8+$)
(x+)                # \11 = tool to make tail = \8
)
\7*(?=\7$) # tail = \7 # Assert that \7, \8, and \9 are all prime (?! (\11\10?)? # tail = either \7, \8, or \9 (xx+)\13+$          # Assert tail is not composite; no need to test against
# 0 and 1, because we captured \9 > 1
)
)

# Step 2: Require that the result of repeatedly dividing N by \2 is 1.
(
# In this division calculation \15 = N / \2, we can skip the test for N-\15 being
# divisible by \2-1. (Actually doing that test would make the regex 8 bytes larger.)
# The quotient matched by this abbreviated algorithm is guaranteed to be correct if
# the real quotient is 1 or is not less than the divisor, or if the divisor is a
# prime power. Our input here is not guaranteed to satisfy any of those constraints.
# But if none of those constraints are met, and the quotient matched by the algorithm
# is incorrect, the incorrect quotient is guaranteed to be less than the divisor and
# greater than 1. And once that happens, this loop will terminate, as the current
# value will no longer be divisible by our divisor \2.
(?=
(                   # \15 = tail / \2
(x*)            # \16 = \15-1
(?=\2\16*$) x ) (\15*$)             # \17 = tool to make tail = \15
)
\17                     # tail = \15
)*
x$# Require that the end result is 1  Regex (ECMAScript+(?*), 195187183 bytes Obsoleted by no longer treating N being square-free as a special case (i.e. is now identical to the above). Regex (ECMAScript 2018), 198195194186182 bytes Obsoleted by no longer treating N being square-free as a special case and changing (?=A)(?!B) to (?!B)A (i.e. is now identical to the above). • You can replace .*$\2 with \2^ Apr 22, 2019 at 5:36
• Although, I believe this is valid: ^(?=(|(x+)\2*(?=\2$))(((?=(xx+?)\5*$)(?=(x+)(\6+$))\7(?!\5*$))*x$))(?!(((xx+)(?=\10+$)(x+))(?!\9+$)(x+))\8*(?=\8$)(?!(\12\11?)?(xx+)\14+$))((?=((x*)(?=\2\17+$)x)(\16*$))\19)*\3$ Apr 22, 2019 at 6:28
• Doesn't seem close to optimal though Apr 22, 2019 at 6:29
• @H.PWiz thanks! Note, the \19 at the end is supposed to be a \18. Mar 29, 2021 at 19:00

Jelly, 136 5 bytes

ÆEt0E


Try it online!

Still outgolfed... (thanks Erik for -1 byte)

Explanation

ÆE     # get a list of prime exponents (noooo long builtin name)
t0   # remove zeroes on both sides (leading or trailing)
E  # all remaining elements are equal

• œl -> t. There's no reason for trailing zeroes to be present in ÆE's output. May 16, 2018 at 16:01
• May 16, 2018 at 17:44
• @dylnan That fails for 2250. May 16, 2018 at 17:49
• @Dennis Thanks, I had realized it wouldn't work but I am hoping it will inspire a four byte solution May 16, 2018 at 18:04

JavaScript (ES6), 87 bytes

Returns 0 for truthy or a non-zero integer for falsy.

f=(n,k=2,j,i)=>n%k?j*(P=d=>k%--d?P(d):d==!i)(k)|j-i|(n>1&&f(n,k+1,j||i)):f(n/k,k,j,-~i)


Try it online!

Commented

f = (                     // f() = recursive function taking:
n,                      //   n = input
k = 2,                  //   k = current factor
j,                      //   j = reference exponent, initially undefined
i                       //   i = current exponent, undefined each time we start testing
) =>                      //       the next factor
n % k ?                 // if k is not a divisor of n:
j * (                 //   ignore the primality of k if j is still undefined
P = d =>            //     P() = function testing if k is prime:
k % --d ?         //       decrement d; if d is not a divisor of k:
P(d)            //         do a recursive call until it is
:                 //       else:
d == !i         //         unless i is already defined: d must not be equal to 1
//         (if it is: k is the next prime but does not divide n)
)(k) |                //   initial call to P() with d = k
j - i | (             //   if both i and j are defined, they must be equal
n > 1 &&            //   if n is not yet equal to 1,
f(n, k + 1, j || i) //   go on with k + 1; if j is undefined, set it to i
)                     //   (otherwise, stop recursion and return what we have)
:                       // else:
f(n / k, k, j, -~i)   //   increment the current exponent and go on with n / k

• This was broken by the change of j||i to i. It now yields lots of false positives. Feb 11, 2019 at 18:03
• @Deadcode I can't check or correct this for the moment, so I've just rollbacked for now. Feb 11, 2019 at 19:14

CJam, 30 29 bytes

{mFz~)-!\__W=,\0=>\-:mp1#W=&}


Try it online!

My first answer after a nearly 2(!)-year break, so it can probably be golfed more. This is a block that takes input as an integer (can also be mapped over for arrays of integers).

Explanation

{        e# Begin block
mF      e# Factor input, giving an array of primes and their powers
z~      e# Transpose and dump, giving an array of primes and an array of powers
)-      e# Check that the powers are the same: subtract each power from the last element
!       e# Negate to make sure they're all 0
\__W=,  e# Get the range from 0 to the largest prime minus one
\0=>    e# Slice that array so it only includes everything larger than the smallest prime
\-      e# Remove the original primes from the range array
:mp     e# Check each element for primality. If the input's primes are consecutive,
e# this will contain no primes
&       e# If the powers are the same AND primes are consecutive, return 1. Otherwise, 0.
}


Stax, 5 6 bytes

╣♥qJ╬c


Run and debug it

Unpacked, ungolfed, and commented, it looks like this.

|n    get the exponents of the prime factorization
:u    array has exactly a single distinct element


Edit: This doesn't work on 512. I'll give it some thought and hopefully a fix later. It works now.

Stax, 9 bytes

1 is truthy, 0 is falsy

αAG<└\{┬⌠


Run and debug it

Explanation

|nX0-u%x:^=      # Full Program, unpacked, implicit input
|n               # Exponents of sequential primes in factorization. (eg. 20 -> [2 0 1])
X              # Save to X register
0-            # Remove all '0' from array
u%          # Get unique numbers and get length of array
x         # Copy back the array saved to X
:^       # Is it ascending
=       # Are the two comparisons equal? implicit output


Can probably be golfed more, but it covers the cases I was missing in the last solution.

Husk, 11 bytes

Λ§*≈¤≈oṗ←gp


Try it online!

Outputs 0 if not a consecutive-prime/constant-exponent number.

Java 10, 223191178176 168 bytes

n->{var s=new java.util.HashSet();for(int f=1,i=1,x,j;n>1;){for(x=++i,j=2;j<x;)x=x%j++<1?1:x;if(x>1){for(j=0;n%i<1&&n>(f=0);n/=i)j++;if(f<1)s.add(j);}}return s.size();}


Returns 1 as truthy, and >=2 as falsey.

Try it online.

Explanation:

n->{                   // Method with integer parameter and boolean return-type
var s=new java.util.HashSet();
//  Set to keep track of the prime exponents
for(int f=1,         //  Prime-flag, starting at 1
i=1,x,j;     //  Index and temp integers
n>1;){       //  Loop as long as n is still larger than 1
for(x=++i,         //   Set x to i, after we've increased i by 1 first with ++i
j=2;           //   Set j to 2 (first prime)
j<x;)          //   Inner loop as long as j is still smaller than x
x=x%j++<1?       //    If x is divisible by j:
1             //     Set x to 1
:              //    Else:
x;            //     Leave x unchanged
if(x>1){           //    If x is larger than 1 (if i is a prime):
for(j=0;         //     Use j as counter, and set it to 0
n%i<1        //     If n is divisible by i:
//      And loop as long as n is still divisible by i,
&&n>         //      and n is larger than 0
(f=0);   //      (and set f to 0 at the same time)
n/=i)        //       Divide n by i
j++;           //       And increase j by 1
if(f<1)          //     If the flag f is now/still 0:
s.add(j);}}    //      Add counter j to the Set
return s.size();}    //  Return the amount of items in the Set
//  (1 being true, >=2 being false)


Some example inputs:

n=15:

• Flag remains 1 for the first prime 2 (because 15 is not divisible by 2).
• Flag goes from 1 to 0 as soon as we're at the prime 3. Since 15 is divisible by 3, n becomes 5 (15/31), and the Set becomes [] → [1].
• Then we check the next prime 5. Since 5 is divisible by 5, n becomes 1 (5/51), and the Set remains the same ([1] → [1]).
• Now n=1, so we stop the outer loop. The Set ([1]) only contains one item, the 1 from both adjacent primes 3 and 5, so we return true.

n=14:

• Flag goes from 1 to 0 for the first prime 2 (because 14 is divisible by 2). n becomes 7 (14/21), and the Set becomes [] → [1].
• Then we check the next prime 3. Since 7 is not divisible by 3, n remains the same, and the Set becomes [1] → [1,0].
• Then we check the next prime 5. Since 7 is also not divisible by 5, n remains the same, and the Set remains the same as well ([1,0] → [1,0]).
• Then we check the next prime 7. Since 7 is divisible by 7, n becomes 1 (7/71), and the Set remains the same ([1,0] → [1,0]).
• Now n=1, so we stop the outer loop. The Set ([1,0]) contains two items, the 1 from the non-adjacent primes 2 and 7, and the 0 from the primes 3 and 5, so we return false.

n=72:

• Flag goes from 1 to 0 for the first prime 2, because 72 is divisible by 2 (multiple times). So n becomes 9 (72/23), and the Set becomes [] → [3].
• Then we check the next prime 3. Since 9 is divisible by 3 (multiple times), n becomes 1 (9/32), and the Set becomes [3] → [3,2].
• Now n=1, so we stop the outer loop. The Set ([3,2]) contains two items, the 3 from prime 2, and the 2 from prime 3, so we return false.
• You can remove <2 and return an int (specify that you return 1 for truthy). May 20, 2018 at 18:37
• @wastl Ah, missed the rule about only one of the two values being consistent. In that case 1 is truthy and 2 or higher is falsey. Thanks. May 20, 2018 at 20:27
• Thanks to whoever gave me the bounty, but why? Feb 12, 2019 at 7:05
• I had started a "Reward existing answer" bounty to bring more attention to my ECMAScript answer, which I still believe deserves more than it's gotten (I'd consider the bounty to have failed). When the week was over, I had to choose an answer other than my own to award the bounty to, or leave it to default to the highest-voted. I did not think any deserved it, but your answer had the best explanation, and that's why I awarded it to you; good explanations are way too rare on PPCG. As for my answer, I figure I have to give it a better writeup, which I plan on when I have time. Feb 15, 2019 at 20:35
• @Deadcode Ah, so that's why. I thought maybe someone started the bounty, but they accidentally let it expire and it came to me instead. Still a bit confused why my answer then and not the highest voted one. I must say I'm impressed by all your Regex answers. I've seen some of them, and every time I'm amazed. Especially when I come back later to the same answer and you golfed it even loads more. :D I just noticed I hadn't seen nor upvoted the one for this challenge, so I just did. You know, I will add a bounty to this answer of yours. :) Feb 15, 2019 at 21:16

MATL, 12 11 10 bytes

YFtgYsg)Zs


Try it at MATL Online!

Thanks to Luis Mendo for the remove-leading-zeroes part. He also pointed out that swapping truth values is allowed, so this returns 0 for numbers that satisfy the challenge requirements and any positive value otherwise.

Grosso Modo, this generates the exponents of the sequential prime factorization, removes leading zeroes and calculates the standard deviation.

• I think 0iYFhdz works for 7 bytes: prepend a 0 to the exponents of the sequential factorization, consecutive differences, number of non-zeros. The result is 1 iff input satisfies the requirement May 18, 2018 at 10:48
• @LuisMendo Sorry for the delayed response, but you can post that as a separate answer. It is most definitely very different. May 20, 2018 at 7:04
• OK, I posted it as an answer May 20, 2018 at 12:39

J, 16 bytes

Big thanks to FrownyFrog for -8 bytes!

(=&#+/\=@#])_&q:


Try it online!

J, 24 bytes

[:(1=[:#@~.{.@I.}.])_&q:


Try it online!

Explanation:

_&q: prime exponents

{.@I.}.] removes the leading zeros by finding the first non-zero element:

     }.   drop
]  from the list of exponents
{.@       as much items as the first of the
I.     indices of non-zero elements


1=[:#@~. tests if all remaining numbers are equal:

  [:#@~.  finds the length of the list after removing the duplicates
1=        is it 1?


MATL, 7 bytes

0iYFhdz


The result is 1 iff input satisfies the requirement.

Explanation

0    % Push 0
i    % Push input number
YF   % Exponents of consecutive prime factors
h    % Concatenate horizontally
d    % Consecutive differences
z    % Number of nonzeros. Implicitly display


Octave, 67 bytes

@(x)~any(diff(find(h=histc(factor(x),primes(x))))-1)&h(h>0)==max(h)


Try it online!

I believe this is the only solution that uses a histogram.

Explanation:

This makes a histogram, where the variable to be counted are the factors of the input, and placed in the bins primes(x), which is all primes less than the input. We then find the location of the prime factors, takes the difference between each of the indices and subtract one. If there are any elements that aren't zero (i.e. the difference of the indices of prime numbers is not 1), then this will result in a falsy value, otherwise it will return a truthy value.

We then chech that all non-zero elements in the histogram are equal to the maximum element. If there are values that aren't equal then this will result in a falsy value, otherwise it will return a truthy value.

If both those blocks are truthy then our input is a consecutive prime constant exponent number!

APL (Dyalog Extended), 28 bytes

{f p←↓⍭⍵⋄(1=≢∪p)∧∨/f⍷⍸1⍭⍳⍵}


Try it online!

How:

{f p←↓⍭⍵⋄(1=≢∪p)∧∨/f⍷⍸1⍭⍳⍵} ⍝ Monadic function, takes an argument ⍵
⍭⍵                     ⍝ Prime factors and exponents of ⍵
↓                        ⍝ split the resulting matrix in 2 vectors
f p←                          ⍝ assign the factors to f and the powers to p
⋄                     ⍝ then
⍳⍵   ⍝ range [1..⍵]
1⍭     ⍝ primality check for each element in the vector
⍸        ⍝ where; returns the indices of truthy values
f⍷         ⍝ find the factors; returns a boolean vector
∨/           ⍝ logical OR reduction
∧             ⍝ logical AND
(   ∪p)              ⍝ unique members of the powers
≢                 ⍝ tally; returns the number of elements in the vector
1=                  ⍝ check if there's only one element


Wolfram Language (Mathematica), 65 bytes

Differences[{PrimePi@#,#2}&@@@FactorInteger@#]~MatchQ~{{1,0}...}&


Try it online!

Pari/GP, 63 bytes

n->#Set(vector(#(a=factor(n)~),i,[primepi(a[1,i])-i,a[2,i]]))<2


Try it online!

J, 14 bytes

1#.2~:/\0,_&q:


1 in the output indicates consecutive constant exponent.

Try it online!

        0,_&q:   zero followed by the prime exponents of input
2~:/\         for every two consecutive values, 1 if they are different
1#.              convert from base-1, just add them up


Clean, 127 bytes

import StdEnv
@n=[]== $n ?n#j=$n
= @n||j==filter@[hd j..last j]&&any(\p=(prod j)^p==n)[1..n]
$n=[i\\i<-[2..n-1]|n/i*i==n&& @i]  Try it online! Defines the function ? :: Int -> Bool using $ :: Int -> [Int] to factorize and @ :: Int -> Bool to check primality.

APL(NARS) 41 chars, 82 bytes

{(1=≢∪+/¨{v=⍵⊃v}¨⍳≢v)∧(1↓w)≡¯1↓1πw←∪v←π⍵}


{π⍵} is the function factorization of argument ⍵ in the list of prime factors (repeat if one prime appear more time);
{1π⍵} is the function next prime (note that in this case its argument is not a scalar but one array of integers). test:

  h←{(1=≢∪+/¨{v=⍵⊃v}¨⍳≢v)∧(1↓w)≡¯1↓1πw←∪v←π⍵}
(2..30)/⍨h¨2..30
2 3 4 5 6 7 8 9 11 13 15 16 17 19 23 25 27 29 30
h¨27000 456533 72 10000000
1 1 0 0
`