Publications by xi'an

Le Monde puzzle [#1139]

23.04.2020

A weekly Monde current mathematical puzzle that reminded me of an earlier one (but was too lazy to check): The integer n=36 enjoys the property that all the differences between its ordered divisors are also divisors of 36. Find the only 18≤m≤100 that enjoys this property such that all its prime dividers areof multiplicity one. Are there othe...

956 sym R (124 sym/1 pcs) 2 img

another surmortaliy graph

27.04.2020

Another graph showing the recent peak in daily deaths throughout France as recorded by INSEE and plotted by Baptiste Coulmont from Paris 8 Sociology Department. And further discussed by Arthur Charpentier on Freakonometrics. With a few days off due to reporting, this brings an objective perspective on the impact of the epidemics (and of the quara...

933 sym 2 img

Le Monde puzzle [#1141]

03.05.2020

The weekly puzzle from Le Monde is in honour of John Conway, who just passed away, ending up his own game of life: On an 8×8 checker-board, Alice picks n squares as “infected”. She then propagates the disease by having each square with least two infected neighbours to become infected as well. What is the minimal value of n for the entire bo...

1447 sym 8 img

more games of life

04.05.2020

Another puzzle in memoriam of John Conway in The Guardian: Find the ten digit number, abcdefghij. Each of the digits is different, and a is divisible by 1 ab is divisible by 2 abc is divisible by 3 abcd is divisible by 4 abcde is divisible by 5 abcdef is divisible by 6 abcdefg is divisible by 7 abcdefgh is divisible by 8 abcdefghi is divisible b...

964 sym R (408 sym/2 pcs) 2 img

another Bernoulli factory

17.05.2020

A question that came out on X validated is asking for help in figuring out the UMVUE (uniformly minimal variance unbiased estimator) of (1-θ)½ when observing iid Bernoulli B(θ). As it happens, there is no unbiased estimator of this quantity and hence not UMVUE. But there exists a Bernoulli factory producing a coin with probability (1-θ)½ fro...

1559 sym 2 img

not a Bernoulli factory

19.05.2020

A Riddler riddle I possibly misunderstood: Four isolated persons are given four fair coins, which can be either flipped once or returned without being flipped. If all flipped coins come up heads, the team wins! Else, if any comes up tails, or if no flip at all is done, it looses. Each person is further given an independent U(0,1) realisation. Wh...

2329 sym 4 img

a new Monty Hall riddle

21.05.2020

The Riddler was sort of feeling the rising boredom of being under lockdown when proposing the following variant to the Monty Hall puzzle: There are zero to three goats, with a probability ¼ each, and they are allocated to different doors uniformly among the three doors of the show. After the player chooses a door, Monty opens another door hiddi...

1135 sym Python (64 sym/1 pcs) 2 img

minimax, maximin or plain

31.05.2020

A simple riddle from The Riddler on choosing between the maximum between two minima of two throws of an N-face dice, the minimum between two maxima of two throws of an N-face dice, and a single throw. Since maximin is always less than maximin, second choice is always worse than first and a stochastic domination version of the argument shows the s...

784 sym 2 img

Le Monde puzzle [#1146]

04.06.2020

The weekly puzzle from Le Monde is once more disappointing. Everyday of the month, take 0, 1 or 2 units. If one unit taken past day, next day none can be taken. If two units taken two day ago, none can be taken the current day. What is the strategy maximising the number of units  for February? March? April? Generalise to 0,..,3 units taken each...

1638 sym R (443 sym/2 pcs) 2 img

the strange occurrence of the one bump

07.06.2020

When answering an X validated question on running an accept-reject algorithm for the Gamma distribution by using a mixture of Beta and drifted (bt 1) Exponential distributions, I came across the above glitch in the fit of my 10⁷ simulated sample to the target, apparently displaying a wrong proportion of simulations above (or below) one. a=.9 g...

959 sym R (192 sym/2 pcs) 4 img