Permutationen, bei denen kein k + 2 Punkte auf ein Polynom vom Grad k fallen

16

Beschreibung

Lassen Sie eine Permutation der ganzen Zahlen {1, 2, ..., n}als minimal interpolierbar bezeichnet werden, wenn keine k+2Punktmenge (zusammen mit ihren Indizes) auf ein Polynom vom Grad fällt k. Das ist,

  1. Keine zwei Punkte fallen auf eine horizontale Linie (0-Grad-Polynom)
  2. Es fallen keine drei Punkte auf eine Linie (1-Grad-Polynom)
  3. Keine vier Punkte fallen auf eine Parabel (2-Grad-Polynom)
  4. Und so weiter.

Herausforderung

Ein Programm schreiben , die OEIS Sequenz berechnet A301802 (n) , die Anzahl der Permutationen von minimal interpolierbar {1, 2, ..., n}für nals groß wie möglich.


Wertung

Ich werde Ihren Code auf meinem Computer (2,3 GHz Intel Core i5, 8 GB RAM) mit zunehmenden Eingaben timen. Ihre Punktzahl ist die größte Eingabe, die weniger als 1 Minute dauert, um den richtigen Wert auszugeben.


Beispiel

Beispielsweise ist die Permutation [1, 2, 4, 3]wegen minimal interpolierbar

the terms together with their indices 
[(1, 1), (2, 2), (3, 4), (4, 3)] 
have the property that
  (0) No two points have the same y-value.
  (1) No three points lie on a line.
  (2) No four points lie on a parabola.

Beispiel zur Veranschaulichung, dass [1,2,4,3] minimal interpolierbar ist. In der Abbildung ist zu sehen, dass die horizontalen Linien (rot) höchstens einen Punkt aufweisen, die Linien (blau) höchstens zwei Punkte aufweisen und die Parabeln (grün) drei Punkte aufweisen.


Daten

Hier sind die minimal interpolierbar Permutationen für n=3, n=4und n=5:

n = 3: [1,3,2],[2,1,3],[2,3,1],[3,1,2]
n = 4: [1,2,4,3],[1,3,2,4],[1,3,4,2],[1,4,2,3],[2,1,3,4],[2,1,4,3],[2,3,1,4],[2,4,1,3],[2,4,3,1],[3,1,2,4],[3,1,4,2],[3,2,4,1],[3,4,1,2],[3,4,2,1],[4,1,3,2],[4,2,1,3],[4,2,3,1],[4,3,1,2]
n = 5: [1,2,5,3,4],[1,3,2,5,4],[1,3,4,2,5],[1,4,2,3,5],[1,4,3,5,2],[1,4,5,2,3],[1,4,5,3,2],[1,5,3,2,4],[2,1,4,3,5],[2,3,1,4,5],[2,3,5,1,4],[2,3,5,4,1],[2,4,1,5,3],[2,4,3,1,5],[2,4,5,1,3],[2,5,1,3,4],[2,5,1,4,3],[2,5,3,4,1],[2,5,4,1,3],[3,1,4,5,2],[3,1,5,2,4],[3,1,5,4,2],[3,2,5,1,4],[3,2,5,4,1],[3,4,1,2,5],[3,4,1,5,2],[3,5,1,2,4],[3,5,1,4,2],[3,5,2,1,4],[4,1,2,5,3],[4,1,3,2,5],[4,1,5,2,3],[4,1,5,3,2],[4,2,1,5,3],[4,2,3,5,1],[4,2,5,1,3],[4,3,1,2,5],[4,3,1,5,2],[4,3,5,2,1],[4,5,2,3,1],[5,1,3,4,2],[5,2,1,3,4],[5,2,1,4,3],[5,2,3,1,4],[5,2,4,3,1],[5,3,2,4,1],[5,3,4,1,2],[5,4,1,3,2]

Wenn mein Programm korrekt ist, die ersten paar Werte von a(n), die Anzahl der minimal interpolierbaren Permutationen von {1, 2, ..., n}:

a(1) = 1
a(2) = 2
a(3) = 4
a(4) = 18
a(5) = 48
a(6) = 216
a(7) = 584
a(8) = 2870
Peter Kagey
quelle
Schöne Sequenznummer! | Obwohl Sie den schnellsten Code angegeben haben , haben Sie nicht angegeben, auf welchem ​​Computer er am schnellsten ist. Was ist genau das Gewinnkriterium?
user202729
3
Um den Kommentar von user202729 zu ergänzen, schlage ich einige Tags vor, mit denen Sie die Gewinnkriterien bestimmen können: Für den schnellste Code müssen die Einsendungen auf demselben Computer getestet werden, um die Laufzeit zu vergleichen (in der Regel wird dies vom OP der Challenge durchgeführt). Der schnellste Algorithmus fordert die Antwortenden auf, Code mit der geringstmöglichen zeitlichen Komplexität zu entwickeln. Code-Golf fordert die Benutzer auf, Code mit dem kürzesten Quellcode (oder einem gleichwertigen Code) zu entwickeln. Davon abgesehen ist dies in der Tat eine schöne Herausforderung.
JungHwan Min
Ihr Beispieltext verwendet eine Indexierung von Null, obwohl das Bild eine Indexierung verwendet.
Jonathan Frech
Ist es nicht unmöglich, dass zwei Punkte dieselbe Höhe einnehmen, da alle Punkte durch Permuationen der ersten natürlichen Zahlen definiert sind?
Jonathan Frech
@ JonathanFrech, in der Tat sollte es 1-indiziert sein, da dies Permutationen sind. Und du hast recht! Da es sich um Permutationen handelt, ist die 0-Grad-Polynombedingung kostenlos.
Peter Kagey

Antworten:

5

C #

using System;
using System.Diagnostics;
using BigInteger = System.Int32;

namespace Sandbox
{
    class PPCG160382
    {
        public static void Main(params string[] args)
        {
            if (args.Length != 0)
            {
                foreach (var arg in args) Console.WriteLine(CountValidPerms(int.Parse(arg)));
            }
            else
            {
                int[] smallValues = new int[] { 1, 1, 2, 4, 18, 48 };
                for (int n = 0; n < smallValues.Length; n++)
                {
                    var observed = CountValidPerms(n);
                    var expected = smallValues[n];
                    Console.WriteLine(observed == expected ? $"{n}: Ok" : $"{n}: expected {expected}, observed {observed}, error {observed - expected}");
                }
                for (int n = smallValues.Length; n < 13; n++)
                {
                    Stopwatch sw = new Stopwatch();
                    sw.Start();
                    Console.WriteLine($"{n}: {CountValidPerms(n)} in {sw.ElapsedMilliseconds}ms");
                }
            }
        }

        private static long CountValidPerms(int n)
        {
            // We work on the basis of exclusion by extrapolation.
            var unused = (1 << n) - 1;
            var excluded = new int[n];
            int[] perm = new int[n];

            // Symmetry exclusion: perm[0] < (n+1) / 2
            if (n > 1) excluded[0] = (1 << n) - (1 << ((n + 1) / 2));

            long count = 0;
            CountValidPerms(ref count, perm, 0, unused, excluded);
            return count;
        }

        private static void CountValidPerms(ref long count, int[] perm, int off, int unused, int[] excluded)
        {
            int n = perm.Length;
            if (off == n)
            {
                count += CountSymmetries(perm);
                return;
            }

            // Quick-aborts
            var completelyExcluded = excluded[off];
            for (int i = off + 1; i < n; i++)
            {
                if ((unused & ~excluded[i]) == 0) return;
                completelyExcluded &= excluded[i];
            }
            if ((unused & completelyExcluded) != 0) return;

            // Consider each unused non-excluded value as a candidate for perm[off]
            var candidates = unused & ~excluded[off];
            for (int val = 0; candidates > 0; val++, candidates >>= 1)
            {
                if ((candidates & 1) == 0) continue;

                perm[off] = val;

                var nextUnused = unused & ~(1 << val);

                var nextExcluded = (int[])excluded.Clone();
                // For each (non-trivial) subset of smaller indices, combine with off and extrapolate to off+1 ... excluded.Length-1
                if (off < n - 1 && off > 0)
                {
                    var points = new Point[off + 1];
                    var denoms = new BigInteger[off + 1];
                    points[0] = new Point { X = off, Y = perm[off] };
                    denoms[0] = 1;
                    ExtendExclusions(perm, off, 0, points, 1, denoms, nextExcluded);
                }

                // Symmetry exclusion: perm[0] < perm[-1] < n - 1 - perm[0]
                if (off == 0 && n > 1)
                {
                    nextExcluded[n - 1] |= (1 << n) - (2 << (n - 1 - val));
                    nextExcluded[n - 1] |= (2 << val) - 1;
                }

                CountValidPerms(ref count, perm, off + 1, nextUnused, nextExcluded);
            }
        }

        private static void ExtendExclusions(int[] perm, int off, int idx, Point[] points, int numPoints, BigInteger[] denoms, int[] excluded)
        {
            if (idx == off) return;

            // Subsets without
            ExtendExclusions(perm, off, idx + 1, points, numPoints, denoms, excluded);

            // Just add this to the subset
            points[numPoints] = new Point { X = idx, Y = perm[idx] };
            denoms = (BigInteger[])denoms.Clone();
            // Update invariant: denoms[s] = prod_{t != s} points[s].X - points[t].X
            denoms[numPoints] = 1;
            for (int s = 0; s < numPoints; s++)
            {
                denoms[s] *= points[s].X - points[numPoints].X;
                denoms[numPoints] *= points[numPoints].X - points[s].X;
            }
            numPoints++;

            for (int target = off + 1; target < excluded.Length; target++)
            {
                BigInteger prod = 1;
                for (int t = 0; t < numPoints; t++) prod *= target - points[t].X;

                Rational sum = new Rational(0, 1);
                for (int s = 0; s < numPoints; s++) sum += new Rational(prod / (target - points[s].X) * points[s].Y, denoms[s]);

                if (sum.Denom == 1 && sum.Num >= 0 && sum.Num < excluded.Length) excluded[target] |= 1 << (int)sum.Num;
            }

            // Subsets with
            ExtendExclusions(perm, off, idx + 1, points, numPoints, denoms, excluded);
        }

        private static int CountSymmetries(int[] perm)
        {
            if (perm.Length < 2) return 1;

            int cmp = 0;
            for (int i = 0, j = perm.Length - 1; i <= j; i++, j--)
            {
                cmp = perm.Length - 1 - perm[i] - perm[j];
                if (cmp != 0) break;
            }

            return cmp > 0 ? 4 : cmp == 0 ? 2 : 0;
        }

        public struct Point
        {
            public int X;
            public int Y;
        }

        public struct Rational
        {
            public Rational(BigInteger num, BigInteger denom)
            {
                if (denom == 0) throw new ArgumentOutOfRangeException(nameof(denom));

                if (denom < 0) { num = -num; denom = -denom; }

                var g = _Gcd(num, denom);
                Num = num / g;
                Denom = denom / g;
            }

            private static BigInteger _Gcd(BigInteger a, BigInteger b)
            {
                if (a < 0) a = -a;
                if (b < 0) b = -b;
                while (a != 0)
                {
                    var tmp = b % a;
                    b = a;
                    a = tmp;
                }
                return b;
            }

            public BigInteger Num;
            public BigInteger Denom;

            public static Rational operator +(Rational a, Rational b) => new Rational(a.Num * b.Denom + a.Denom * b.Num, a.Denom * b.Denom);
        }
    }
}

Nimmt Werte von nals Kommandozeilenargumente, oder wenn ohne Argumente ausgeführt, mal sich bis zu n=10. Als "Release" in VS 2017 kompiliert und auf einem Intel Core i7-6700 ausgeführt, rechne ich n=9in 1,2 Sekunden und n=10in 13,6 Sekunden. n=11ist etwas mehr als 2 Minuten.

FWIW:

n    a(n)
9    10408
10   45244
11   160248
12   762554
Peter Taylor
quelle