Winner the cyclist (typ2)

2006-04-12 09:00:00 UTC

# Cluster-based solver 5c.

Status: Passed
Results: 85367
CPU Time: 78.3
Score: 871.48
Submitted at: 2006-04-12 03:23:01 UTC
Scored at: 2006-04-12 06:21:02 UTC

Current Rank: 1236th
Based on: newbie3 (diff)

Markus Buehren
12 Apr 2006
I have optimized the parameters of a modification of Vivat Profiler over night and put them into newbie3.
Code
```function moves=solver(board,nMoves)

[rows, cols] = size(board);
moves = ones(nMoves,3);
count=0;

% parameters of optimization over night
next_move_count = 26;
if rows*cols >= 350;
choke_factor = 0.922;
depth_factor = 2.82;
rot          = 1.49;
else
choke_factor = 0.92;
depth_factor = 2.97;
rot          = 1.35;
end

board = nthroot(board,rot);

depth = ceil(depth_factor*nMoves/10);
flip = floor(0.66*nMoves);
while 1
% restore original colour values half-way through game
if count==flip, board = nthroot(board,2/(1+rot)); end

% find highest value moves
if isempty(cell_list),
[swapmv,board]=CheckForSwap(board);
if swapmv(1)
count = count+1;
moves(count,:) = swapmv;
if (count==nMoves)
break;
end
continue
end
break
end

% check future move values and re-weight current moves
if count<nMoves-1
[value_list,pos]=sort(value_list,2,'descend');
cell_list=cell_list(pos);
for k=1:min(numel(cell_list),next_move_count)
[new_cell_list,new_value_list] = CalculateMoves(new_board);

if (~isempty(new_cell_list))
depth_ = max(min(depth,nMoves-count-3),1);
if numel(new_value_list)<=depth_
value_list(k) = value_list(k) + sum(new_value_list);
else
new_value_list=sort(new_value_list,2,'descend');
value_list(k) = value_list(k) + sum(new_value_list(1:depth_));
end
end
if (value_list(k) < choke_factor * value_list(1))
break;
end
end
end

% take new most-valuable move
[max_val,pos] = max(value_list); %#ok
max_cell = cell_list(pos);
count = count+1;
moves(count,:) = [mod(max_cell-1,rows)+1,ceil(max_cell/rows),0];

if (count==nMoves)
break;
end
end
moves = moves(1:count,:);

[rows,cols] = size(board);

% find first row from top that has non-zeros
first_row = find(any(board,2),1);

if (first_row == 1)
% first row is non-zero so build column weighted match value
if board(1)
end
k=1;
for j=2:cols
m = k;
k = k + rows;
if board(k)
if (board(m) == board(k))
% block to left matches current
end
end
end
end

% Start at first row (or second if first is top and loop to bottom)
for i = max(2,first_row):rows
k=i;
if (board(k))
% first block in row is non-zero
m = k-1;
if (board(m) == board(k))
% block above matches current
n = -1;
while (n ~= m)
n = m;
end
end
end

% loop through columns and build column+row weighted match value
for j=2:cols
m = k;
k = k + rows;   % move over one column
if (board(k))
if (board(m)==board(k))
% block to left matches current
end
m = k - 1;
if (board(m)==board(k))
% block above matches current
n = -1;
while (n ~= m)
n = m;
end
end
end
end
end

cell_list = zeros(1,256);
group_count = 0;
group_count_matrix = zeros(rows, cols);

for i=rows:-1:first_row
k = i + rows*(cols-1);
for j=cols:-1:1
if board(k)
if (m == k)
group_count_matrix(k) = 1;
else
while group_count_matrix(m) == 0
end
group_block_count = group_count_matrix(m) + 1;
group_count_matrix(m) = group_block_count;
if (group_block_count == 2)
group_count = group_count + 1;
cell_list(1,group_count) = m;
end
end
end
k = k - rows;
end
end

cell_list = cell_list(1:group_count);
value_list = zeros(1,group_count);

for k = 1:group_count
m = cell_list(k);
value_list(k) = group_count_matrix(m) * board(m);
end

function B=ProcessMove(B,dj,y)
R=size(B,1);B(y)=0;
col=ceil(y/R);
minCol=col;maxCol=col;
while dj(y)~=y
y=dj(y);
B(y)=0;
col=ceil(y/R);
minCol=min(minCol, col);
maxCol=max(maxCol, col);
end;
for i=minCol:maxCol,
k=B(B(:,i)>0,i);
B(:,i)=0;B(R-numel(k)+1:R,i)=k;
end

function [swapmv,outboard]=CheckForSwap(board)
[rows cols]=size(board);
swapper=[-1 rows 1 -rows];
lentmp=rows*cols;
swapmv=zeros(1,3);
outboard=board;
maxswap=0;
list=find(board~=0)';
nList=numel(list);
for k=1:nList
jj=list(k);
for swapdir=1:4
np=jj+swapper(swapdir);
if np>0 && np<=lentmp && board(np)>0 && board(np)~=board(jj)
tboard=DoSwap(board,jj,np);
[s_ceil,s_value]=CalculateMoves(tboard); %#ok
if (~isempty(s_value))
if (max(s_value)> maxswap)
maxswap=max(s_value);
idxq=jj;idxr=swapdir;idxs=np;
end
else
for k2=max(1,k-20):min(nList,k+20)
jj2=list(k2);
for swapdir2=1:4
np2=jj2+swapper(swapdir2);
if np2>0 && np2<=lentmp && tboard(np2)>0 && tboard(np2)~=tboard(jj2)
tboard2=DoSwap(tboard,jj2,np2);
if (~isempty(s_value2) && max(s_value2)>maxswap)
isok=false;
[s_max,s_idx]=max(s_value2); %#ok
Y=s_ceil2(s_idx);
if (Y==np || Y==jj)
isok=true;
end
end
if (Y==np || Y==jj)
isok=true;
end
if (isok)
maxswap=max(s_value2);
idxq=jj;idxr=swapdir;idxs=np;
end
end
end
end
end
end
end
end
end
if (maxswap>0)
swapmv = [mod(idxq-1,rows)+1,ceil(idxq/rows),idxr];
outboard=DoSwap(board,idxq,idxs);
end

function board = DoSwap(board,el1,el2)
board([el1,el2]) = board([el2,el1]);
```