/*
* A160426.c
* Calculates the entries of Sloane's A160426
*
* Created by Nathaniel Johnston (nathaniel@nathanieljohnston.com)
* March 30, 2011
*/

#include <math.h>
#include <stdio.h>


unsigned int prevPts = 5, numPts = 5, pt[268435455];
FILE *file;

unsigned int getCoord (unsigned int pt, unsigned short coord)
{
    unsigned short shift = 14*(coord-1);
    return (pt & (16383*(1 << shift))) >> shift;
}

unsigned char getOrien (unsigned int pt)
{
    return (pt & (1 << 28)) >> 28;
}

unsigned int getPt (unsigned short x, unsigned short y, unsigned short orien)
{
    return x | (y << 14) | (orien << 28);
}

unsigned char addOrien (unsigned short x, unsigned short y, unsigned short orien)
{
    unsigned int j, nn = 0;
    unsigned short px, py, po;
    for(j=0;j<prevPts;j++)
    {
        px = getCoord(pt[j],1);
        py = getCoord(pt[j],2);
        po = getOrien(pt[j]);

        if(po == 1 && ((px == x && py == y+1) || (px == x && py == y-1))){
            if(orien == 1)return 0;nn++;
        }
        else if(po == 0 && ((px == x-1 && py == y) || (px == x+1 && py == y))){
            if(orien == 0)return 0;nn++;
        }
    }
    return (nn == 1);
}

unsigned char isOccupied (unsigned short x, unsigned short y)
{
    unsigned int j;
    for(j=0;j<prevPts;j++)
    {
        if(x == getCoord(pt[j],1) && y == getCoord(pt[j],2))return 1;
    }
    return 0;
}

int main ()
{
    unsigned short xi, yi, oi, tOr, gen, maxGens;
    pt[0] = getPt(8192,8192,0);
    pt[1] = getPt(8194,8192,0);
    pt[2] = getPt(8196,8192,0);
    pt[3] = getPt(8195,8192,1);
    pt[4] = getPt(8195,8190,1);

    printf("This tool will calculate the entries of Sloane's A160426.\nPlease enter the number of terms to compute (an integer from 1 to 8000): ");
    scanf("%d",&maxGens);

    file = fopen("b160426.txt","w");
    fprintf(file,"0 0\n1 5\n");
    fclose(file);
    printf("0 0\n1 5\n");

    for(gen=1;gen<=maxGens;gen++){
        for(xi=8191-(int)(gen/2);xi<=8197+(int)(gen/2);xi++){
            for(yi=8189-(int)(gen/2);yi<=8193+(int)(gen/2);yi++){
                tOr = isOccupied(xi,yi);
                for(oi=0;oi<2;oi++){
                    if(!tOr && addOrien(xi,yi,oi))pt[numPts++]=getPt(xi,yi,oi);
                }
            }
        }
        file = fopen("b160426.txt","a");
        fprintf(file,"%d %d\n",gen+1,numPts);
        fclose(file);
        printf("%d %d\n",gen+1,numPts);
        prevPts = numPts;
    }

    printf("Done!");
    getchar();
    return 0;
}