d13.c

#include "a.h"
#define P 320
#define D(i,j,k,l) a[i] * a[j] - a[k] * a[l]
L f(L a[6], I p)_(I(p, a[4] += 1e13; a[5] += 1e13)
  L m = D(0,3,1,2), c = D(4,3,5,2), d = D(5,0,4,1);
  I(c % m || d % m || c / m < 0 || d / m < 0, R 0); (3 * c + d) / m)
L s[2], i, a[P][6]; M(Fp fp = fget("d13.txt"); LINESF(99, fp, W(*p,
  I(isdigit(*p++), a[i/6][i++%6] = strtol(p-1, &p, 10))))
  Fi(P, Fx(2, s[x] += f(a[i], x))) resl(s[0], s[1]))