%% Source code:
c=@(x,p)char(mod(idivide(x,1000*26^p),26)+65);
s=@(n)[c(n,3),c(n,2),num2str(mod(n,1000),'-%03d-'),c(n,1),c(n,0)];
%% chose some test cases
n = uint32([0;1;999;1000;675999;676000;456975999]) ;
% work out their plate numbers
s( n )
JSUgU291cmNlIGNvZGU6CmM9QCh4LHApY2hhcihtb2QoaWRpdmlkZSh4LDEwMDAqMjZecCksMjYpKzY1KTsKcz1AKG4pW2MobiwzKSxjKG4sMiksbnVtMnN0cihtb2QobiwxMDAwKSwnLSUwM2QtJyksYyhuLDEpLGMobiwwKV07CgolJSBjaG9zZSBzb21lIHRlc3QgY2FzZXMKbiA9IHVpbnQzMihbMDsxOzk5OTsxMDAwOzY3NTk5OTs2NzYwMDA7NDU2OTc1OTk5XSkgOwoKJSB3b3JrIG91dCB0aGVpciBwbGF0ZSBudW1iZXJzCnMoIG4gKQo=