Go: package main

import (
	"bufio"
		"os"
	"strings"
	"strconv"
	"fmt"
	"math"
)

var scanner = bufio.NewScanner(os.Stdin)

type Point struct {
	x, y float64
}
func main() {
	for scanner.Scan() {
		buf1 := strings.Split(scanner.Text(), " ")
		line1 := make([]float64, len(buf1))
		for i := range line1 {
			line1[i], _ = strconv.ParseFloat(buf1[i], 64)
		}
		if line1[0] == 0 && line1[1] == 0 &&
			line1[2] == 0 && line1[3] == 0 {
				break
		}
		scanner.Scan()
		buf2 := strings.Split(scanner.Text(), " ")
		line2 := make([]float64, len(buf2))
		for i := range line2 {
			line2[i], _ = strconv.ParseFloat(buf2[i], 64)
		}
		scanner.Scan()
		buf3 := strings.Split(scanner.Text(), " ")
		line3 := make([]float64, len(buf3))
		for i := range line2 {
			line3[i], _ = strconv.ParseFloat(buf3[i], 64)
		}
		if !intersect(line1[0], line1[1], line1[2], line1[3], line2[0], line2[1], line2[2], line2[3]) ||
			!intersect(line2[0], line2[1], line2[2], line2[3], line3[0], line3[1], line3[2], line3[3]) ||
			!intersect(line3[0], line3[1], line3[2], line3[3], line1[0], line1[1], line1[2], line1[3]) {
				fmt.Println("kyo")
		} else {
			p1 := intersection(line1[0], line1[1], line1[2], line1[3], line2[0], line2[1], line2[2], line2[3])
			p2 := intersection(line2[0], line2[1], line2[2], line2[3], line3[0], line3[1], line3[2], line3[3])
			p3 := intersection(line3[0], line3[1], line3[2], line3[3], line1[0], line1[1], line1[2], line1[3])
			area := math.Abs((p2.x - p1.x) * (p3.y - p1.y) - (p2.y - p1.y) * (p3.x - p1.x))/2
			if area < 0.00000001 {
				fmt.Println("kyo")
			} else if 1899999.99999999 < area {
				fmt.Println("dai-kichi")
			} else if 999999.99999999 < area {
				fmt.Println("chu-kichi")
			} else if 99999.99999999 < area {
				fmt.Println("kichi")
			} else {
				fmt.Println("syo-kichi")
			}
		}
	}
}

func cross(x1, y1, x2, y2 float64) int {
	z :=  x1 * y2 - y1 * x2
	if 0 < z {
		return 1
	} else if z < 0 {
		return -1
	} else {
		return 0
	}
}

func intersect(x1, y1, x2, y2, x3, y3, x4, y4 float64) bool {
	if 0 <= cross(x2-x1, y2-y1, x3-x1, y3-y1) *
		cross(x2-x1, y2-y1, x4-x1, y4-y1) {
			return false
	}
	if 0 <= cross(x4-x3, y4-y3, x1-x3, y1-y3) *
		cross(x4-x3, y4-y3, x2-x3, y2-y3) {
			return false
	}
	return true
}

func intersection(x1, y1, x2, y2, x3, y3, x4, y4 float64) Point {
	if x3 == x4 {
		x1, x3 = x3, x1
		y1, y3 = y3, y1
		x2, x4 = x4, x2
		y2, y4 = y4, y2
	}
	ret := Point{}
	if x1 == x2 {
		// y軸に平行な直線の場合
		c := (y4 - y3) / (x4 - x3)
		d := y3 - c * x3
		x := x1
		y := c * x + d
		ret.x, ret.y = x, y
	} else {
		a := (y2 - y1) / (x2 - x1)
		b := y1 - a * x1
		c := (y4 - y3) / (x4 - x3)
		d := y3 - c * x3
		x := (d - b) / (a - c)
		y := a * x + b
		ret.x, ret.y = x, y
	}
	return ret
}

C++: #include<iostream>
#include<cmath>
#include<algorithm>
#include<vector>
#include<cstdio>
#define EPS (1e-10)
#define COUNTER_CLOCKWISE 1
#define CLOCKWISE -1
#define ONLINE_BACK 2
#define ONLINE_FRONT -2
#define ON_SEGMENT 0
using namespace std;


class Point
{
public:
  double x,y;

  Point(double x =0,double y=0):x(x),y(y){}

  Point operator + (Point p){return Point(x+p.x,y+p.y); }
  Point operator - (Point p){return Point(x-p.x,y-p.y); }
  Point operator * (Point p){return Point(x*p.x,y*p.y); }

  bool operator < (const Point &p) const {
    return x != p.x ? x < p.x : y < p.y;
  }  

  bool operator == (const Point &p) const {
    return fabs(x-p.x) < EPS && fabs(y-p.y) < EPS;
  }

};

double norm(Point a)
{
  return a.x*a.x+a.y*a.y;
}

double dot(Point a,Point b)
{
  return a.x*b.x + a.y*b.y;
}

double cross(Point a,Point b)
{
  return a.x*b.y - a.y*b.x;
}

int ccw(Point p0,Point p1,Point p2)
{
  Point a = p1-p0;
  Point b = p2-p0;
  if(cross(a,b) > EPS)return COUNTER_CLOCKWISE;
  if(cross(a,b) < -EPS)return CLOCKWISE;
  if(dot(a,b) < -EPS)return ONLINE_BACK;
  if(norm(a) < norm(b))return ONLINE_FRONT;
  return ON_SEGMENT;
}

bool isIntersect(Point p1,Point p2,Point p3,Point p4)
{
  return (ccw(p1,p2,p3) * ccw(p1,p2,p4) <= 0 &&
          ccw(p3,p4,p1) * ccw(p3,p4,p2) <= 0 );
}

Point getCrossPoint(Point a1,Point a2,Point b1,Point b2)
{
  if(!isIntersect(a1,a2,b1,b2))
    {
      return Point(-9993999,-9993999);
    }
  Point base = b2 - b1;
  double d1 = abs(cross(base,a1-b1));
  double d2 = abs(cross(base,a2-b1));
  double t = d1/(d1+d2);
  return a1 + (a2 - a1)*t;
}

double getArea(Point p1,Point p2,Point p3)
{
  vector<Point> p;
  p.resize(4);
  p[0] = p1,p[1] = p2,p[2] = p3,p[3] = p[0];
  double sum = 0;
  for(int i=0;i<3;i++)
    sum += (p[i].x - p[i+1].x)*(p[i].y+p[i+1].y);
  return fabs(sum/2.0);
}

int main()
{
  int x1,y1,x2,y2;
  while(cin >> x1 >> y1 >> x2 >> y2,(x1|y1|x2|y2))
    {
      Point a(x1,y1);
      Point b(x2,y2);
      cin >> x1 >> y1 >> x2 >> y2;
      Point c(x1,y1);
      Point d(x2,y2);
      cin >> x1 >> y1 >> x2 >> y2;
      Point e(x1,y1);
      Point f(x2,y2);
      Point c1,c2,c3;
      c1 = getCrossPoint(a,b,c,d);
      c2 = getCrossPoint(c,d,e,f);
      c3 = getCrossPoint(e,f,a,b);


      if((c1 == c2 && c2 == c3) || !isIntersect(a,b,c,d) || !isIntersect(c,d,e,f) || !isIntersect(e,f,a,b))
	{
	  cout << "kyo" << endl;
	  continue;
	}
     
      double area = getArea(c1,c2,c3);

      if(0 < area && area < 100000)
	cout << "syo-kichi" << endl;
      else if(100000 <= area && area < 1000000)
	cout << "kichi" << endl;
      else if(1000000 <= area && area < 1900000)
	cout << "chu-kichi" << endl;
      else if(1900000 <= area)
	cout << "dai-kichi" << endl; 

    }
  return 0;
}

C++: #include<iostream>
#include<complex>
#define EPS 1e-10
#define EQ(a,b) (abs((a)-(b))<EPS)
#define EQP(a,b) (EQ((a).real(),(b).real()) && EQ((a).imag(),(b).imag()))
using namespace std;

typedef complex<double> P;

double dot(P a,P b){return real(conj(a)*b);}
double cross(P a,P b){return imag(conj(a)*b);}
double area(P a,P b,P c){return abs(cross(b-a,c-a)/2);}
bool para(P a1,P a2,P b1,P b2){return abs(cross(a1-a2,b1-b2))<EPS;}
bool is_cp(P a1,P a2,P b1,P b2){
  if(cross(a2-a1,b1-a1) * cross(a2-a1,b2-a1) < EPS)
    if(cross(b2-b1,a1-b1) * cross(b2-b1,a2-b1) < EPS)return true;
  return false;
}
P seg_cp(P a1,P a2,P b1,P b2){
  b2-=b1;
  double d = abs(cross(b2,a1-b1));
  return a1 + (a2-a1)*d/( d + abs(cross(b2,a2-b1)) );
}

int main(){
  P p11,p12,p21,p22,p31,p32;
  for(;;){
    cin >> p11.real() >> p11.imag() >> p12.real() >> p12.imag();
    if(!abs(p11) && !abs(p12))break;
    cin >> p21.real() >> p21.imag() >> p22.real() >> p22.imag();
    cin >> p31.real() >> p31.imag() >> p32.real() >> p32.imag();

    if(para(p11,p12,p21,p22) || para(p21,p22,p31,p32) || para(p31,p32,p11,p12))
      cout << "kyo\n";
    else if(!is_cp(p11,p12,p21,p22) || !is_cp(p21,p22,p31,p32) || !is_cp(p31,p32,p11,p12))
      cout << "kyo\n";
    else{
      P cp1 = seg_cp(p11,p12,p21,p22);
      P cp2 = seg_cp(p21,p22,p31,p32);
      P cp3 = seg_cp(p31,p32,p11,p12);
      if(EQP(cp1,cp2))cout << "kyo\n";
      else{
	double s = area(cp1,cp2,cp3);
	if(s+EPS>1900000)cout << "dai-kichi\n";
	else if(s+EPS>1e6)cout << "chu-kichi\n";
	else if(s+EPS>1e5)cout << "kichi\n";
	else cout << "syo-kichi\n";
      }
    }
    
  }
}