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

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

Archive

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