Problem:
最後の扉

秘境の地"ナラマンダ地方"

この奥深い洞窟に眠る黄金の伝説を求め、探検家ボブとマイクは様々な困難に立ち向かいました。ようやくたどり着いた最後の扉。この先には夢に見た黄金が待っています。

その扉には多数の三角形が描かれていました。試しにマイクがその 1 つに触れてみると、その三角形は光を放ちました。光は三角形の底辺から垂直に頂点方向へ長方形を描きながら伸びていき、その先にある別の三角形に触れました。すると、光に触れた三角形も同じように光を放ち、それを繰り返していきました。

扉の三角形は全て二等辺三角形(正三角形は含まない)で、どれも同じ大きさで、どれも重ならず（かつ触れることなく）、またどれも光は一定の距離だけ伸びていました。最後の三角形が光り終えると、そこには何事もなかったかのような静けさが戻ってきました。

ボブは自分の足下に奇妙なメッセージを見つけました。

「むやみに触れるべからず。汝正しくこの扉に触れたとき、この扉は放たれん。 」

ボブとマイクはこのメッセージの意味を考えながら、 何度も三角形に触れてみましたがなかなか扉は開きません。 しかし、それを繰り返していくうちに、二人はメッセージの意味がだんだん分かってきました。

「"むやみに触れるべからず"とは、扉に触れる回数は最小ということ。

できるだけ少ない回数で触れて、全ての三角形を光らせれば、きっと扉は開くはずだ!」

例えば、扉には図 1 のように複数の二等辺三角形が描かれています。

図 1

最初に左下の三角形に触れると、図 2 に示すように、その三角形の底辺から垂直に頂点方向へ光が放たれ、その光が触れる三角形が次々に光ります。

図 2

次に左端の三角形に触れると、図 3 に示すように、全ての三角形が光ります。この手順が全ての三角形を光らせるための最少の手数(2 手)となります。

図 3

さあ、ボブとマイクが最後の扉を開ける手助けをしよう。

扉に描かれている三角形の個数、光の伸びる長さ、それぞれの三角形の頂点の座標を入力とし、全ての三角形を光らせるために扉に触れる最少の回数を出力するプログラムを作成して下さい。光の伸びる長さは、底辺からの長さであるものとします。また、光はそれを放つ三角形の頂点を含むような十分な長さを持ちます。

なお、この問題を解くにあたっては、2 点の距離が 0.01 以下の場合には、同一の点と処理します。また、1 点とひとつの直線の距離が 0.01 以下の場合には、この点はこの直線上にあるとして処理します。

Input

複数のデータセットの並びが入力として与えらる。入力の終わりはゼロふたつの行で示される。各データセットは以下の形式で与えらる。

n d
T1
T2
:
Tn

１行目に扉に描かれている三角形の個数 n (1 ≤ n ≤ 100) と光の伸びる長さを表す整数 d (1 ≤ d ≤ 1000) が与えられる。

続く n 行に各三角形の情報 Ti が以下の形式で与えられる。

x1 y1 x2 y2 x3 y3

三角形の３つの頂点の座標が整数で与えられる(-1000 ≤ x1, y1, x2, y2, x3, y3 ≤ 1000)。

データセットの数は 20 を超えない。

Output

データセットごとに、扉に触れる最少の回数を1行に出力する。

Sample Input

3 4
1 0 3 0 2 1
2 3 2 5 3 4
5 3 5 5 6 4
3 2
1 0 3 0 2 1
2 3 2 5 3 4
5 3 5 5 6 4
0 0

Output for the Sample Input

1
3

Solution in C++ (judge verdict: Wrong Answer)

#include<iostream>
#include<vector>
#include<complex>
#include<cstring>
#include<algorithm>

#define rep(i,a) for(int i=0;i<a;i++)
#define EQ(a,b) (abs((a)-(b))<EPS)
#define pb push_back
#define fs first
#define sc second
using namespace std;

typedef double D;
typedef complex<D> P;
typedef pair<P,P> L;
typedef vector<P> Poly;

namespace std{
  bool operator<(const P &a,const P &b){
    return real(a)==real(b)?imag(a)<imag(b):real(a)<real(b);
  }
}

const D EPS = 1e-8;

vector<int> g[110];
vector<int> rg[110];
vector<int> vs;
bool use[110];
int ord[110];
int n,d;
Poly tri[110];
Poly rect[110];

P unit(P p){return p/abs(p);}
pair<P,P> norm(P p){return make_pair(p*P(0,1),p*P(0,-1));}

D dot(P x,P y){return real(conj(x)*y);}
D cross(P x,P y){return imag(conj(x)*y);}

vector<P> convex_hull(vector<P> v){
  int n = v.size(), k = 0;
  sort(v.begin(),v.end());
  vector<P> r(2*n);
  for(int i=0;i<n;i++){
    while(k>1 && cross(r[k-1]-r[k-2],v[i]-r[k-1]) <= EPS)k--;
    r[k++] = v[i];
  }
  for(int i=n-2,t=k;i>=0;i--){
    while(k>t && cross(r[k-1]-r[k-2],v[i]-r[k-1]) <= EPS)k--;
    r[k++] = v[i];
  }
  r.resize(k-1);
  return r;
}

D seg_p_dis(L a,P x){
  if(dot(a.sc-a.fs,x-a.fs)<EPS)return abs(x-a.fs);
  if(dot(a.fs-a.sc,x-a.sc)<EPS)return abs(x-a.sc);
  return abs(cross(a.sc-a.fs,x-a.fs))/abs(a.sc-a.fs);
}

int ccw(P a,P b,P c){
  b -= a; c -= a;
  if(cross(b,c)>EPS)return 1;
  if(cross(b,c)<-EPS)return -1;
  if(dot(b,c)<-EPS)return 2;
  if(abs(b)<abs(c))return -2;
  return 0;
}

bool is_cp(L a,L b){
  if(ccw(a.fs,a.sc,b.fs)*ccw(a.fs,a.sc,b.sc)<=0)
    if(ccw(b.fs,b.sc,a.fs)*ccw(b.fs,b.sc,a.sc)<=0)return true;
  return false;
}

bool inter_pos(Poly p,P x){
  int s = p.size();
  if(!s)return false;
  rep(i,s)if(ccw(p[i],p[(i+1)%s],x)==-1)return false;
  return true;
}

bool cross_poly(Poly a, Poly b){
  int n1 = a.size(), n2 = b.size();
  rep(i,n2)if(inter_pos(a,b[i]))return true;
  rep(i,n1)if(inter_pos(b,a[i]))return true;
  rep(i,n1)rep(j,n2)
    if(is_cp(L(a[i],a[(i+1)%n1]),L(b[j],b[(j+1)%n2])))return true;
  return false;
}

void fdfs(int v){
  use[v] = true;
  for(int i=0;i<g[v].size();i++){
    if(!use[g[v][i]])fdfs(g[v][i]);
  }
  vs.pb(v);
}

void rdfs(int v,int k){
  use[v] = true;
  ord[v] = k;
  for(int i=0;i<rg[v].size();i++){
    if(!use[rg[v][i]])rdfs(rg[v][i],k);
  }
}

int scc(){
  memset(use,0,sizeof(use));
  vs.clear();
  for(int i=0;i<n;i++){
    if(!use[i])fdfs(i);
  }

  memset(use,0,sizeof(use));
  int k = 0;
  for(int i=vs.size()-1;i>=0;i--){
    if(!use[vs[i]])rdfs(vs[i],k++);
  }
  return k;
}

int main(){
  while(cin >> n >> d, n+d){
    rep(i,n){
      tri[i].resize(3);
      cin >> tri[i][0].real() >> tri[i][0].imag();
      cin >> tri[i][1].real() >> tri[i][1].imag();
      cin >> tri[i][2].real() >> tri[i][2].imag();
      tri[i] = convex_hull(tri[i]);

      rect[i].clear();
      rep(j,3){
	int n1 = (j+1)%3, n2 = (j+2)%3;
	if(EQ(abs(tri[i][n1]-tri[i][j]),abs(tri[i][n2]-tri[i][j]))){
	  rect[i].pb(tri[i][n1]); rect[i].pb(tri[i][n2]);
	  if(ccw(tri[i][n1],tri[i][n2],tri[i][j])!=1){
	    swap(rect[i][0],rect[i][1]); break;
	  }
	}
      }

      P ver = (D)d*unit(norm(rect[i][0]-rect[i][1]).sc);
      rect[i].pb(ver+rect[i][1]); rect[i].pb(ver+rect[i][0]);
      rect[i] = convex_hull(rect[i]);
    }

    rep(i,n){
      g[i].clear(); rg[i].clear();
      rep(j,n){
	if(i==j)continue;
	if(cross_poly(rect[i],tri[j])){g[i].pb(j); rg[j].pb(i);}
      }
    }

    int cmp = scc();

    int in[110] = {0};
    rep(i,n){
      rep(j,g[i].size())if(ord[g[i][j]] != ord[i])in[ord[g[i][j]]]++;
    }

    int res = 0;
    rep(i,cmp)if(!in[i])res++;
    cout << res << endl;
      
  }
}