function a = chain(n, a)
if nargin == 1
a = [];
end
if n <= 1
a = [a n];
elseif mod(n, 2)
a = chain(n * 3 + 1, [a n]);
else
a = chain(n / 2, [a n]);
end
end
f884 = @(a, b, c) sum(c <= arrayfun(@(n) numel(chain(n)), a:b));
chain(10), chain(30), f884(1, 100, 15)
ZnVuY3Rpb24gYSA9IGNoYWluKG4sIGEpCglpZiBuYXJnaW4gPT0gMQoJCWEgPSBbXTsKCWVuZAoJaWYgbiA8PSAxCgkJYSA9IFthIG5dOwoJZWxzZWlmIG1vZChuLCAyKQoJCWEgPSBjaGFpbihuICogMyArIDEsIFthIG5dKTsKCWVsc2UKCQlhID0gY2hhaW4obiAvIDIsIFthIG5dKTsKCWVuZAplbmQKZjg4NCA9IEAoYSwgYiwgYykgc3VtKGMgPD0gYXJyYXlmdW4oQChuKSBudW1lbChjaGFpbihuKSksIGE6YikpOwpjaGFpbigxMCksIGNoYWluKDMwKSwgZjg4NCgxLCAxMDAsIDE1KQ==