SET VARIABLE start_year = '2000';
SET VARIABLE end_year = '2030-12-31'::date;
SET VARIABLE first_jan = concat(getvariable('start_year'), '-01-01')::date;
SET VARIABLE first_march = concat(getvariable('start_year'), '-03-01')::date;

WITH t_january AS (
    SELECT
        YEAR(generate_series) AS current_year,
        DAYOFWEEK(generate_series) AS wday
    FROM generate_series(
        getvariable('first_jan'),
        getvariable('end_year'),
        INTERVAL 1 year
    )
),
t_march AS (
    SELECT
        YEAR(generate_series) AS current_year,
        DAYOFWEEK(generate_series) AS wday
    FROM generate_series(
        getvariable('first_march'),
        getvariable('end_year'),
        INTERVAL 1 year
    )
),
key_by_year AS (
    SELECT
        t_january.current_year,
        concat(t_january.wday, t_march.wday) AS wday_key
    FROM t_january
    JOIN t_march
        ON t_january.current_year = t_march.current_year
),
calendars_joined AS (
    SELECT
        k1.current_year,
        k1.wday_key,
        k2.current_year AS compatible_year,
    FROM key_by_year AS k1
    JOIN key_by_year AS k2
        ON k1.wday_key = k2.wday_key
        AND k1.current_year != k2.current_year
)
SELECT
    current_year,
    array_agg(compatible_year ORDER BY compatible_year) AS same_wday_years
FROM calendars_joined
GROUP BY current_year
ORDER BY current_year;