mscroggs.co.uk
mscroggs.co.uk

subscribe

Comment

Comments

Comments in green were written by me. Comments in blue were not written by me.
@Lord Sméagol: Hello and happy New Year!
9 minutes is cool!
The answer is wrong because of _lzcnt instruction, as you suspected, as turns out it works differently on different cpus: https://nextmovesoftware.com/blog/2017...
With this error, solutions having 1x1 square directly in the center are not counted.

I guess, gcc/clang do it correctly because I specify -march=native (so it checks cpu and generates correct instruction), and run where I compile. But it's a potential problem I probably need to add some assertions to the code.

Maybe on your hardware you can either use WSL and clang compiler, or set constexpr bool USE_SSE_QUADRANT_FILL=false, to fall back to slower.
You could also try to use BitScanReverse instead of __lzcnt, but it has different input/output so I'm not sure how hard would that be to fix it.
Oleg
on /blog/119
               
@Oleg:

The includes got filtered:

Util.h

//#include [bits/stdc++.h] // not including all
#include [filesystem] // just include what's needed
#include [array] // just include what's needed
#include [mutex] // just include what's needed

:)
Lord Sméagol
on /blog/119
               
@Oleg:

I removed my macros:
#define __tzcnt_u32(v) ((v) ? (_tzcnt_u32(v)) : (32))
#define __lzcnt32(v) ((v) ? (_lzcnt_u32(v)) : (32))
replacing them with simple inline code


Util.h

//#include // not including all
#include // just include what's needed
#include // just include what's needed
#include // just include what's needed

#if 1 // use safe localtime
struct tm buf; // use safe localtime
auto err = localtime_s(&buf, &cur_time); // use safe localtime
return std::put_time(&buf, "%F %T"); // use safe localtime
#else // use safe localtime
return std::put_time(std::localtime(&cur_time), "%F %T");
#endif // use safe localtime


State.h

changed _mm_set_epi8(0x80 to -0x80 to stop warnings

inline replacement:
//int i = __tzcnt_u32(mask); // for no BMI; without zero test, as not needed here
int i = _tzcnt_u32(mask); // for no BMI; without zero test, as not needed here

inline replacement:
//int last_idx_before_mid = 31 - __lzcnt32(off_mask); // for no BMI; without zero test, as not needed here
int last_idx_before_mid = _lzcnt_u32(off_mask); // for no BMI; without zero test, as not needed here


Solver.h

inline replacement:
//return ini.size(); // to stop warning
return (int)ini.size(); // to stop warning

inline replacement:
//const int dim = __tzcnt_u32(mask); // for no BMI; without zero test, as not needed here
const int dim = _tzcnt_u32(mask); // for no BMI; without zero test, as not needed here


I tried '9' runs: with asserts: 10:31, without: 10:18 (saved 2%)
A minute slower than the faulty version, but still not too bad for a 2013 (Q3) CPU :)
Lord Sméagol
on /blog/119
               
@Oleg: Happy new year!

I just added this:

#if 0
int last_idx_before_mid = 31 - __lzcnt32(off_mask); // 31 - LZCNT ==> index of MSb
#else

// if off_mask can never be zero, no need for check to override BSR result
assert(off_mask);
// a '9' run didn't reveal any 0 [you would know for sure for other sizes]

// need unsigned long result
unsigned long last_idx_before_mid;

// get index of MSb [no need for adjustment if off_mask can never be zero]
_BitScanReverse(&last_idx_before_mid, off_mask);
#endif

a run of '9' now produces the correct result: 1,730,280 :)
Lord Sméagol
on /blog/119
               

Archive

Show me a random blog post
 2026 

May 2026

World Cup stickers 2026

Apr 2026

A new puzzle every day
Mixing Wordle with other games

Feb 2026

Christmas (2025) is over
 2025 

Dec 2025

Christmas card 2025

Nov 2025

Christmas (2025) is coming!

Sep 2025

The partridge puzzle

Aug 2025

TMiP 2025 puzzle hunt

Jun 2025

A nonogram alphabet

Mar 2025

How to write a crossnumber

Jan 2025

Christmas (2024) is over
Friendly squares
 2024 

Dec 2024

A regular expression Christmas puzzle
Christmas card 2024

Nov 2024

Christmas (2024) is coming!

Feb 2024

Zines, pt. 2

Jan 2024

Christmas (2023) is over
 2023 
▼ show ▼
 2022 
▼ show ▼
 2021 
▼ show ▼
 2020 
▼ show ▼
 2019 
▼ show ▼
 2018 
▼ show ▼
 2017 
▼ show ▼
 2016 
▼ show ▼
 2015 
▼ show ▼
 2014 
▼ show ▼
 2013 
▼ show ▼
 2012 
▼ show ▼

Tags

preconditioning pizza cutting approximation a gamut of games 24 hour maths puzzles christmas javascript golden ratio oeis bempp games standard deviation advent calendar crosswords electromagnetic field royal institution partridge puzzle pythagoras wave scattering nonograms mathsjam nine men's morris coins gerry anderson pi bodmas sobolev spaces countdown data visualisation books game of life data determinants anscombe's quartet crossnumbers kenilworth ternary final fantasy realhats simultaneous equations signorini conditions manchester science festival estimation folding tube maps matrix of cofactors london reddit crochet arrangement puzzles reuleaux polygons gaussian elimination zines weather station people maths football london underground stirling numbers thirteen manchester edinburgh dates inline code speed bluesky statistics php finite group european cup rugby ucl hats matrix of minors mathslogicbot fonts plastic ratio pac-man captain scarlet tennis noughts and crosses inverse matrices hannah fry triangles rust dragon curves golden spiral numbers wool logo interpolation binary geogebra world cup asteroids weak imposition turtles pi approximation day bots dinosaurs recursion chebyshev coventry polynomials light go platonic solids arithmetic crossnumber cross stitch craft big internet math-off mathsteroids braiding talking maths in public trigonometry tetris christmas card runge's phenomenon youtube sport menace computational complexity misleading statistics stickers finite element method datasaurus dozen matt parker squares cambridge error bars bubble bobble frobel guest posts pokémon raspberry pi palindromes regular expressions dataset warwick newcastle martin gardner news correlation harriss spiral draughts chalkdust magazine sound friendly squares royal baby logic hyperbolic surfaces errors fractals matrix multiplication programming phd python logs matrices folding paper video games national lottery kings sorting propositional calculus boundary element methods radio 4 rhombicuboctahedron gather town wordle geometry map projections numerical analysis curvature flexagons quadrilaterals the aperiodical accuracy chess pokémon wordle mean databet probability game show probability alphabets tmip pascal's triangle convergence hexapawn machine learning graphs exponential growth graph theory fence posts live stream latex

Archive

Show me a random blog post
▼ show ▼
© Matthew Scroggs 2012–2026