Here you will find daily news and tutorials about R, contributed by over 573 bloggers.
There are many ways to follow us - By e-mail:On Facebook: If you are an R blogger yourself you are invited to add your own R content feed to this site (Non-English R bloggers should add themselves- here)

The puzzle in Le Monde this week[end] is called the “square five” (sic!):

Two players each have twenty-five cards with five times each of the digits 1,2,3,4,5. They alternate putting one card on top of the pile, except that they can instead take an arbitrary number of consecutive cards from the top of the pile whose product is a squared number. In which case the corresponding player wins! What is the maximum number of moves and who can find a winning strategy?

At any time, a player cannot put a 1 or a 4 because those are squared numbers. The choice is therefore always among 2,3,5, with a prohibition of squares in the sequence. Testing blindly for the longest possible sequence can be done via the following R code

#testing for squares
whole=function(x,tol = .Machine$double.eps^0.5) abs(x - round(x)) < tol
#Monte Carlo loop
maxl=1
for (t in 1:10^4){
counter=matrix(5,ncol=2,nrow=5)
#random start
rounds=1
pile=sample(c(2,3,5),1)
counter[pile,1]=counter[pile,1]-1
cont=TRUE
playr=2
while (cont){
if (max(whole(sqrt(cumprod(pile[rounds:1]))))==1){
cont=FALSE
}else{
if (sum(counter[c(2,3,5),playr])==0){
newt=4 #game over!
}else{
newt=sample(c(2,3,5)[counter[c(2,3,5),playr]>0],1)
}
pile=c(pile,newt)
counter[newt,playr]=counter[newt,playr]-1
rounds=rounds+1
playr=3-playr
}
}
if (rounds>maxl)
sek=pile
maxl=max(maxl,rounds)
}

which gives a maximum length of 8, a figure that makes sense since picking any sequence of length 8 from 2,3,5 leads to a square by a tree decomposition. It does not seem there is a winning strategy for the second player if the first one plays in an optimal manner…