November 15, 2019 10:00 AM
We have a fun little task for your Friday lunch hour:
Given an array of integers, rearrange the integers so no two consecutive integers have a sum that is divisible by 3, or indicate that such an arrangement is not possible.
Your task is to write a program to rearrange an array of integers as described above. When you are finished, you are welcome to read or run a suggested solution, or to post your own solution or discuss the exercise in the comments below.
Posted by programmingpraxis
Categories: Exercises
Tags:
Mobile Site | Full Site
Get a free blog at WordPress.com Theme: WordPress Mobile Edition by Alex King.
As you say, it comes down to taking the integers mod 3, then ensuring that 0 is never next to 0 and 1 is never next to 2. It’s a nice exercise to find all the ways this can be done, given the number of integers with each modulus. Here’s some Haskell, it would benefit from memoization or dynamic programming:
-- sequences with i 0's, j 1's, k 2's with: -- 0 and 0 not adjacent -- 1 and 2 not adjacent -- f0 i j k - sequences beginning 0 -- f1 i j k - sequences beginning 1 -- f2 i j k - sequences beginning 2 -- f i j k - all sequences -- 0 is followed either by nothing or by 1 or by 2 f0 1 0 0 = [[0]] f0 i j k | j + k + 1 < i = [] f0 i j k | i > 0 = map (0:) (f1 (i-1) j k ++ f2 (i-1) j k) f0 i j k = [] -- 1 is followed by nothing or by 0 or by 1 f1 0 1 0 = [[1]] f1 i j k | j + k < i = [] f1 i j k | j > 0 = map (1:) (f0 i (j-1) k ++ f1 i (j-1) k) f1 i j k = [] -- 2 is followed by nothing or by 0 or by 2 f2 0 0 1 = [[2]] f2 i j k | j + k < i = [] f2 i j k | k > 0 = map (2:) (f0 i j (k-1) ++ f2 i j (k-1)) f2 i j k = [] f i j k = f0 i j k ++ f1 i j k ++ f2 i j k main = mapM_ putStrLn (map (map ("012"!!)) (f 4 2 2))By matthew on November 15, 2019 at 7:38 PM
Incidentally, for your example with 9 0’s and 2’s and just 2 1’s, there are 4322 different possible orderings, from 01010202020202020222 to 22202020202020201010 (your ordering and the mirror image of the first).
By matthew on November 15, 2019 at 7:50 PM
Here’s a solution in Python.
from collections import defaultdict, deque def tripartition(l): d = defaultdict(list) for x in l: d[x % 3].append(x) out = deque() if d[0]: out.append(d[0].pop()) while d[1]: out.appendleft(d[1].pop()) if d[0]: out.appendleft(d[0].pop()) while d[2] and (not out or out[-1] != 1): out.append(d[2].pop()) if d[0]: out.append(d[0].pop()) return list(out) if len(out) == len(l) else None tests = [ [0, 0, 0], [0, 0, 1], [0, 0, 2], [0, 1, 1], [0, 2, 2], [1, 1, 1], [2, 2, 2], [1, 1, 2], [0, 1, 2], [2, 67, 65, 71, 86, 83, 12, 17, 9, 30, 42, 61, 33, 68, 57, 30, 74, 42, 41, 51], ] for test in tests: print(f'in: {test}') print(f'out: {tripartition(test)}') print()Output:
By Daniel on November 16, 2019 at 12:17 AM
Here’s a similar solution to the one I posted earlier, but without the dependence on the collections module.
def tripartition(l): d = {0: [], 1: [], 2: []} for x in l: d[x % 3].append(x) out = [] if d[0]: out.append(d[0].pop()) while d[1]: out.append(d[1].pop()) if d[0]: out.append(d[0].pop()) out.reverse() while d[2] and (not out or out[-1] != 1): out.append(d[2].pop()) if d[0]: out.append(d[0].pop()) return out if len(out) == len(l) else NoneBy Daniel on November 16, 2019 at 3:02 AM
Here’s a dynamic programming solution in Python that just computes the number of possible sequences. It’s not to difficult to adapt to compute the sequences themselves. solveseq(s) gives the number of rearrangements (that preserve the ordering for each modulus class) for s. Use numpy for convenient multidimensional arrays – we can also specify the array type, type ‘object’ gives us bignums, or we can use a floats to avoid overflow in bigger problems:
import numpy def solve(I,J,K,atype=int): a = numpy.zeros([3,I+1,J+1,K+1],atype) a[0,1,0,0] = a[1,0,1,0] = a[2,0,0,1] = 1 for n in range(2,I*J*K+1): for i in range(0,min(I+1,n+1)): for j in range(max(0,n-i-K),min(J+1,n-i+1)): k = n-i-j # k = n-i-j <= K so n-i-K <= j if i > 0: a[0,i,j,k] = a[1,i-1,j,k] + a[2,i-1,j,k] if j > 0: a[1,i,j,k] = a[0,i,j-1,k] + a[1,i,j-1,k] if k > 0: a[2,i,j,k] = a[0,i,j,k-1] + a[2,i,j,k-1] return a[0,I,J,K]+a[1,I,J,K]+a[2,I,J,K] def solveseq(s): a = [0,0,0] for i in s: a[i%3] += 1 return solve(*a) s = [2,67,65,71,86,83,12,17,9,30,42,61,33,68,57,30,74,42,41,51] print(solveseq(s)) # 4432 print(solve(9,2,9)) # also 4432 print(solve(10,10,10)) # 41649556 print(solve(50,50,50,object)) # 11648269955901473622920989454259616785721696 print(solve(50,50,50,numpy.float64)) # 1.1648269955901474e+43By matthew on November 16, 2019 at 10:40 AM
We can also use the dynamic programming array to generate uniformly random sample sequences (also fixing a bug with the range for n on line 6 – should be I+J+K+1, and tightening up the bounds on i). Also make the output fruitier:
import numpy import random import time def maketable(I,J,K,atype=int): a = numpy.zeros([3,I+1,J+1,K+1],atype) a[0,1,0,0] = a[1,0,1,0] = a[2,0,0,1] = 1 for n in range(2,I+J+K+1): # n = i+j+k <= i+J+K so n-J-K <= i for i in range(max(0,n-J-K),min(I,n)+1): # k = n-i-j <= K so n-i-K <= j for j in range(max(0,n-i-K),min(J,n-i)+1): k = n-i-j if i > 0: a[0,i,j,k] = a[1,i-1,j,k] + a[2,i-1,j,k] if j > 0: a[1,i,j,k] = a[0,i,j-1,k] + a[1,i,j-1,k] if k > 0: a[2,i,j,k] = a[0,i,j,k-1] + a[2,i,j,k-1] return a def select(*s): x,t,p = random.random(),sum(s),0 for i in range(len(s)-1): p += s[i] if x < p/t: return i return len(s)-1 def randseq(I,J,K): a = maketable(I,J,K,numpy.float64) while True: i,j,k = I,J,K n = [0,1,2][select(a[0,i,j,k],a[1,i,j,k],a[2,i,j,k])] s = [n] for _ in range(i+j+k-1): if n == 0: i -= 1; n = [1,2][select(a[1,i,j,k],a[2,i,j,k])] elif n == 1: j -= 1; n = [0,1][select(a[0,i,j,k],a[1,i,j,k])] elif n == 2: k -= 1; n = [0,2][select(a[0,i,j,k],a[2,i,j,k])] s += [n] yield s s = randseq(9,2,9) for _ in range(10): for n in next(s): print("🍎🍋🍊"[n],end="",flush=True); time.sleep(0.05) print(); time.sleep(0.5)Typical output:
🍎🍊🍎🍊🍎🍊🍊🍊🍎🍊🍎🍋🍎🍊🍊🍎🍋🍎🍊🍎
🍎🍊🍎🍋🍎🍊🍎🍋🍎🍊🍎🍊🍊🍊🍎🍊🍊🍎🍊🍎
🍊🍎🍊🍎🍊🍎🍊🍎🍋🍎🍊🍎🍊🍊🍊🍎🍋🍎🍊🍎
🍋🍎🍊🍎🍊🍎🍊🍎🍊🍎🍊🍎🍊🍊🍎🍋🍎🍊🍎🍊
🍎🍊🍊🍊🍎🍊🍎🍋🍎🍊🍎🍊🍎🍊🍊🍎🍋🍎🍊🍎
🍊🍎🍊🍎🍋🍎🍊🍎🍊🍎🍊🍎🍊🍎🍊🍊🍎🍋🍎🍊
🍎🍋🍎🍊🍊🍊🍎🍊🍎🍊🍊🍎🍊🍎🍊🍎🍋🍎🍊🍎
🍎🍊🍎🍊🍎🍊🍊🍎🍋🍎🍊🍎🍊🍎🍊🍎🍋🍎🍊🍊
🍎🍊🍎🍊🍎🍋🍎🍊🍎🍊🍊🍎🍊🍎🍊🍎🍋🍎🍊🍊
🍎🍊🍊🍎🍊🍊🍎🍋🍎🍊🍎🍋🍎🍊🍎🍊🍊🍎🍊🍎
By matthew on November 16, 2019 at 10:38 PM
Once again, screwed by WordPress code formatting: that should be just ‘print(“🍎🍋🍊”[n],end=””,flush=True)’
By matthew on November 16, 2019 at 10:42 PM
A couple of solutions in Racket.
By r. clayton on November 23, 2019 at 10:31 PM
Another Python solution: this one attempts to rearrange the input list (of 0,1,2) to satisfy the ordering requirements, but making a fairly minimal set of changes. Strategy is to build up result list t, taking first item from s that can be appended to t, also ensuring that a solution is still possible with the remaining items – specifically, that there are enough 1’s and 2’s to keep the 0’s separated (and that if there is only one 0 remaining, it doesn’t get appended when both 1’s and 2’s remain).
def check(n,m,counts): # return True iff m can follow n, given the item counts if n == 0: return m != 0 if n == m: return counts[1]+counts[2] >= counts[0] if m == 0: return counts[0] > 1 or counts[1] == 0 or counts[2] == 0 if n is None: return counts[1]+counts[2] >= counts[0] return False def findnext(s,t,counts): for i in range(len(s)): n = t[-1] if t else None if check(n,s[i],counts): counts[s[i]] -= 1 t.append(s[i]) del s[i] return True return False def isvalid(counts): if counts[0] == 0: return counts[1] == 0 or counts[2] == 0 else: return counts[1]+counts[2] >= counts[0]-1 def f(s): t,counts = [],[0]*3 for i in s: counts[i] += 1 valid = isvalid(counts) while findnext(s,t,counts): pass assert not valid or not s if not s: return "".join("▮○□"[i] for i in t) print(f([0,0,1])) print(f([0,1,2])) print(f([0,0,0,1])) print(f([0,0,0,0,0,0,0,1,2,1,1,1,1])) print(f([1,2,1,1,1,1,0,0,0,0,0,0,0])) print(f([1,2,1,2,1,2,0])) print(f([2,1,2,2,2,2,0,2,0,0,0,1,0,2,0,0,2,0,2,0]))Output (hoping the Unicode chars come out):
By matthew on November 24, 2019 at 1:49 PM
A variation, that takes doesn’t require the mod 3 to already be done:
def check(n,m,counts): n,m = n%3 if n else n,m%3 if n == 0: return m != 0 if n == m: return counts[1]+counts[2] >= counts[0] if m == 0: return counts[0] > 1 or counts[1] == 0 or counts[2] == 0 if n is None: return counts[1]+counts[2] >= counts[0] return False def findnext(s,t,counts): for i in range(len(s)): n = t[-1] if t else None if check(n,s[i],counts): counts[s[i]%3] -= 1 t.append(s[i]) del s[i] return True def f(s): t,counts = [],[0]*3 for i in s: counts[i%3] += 1 while findnext(s,t,counts): pass if not s: return t def show(s): return "".join("▮○□"[i%3] for i in s) if s else None s = [2,67,65,71,86,83,12,17,9,30,42,61,33,68,57,30,74,42,41,51] t = f(s[:]) # Duplicate s print(s) print(t) print(show(s)) print(show(t))Output:
By matthew on November 24, 2019 at 2:33 PM
And here is an in-place version of the above – nice opportunity to use python range assignments & for-else:
def check(n,m,counts): n,m = n%3 if n else n,m%3 if n == 0: return m != 0 if n == m: return counts[1]+counts[2] >= counts[0] if m == 0: return counts[0] > 1 or counts[1] == 0 or counts[2] == 0 if n is None: return counts[1]+counts[2] >= counts[0] return False def g(s): counts = [0]*3 for n in s: counts[n%3] += 1 print([],s) for i in range(len(s)): for j in range(i,len(s)): if check(None if i == 0 else s[i-1],s[j],counts): s[i],s[i+1:j+1] = s[j],s[i:j] counts[s[i]%3] -= 1 break else: return None print(s[:i+1],s[i+1:]) return s g([2,67,65,71,86,83,12,17,9,30,42,61,33,68,57,30,74,42,41,51])This time, we print a trace of the computation:
By matthew on November 24, 2019 at 3:58 PM