t=zeros(1,2);
for ii=1:10
xx = randi( 500, 1, 1000 );
yy = randi( 500, 1, 1000 );
tic;
s1 = sum( A( sub2ind( size(A), yy, xx ) ) == 0 );
t(1) = t(1)+toc;
tic;
s2 = sum( A( [1 size(A,1)] * ( [yy; xx] - 1 ) + 1 ) == 0 );
t(2) = t(2)+toc;
end
t
QSA9IHJhbmQoIDUwMCApID4gLjU7Cgp0PXplcm9zKDEsMik7CmZvciBpaT0xOjEwCgl4eCA9IHJhbmRpKCA1MDAsIDEsIDEwMDAgKTsKCXl5ID0gcmFuZGkoIDUwMCwgMSwgMTAwMCApOwoKCXRpYzsKCXMxID0gc3VtKCBBKCBzdWIyaW5kKCBzaXplKEEpLCB5eSwgeHggKSApID09IDAgKTsKCXQoMSkgPSB0KDEpK3RvYzsKCQoJdGljOwoJczIgPSBzdW0oIEEoIFsxIHNpemUoQSwxKV0gKiAoIFt5eTsgeHhdIC0gMSApICsgMSApID09IDAgKTsKCXQoMikgPSB0KDIpK3RvYzsKCWFzc2VydChzMT09czIpOwplbmQKdA==