import java.awt.*;
import java.applet.*;
import java.awt.image.*;
import java.util.Vector;
import java.awt.event.*;


public class CircleInv extends Applet {
  CirclePanel circlepanel;

  public void init() {
    setLayout(new BorderLayout());
    this.circlepanel=new CirclePanel();
    add("Center", this.circlepanel);
  }
}


class CircleData {
  double xc[],yc[],rad[],xtan[],ytan[];
  int num;
  CircleData() {
    num=4;
    xc=new double[4];
    yc=new double[4];
    rad=new double[4];
    xtan=new double[4];
    ytan=new double[4];
      xc[0]=1.25;yc[0]=.75;rad[0]=.25;xtan[0]=1.5;ytan[0]=.75;
      xc[1]=1.75;yc[1]=.75;rad[1]=.25;xtan[1]=1.75;ytan[1]=1;
      xc[2]=1.75;yc[2]=1.25;rad[2]=.25;xtan[2]=1.5;ytan[2]=1.25;
      xc[3]=1.25;yc[3]=1.25;rad[3]=.25;xtan[3]=1.25;ytan[3]=1;     
  }

  CircleData(int num){
    this.num=num;
    xc=new double[num];
    yc=new double[num];
    rad=new double[num];
    xtan=new double[num];
    ytan=new double[num];
    for (int i=0; i<num; i++){
	xc[i]=1.5+.5*Math.cos(2*Math.PI*i/num);
 	yc[i]=1+.5*Math.sin(2*Math.PI*i/num);
	rad[i]=.5*Math.sin(Math.PI/num);
	xtan[i]=1.5+.5*Math.cos(Math.PI/num)*Math.cos(Math.PI*(2*i+1)/num);
 	ytan[i]=1+.5*Math.cos(Math.PI/num)*Math.sin(Math.PI*(2*i+1)/num);
    }
  }
  
  void shiftCent(int cent,double a, double b){
    double x1,y1,r1,x2,y2,r2,d,xp1,xp2,yp1,yp2,t1,t2;
    x1=xc[(cent+num-1)%num];
    y1=yc[(cent+num-1)%num];
    r1=rad[(cent+num-1)%num];
    x2=xc[(cent+1)%num];
    y2=yc[(cent+1)%num];
    r2=rad[(cent+1)%num];
    d=Math.sqrt((x2-x1)*(x2-x1)+(y2-y1)*(y2-y1));
    yp1=(x2-x1)/d;
    yp2=(y2-y1)/d;
    xp1=yp2;
    xp2=-yp1;
    t1=-(b*yp1-y1*yp1-a*yp2+x1*yp2)/(-xp2*yp1+xp1*yp2);
    t2=(d*d*d-d*(r1-r2)*(r1-r2)+Math.sqrt(d*d-(r1-r2)*(r1-r2))*(r1-r2)*
      Math.sqrt(d*d-(r1-r2)*(r1-r2)+4*t1*t1))/(2*(d*d-(r1-r2)*(r1-r2)));
    rad[cent]=Math.sqrt(t2*t2+t1*t1)-r1;
    xc[cent]=t1*xp1+t2*yp1+x1;
    yc[cent]=t1*xp2+t2*yp2+y1;
    xtan[cent]=rad[cent]/(rad[cent]+rad[(cent+1)%num])*(xc[(cent+1)%num]-xc[cent])+xc[cent];
    ytan[cent]=rad[cent]/(rad[cent]+rad[(cent+1)%num])*(yc[(cent+1)%num]-yc[cent])+yc[cent];
    xtan[(cent+num-1)%num]=rad[(cent+num-1)%num]/(rad[(cent+num-1)%num]+rad[cent])*(xc[cent]-xc[(cent+num-1)%num])+xc[(cent+num-1)%num];
    ytan[(cent+num-1)%num]=rad[(cent+num-1)%num]/(rad[(cent+num-1)%num]+rad[cent])*(yc[cent]-yc[(cent+num-1)%num])+yc[(cent+num-1)%num];
 }

  CCircle invCircle(int num, CCircle C){
    double d=Math.sqrt((xc[num]-C.x)*(xc[num]-C.x)+(yc[num]-C.y)*(yc[num]-C.y));
    double xd=(C.x-xc[num])/d;
    double yd=(C.y-yc[num])/d;
    double r2=(rad[num]*.5*(rad[num]/(d-C.rad)+rad[num]/(d+C.rad)));
    double r3=(rad[num]*.5*(rad[num]/(d-C.rad)-rad[num]/(d+C.rad)));
    return new CCircle(xc[num]+r2*xd,yc[num]+r2*yd,r3,C.c,num);
  }

  TwoPoint invTwoPoint(int num, TwoPoint T){
    double d1=Math.sqrt((xc[num]-T.x1)*(xc[num]-T.x1)+(yc[num]-T.y1)*(yc[num]-T.y1));
    double x1=xc[num]+(T.x1-xc[num])*rad[num]*rad[num]/d1/d1;
    double y1=yc[num]+(T.y1-yc[num])*rad[num]*rad[num]/d1/d1;
    double d2=Math.sqrt((xc[num]-T.x2)*(xc[num]-T.x2)+(yc[num]-T.y2)*(yc[num]-T.y2));
    double x2=xc[num]+(T.x2-xc[num])*rad[num]*rad[num]/d2/d2;
    double y2=yc[num]+(T.y2-yc[num])*rad[num]*rad[num]/d2/d2;
    return new TwoPoint(x1,y1,x2,y2,num);
  }

  void zoom(double d){
    for (int i=0;i<num;i++){
      xc[i]=1.5+(xc[i]-1.5)*d;
      yc[i]=1+(yc[i]-1)*d;
      rad[i]=rad[i]*d;
      xtan[i]=1.5+(xtan[i]-1.5)*d;
      ytan[i]=1+(ytan[i]-1)*d;
    }
  }

  void zoom(double d,double x, double y){
    for (int i=0;i<num;i++){
      xc[i]=1.5+(xc[i]-x)*d;
      yc[i]=1+(yc[i]-y)*d;
      rad[i]=rad[i]*d;
      xtan[i]=1.5+(xtan[i]-x)*d;
      ytan[i]=1+(ytan[i]-y)*d;
    }
  }

  void move(double x, double y){
    for (int i=0;i<num;i++){
      xc[i]=1.5+xc[i]-x;
      yc[i]=1+yc[i]-y;
      xtan[i]=1.5+xtan[i]-x;
      ytan[i]=1+ytan[i]-y;
    }
  }

}

class CCircle{
  double x, y, rad;
  int c,loc;
  CCircle (double x, double y, double r, int c, int loc){
    this.x=x;
    this.y=y;
    this.rad=r;
    this.c=c;
    this.loc=loc;
  }
}

class TwoPoint{
  double x1, y1, x2, y2;
  int loc;
  TwoPoint (double x1, double y1, double x2, double y2, int loc){
    this.x1=x1;
    this.y1=y1;
    this.x2=x2;
    this.y2=y2;
    this.loc=loc;
  }
}

class CCircleQueue{
  double x[], y[], r[];
  int c[], l[];
  int start, insert;
  int num;

  CCircleQueue(int num){
    this.num=num;
    x=new double[num];
    y=new double[num];
    r=new double[num];
    c=new int [num];
    l=new int [num];
    start=0;
    insert=0;
  }

  boolean isFull(){
     return ((insert+1)%num==start);
  } 

  boolean isEmpty(){
      return (start==insert);
  }
  boolean add(double xn, double yn, double rn, int cn, int ln){
    if (isFull())
      return false;
    x[insert]=xn;
    y[insert]=yn;
    r[insert]=rn;
    c[insert]=cn;
    l[insert]=ln;
    insert=((insert+1)%num);
    return true;
  }

  boolean add(CCircle C){
    if (isFull())
      return false;
    x[insert]=C.x;
    y[insert]=C.y;
    r[insert]=C.rad;
    c[insert]=C.c;
    l[insert]=C.loc;
    insert=((insert+1)%num);
    return true;
  }    

  CCircle get(){
    CCircle cc=new CCircle(x[start],y[start],r[start],c[start],l[start]);
    start=((start+1)%num);
    return cc;
  }

  void clear(){
    start=0;
    insert=0;
  }
}


class CirclePanel extends Panel implements MouseListener, MouseMotionListener {
  int mult;
  CircleData circledata;
  boolean movingCenter=false;
  boolean zoomin=false,moveit=false;
  int cent;
  int drawmode=-1,prevmode;
  int blankx, blanky, blankr;
  CCircleQueue ccq;
  int num;
  int detail=40000;
  boolean filled=true;
  CirclePanel() {
    ccq=new CCircleQueue(1000);
    circledata=new CircleData(4); num=4;
    addMouseMotionListener(this);
    addMouseListener(this);
  }

  public void paint(Graphics g) {
    update(g);
  }

  void setColor(Graphics g, int i) {
      switch (i%9) {
      case (0):g.setColor(Color.red); break;
      case (1):g.setColor(Color.blue); break;
      case (2):g.setColor(Color.green); break;
      case (3):g.setColor(Color.yellow); break;
      case (4):g.setColor(Color.cyan); break;
      case (5):g.setColor(Color.magenta); break;
      case (6):g.setColor(Color.orange); break;
      case (7):g.setColor(Color.pink); break;
      case (8):g.setColor(Color.gray); break;
      }       
  }

  public void update(Graphics g) {
    mult=Math.min(getSize().width/3,getSize().height/2);
    switch (drawmode){
    case (-1):g.setColor(Color.white);g.fillRect(0,0,getSize().width-1,getSize().height-1);drawmode=0;
    case (0): paintInside(g); break;
    case (1): paintMoving(g); break;
    case (2): paintCurve2(g); break;
   }
  }

  void paintMenu(Graphics g){
    g.setFont( new Font("TimesRoman", Font.PLAIN, 10));
    g.setColor(Color.black);
    g.fillRect(1,1,60,12);
    g.setColor(Color.white);
    if (drawmode==0)
	g.drawString("View Curve",4,10);
    if (drawmode==2)
      g.drawString("View Circles",4,10);


    g.setColor(Color.black);
    g.fillRect(1,14,60,12);    
    g.setColor(Color.white);
    if (filled)
	g.drawString("Not Filled",4,23);
    else
	g.drawString("Filled",4,23);

    g.setColor(Color.black);
    g.fillRect(1,getSize().height-15,60,12);    
    g.setColor(Color.white);
    g.drawString("Reset",4,getSize().height-6);

    g.setColor(Color.black);
    g.fillRect(1,27,60,12);    
    g.setColor(Color.white);
    g.drawString("Zoom Out",4,36);

    g.setColor(Color.black);
    g.fillRect(1,40,60,12);    
    if (zoomin){
	g.setColor(Color.cyan);
	g.fillRect(2,41,58,10);    
	g.setColor(Color.black);
	g.drawString("Zoom In",4,49);
    } else {
	g.setColor(Color.white);
	g.drawString("Zoom In",4,49);
    }

    g.setColor(Color.black);
    g.fillRect(1,53,60,12);
    if (moveit){
	g.setColor(Color.cyan);
	g.fillRect(2,54,58,10);    
	g.setColor(Color.black);
	g.drawString("Move",4,62);
    } else {
	g.setColor(Color.white);
	g.drawString("Move",4,62);
    }
    
    g.setColor(Color.black);
    g.fillRect(1,79,60,12);    
    g.setColor(Color.white);
    g.drawString("Add Circle",4,88);
    
    if (num>3){
      g.setColor(Color.black);
      g.fillRect(1,92,60,12);    
      g.setColor(Color.white);
      g.drawString("Remove Circle",4,101);
    }

    if (drawmode==2){
      g.setColor(Color.black);
      g.fillRect(1,118,60,12);    
      g.setColor(Color.white);
      g.drawString("More Detail",4,127);

      g.setColor(Color.black);
      g.fillRect(1,131,60,12);    
      g.setColor(Color.white);
      g.drawString("Less Detail",4,140);
    }
  }      
    
  boolean menuClicked(int x, int y){
    if ((x>=1)&&(x<1+60)&&(y>=1)&&(y<1+12)){
	if (drawmode==0)
	    drawmode=2;
	else
	    drawmode=0;
	zoomin=false;
	moveit=false;
	repaint();
	return true;
    }
    if ((x>=1)&&(x<1+60)&&(y>=14)&&(y<14+12)){
	filled=!filled;
	zoomin=false;
	moveit=false;
	repaint();
	return true;
    }      
    if ((x>=1)&&(x<1+60)&&(y>=getSize().height-15)&&(y<getSize().height-15+12)){
      circledata=new CircleData(num);
      detail=40000;
      zoomin=false;
      moveit=false;
      repaint();
      return true;
    }
    if ((x>=1)&&(x<1+60)&&(y>=27)&&(y<27+12)){
      circledata.zoom(.5);
      zoomin=false;
      moveit=false;
      repaint();
      return true;
    }
    if ((x>=1)&&(x<1+60)&&(y>=40)&&(y<40+12)){
      zoomin=(!zoomin);
      moveit=false;
      repaint();
      return true;
    }

    if ((x>=1)&&(x<1+60)&&(y>=53)&&(y<53+12)){
      moveit=(!moveit);
      zoomin=false;
      repaint();
      return true;
    }

    if ((x>=1)&&(x<1+60)&&(y>=79)&&(y<79+12)){
      num++;
      circledata=new CircleData(num);
      zoomin=false;
      moveit=false;
      repaint();
      return true;
    }

    if (num>3){
      if ((x>=1)&&(x<1+60)&&(y>=92)&&(y<92+12)){
	num--;
	circledata=new CircleData(num);
	zoomin=false;
	moveit=false;
	repaint();
	return true;
      }
    }

    if (drawmode==2){
      if ((x>=1)&&(x<1+60)&&(y>=118)&&(y<118+12)){
	detail=detail*3;
	zoomin=false;
	moveit=false;
	repaint();
	return true;
      }

      if ((x>=1)&&(x<1+60)&&(y>=131)&&(y<131+12)){
	detail=detail/3;
	zoomin=false;
	moveit=false;
	repaint();
	return true;
      }
    }

    return false;
  }

  void paintMoving(Graphics g){
     g.setColor(Color.white);
     g.fillOval(blankx-blankr-2,blanky-blankr-2,2*blankr+4,2*blankr+4);
     for (int i=0; i<num; i++){
	 setColor(g,i);
	 if (filled)
	     g.fillOval(conv(circledata.xc[i]-circledata.rad[i]),conv(circledata.yc[i]-circledata.rad[i]),
		    conv(2*circledata.rad[i]),conv(2*circledata.rad[i]));
	 else
	     g.drawOval(conv(circledata.xc[i]-circledata.rad[i]),conv(circledata.yc[i]-circledata.rad[i]),
		    conv(2*circledata.rad[i]),conv(2*circledata.rad[i]));
     }
    g.setColor(Color.black);
    for (int i=0;i<num; i++){
    	plotPoint(g,conv(circledata.xc[i]),conv(circledata.yc[i]));
    	//plotPoint(g,conv(circledata.xtan[i]),conv(circledata.ytan[i]));
    } 
  }
  
  void paintInside(Graphics g){
    CCircle cc;
    mult=Math.min(getSize().width/3,getSize().height/2);
    g.setColor(Color.white);
    g.fillRect(0,0,getSize().width-1,getSize().height-1);
    ccq.clear();
    for (int i=0;i<num; i++)   
      ccq.add(circledata.xc[i],circledata.yc[i],circledata.rad[i],i,i);
    int count=0;
    while ((!(ccq.isEmpty()))&&(count<100000)){
      count++;
      cc=ccq.get();
      setColor(g,cc.c);
      if (filled)
	  g.fillOval(conv(cc.x-cc.rad),conv(cc.y-cc.rad),conv(2*cc.rad),conv(2*cc.rad));
      else
	  g.drawOval(conv(cc.x-cc.rad),conv(cc.y-cc.rad),conv(2*cc.rad),conv(2*cc.rad));	  
      if (cc.rad>.005){
	  for (int i=1;i<num;i++)
	      ccq.add(circledata.invCircle((cc.loc+i)%num,cc));
      }         
    } 
    g.setColor(Color.black);
    for (int i=0;i<num; i++){
    	plotPoint(g,conv(circledata.xc[i]),conv(circledata.yc[i]));
    	//plotPoint(g,conv(circledata.xtan[i]),conv(circledata.ytan[i]));
    } 
    paintMenu(g);
  }

  void paintCurve2(Graphics g){
    g.setColor(Color.white);
    g.fillRect(0,0,getSize().width-1,getSize().height-1);
    for (int i=0; i<num; i++){
      setColor(g,i);
      if (filled)
	  g.fillOval(conv(circledata.xc[i]-circledata.rad[i]),conv(circledata.yc[i]-circledata.rad[i]),
		     conv(2*circledata.rad[i]),conv(2*circledata.rad[i]));
      else
		  g.drawOval(conv(circledata.xc[i]-circledata.rad[i]),conv(circledata.yc[i]-circledata.rad[i]),
		     conv(2*circledata.rad[i]),conv(2*circledata.rad[i]));  
    }
    g.setColor(Color.black);
    for (int i=0;i<num; i++)
      paintRepeat(g,new TwoPoint(circledata.xtan[(i+num-1)%num],circledata.ytan[(i+num-1)%num],
				 circledata.xtan[i],circledata.ytan[i],i),
		  (int)(Math.floor(Math.log(detail/num)/Math.log(num-1))));
    for (int i=0;i<num; i++){
    	plotPoint(g,conv(circledata.xc[i]),conv(circledata.yc[i]));
    	//plotPoint(g,conv(circledata.xtan[i]),conv(circledata.ytan[i]));
    } 
    paintMenu(g);
  }

  void paintRepeat(Graphics g,TwoPoint T, int n){
     if (n>0) {
       for (int i=1;i<num; i++)
	 paintRepeat(g,circledata.invTwoPoint((T.loc+i)%num,T),n-1);
     } else
       g.drawLine(conv(T.x1),conv(T.y1),conv(T.x2),conv(T.y2));
  }

  int conv(double x){
    return (int)(x*mult);
  }
  
  double inv(int x){
    return (double)(x)/mult;
  }

  void plotPoint(Graphics g, int x, int y){
    g.drawRect(x-2,y-2,5,5);
  }
  
  boolean isInPoint(int x1, int y1, int x2, int y2){
    return ((Math.abs(x1-x2)<=3)&&(Math.abs(y1-y2)<=3));
  }
  
  public void mousePressed(MouseEvent e) {
    e.consume();
    for (int i=0; i<num; i++){
      if (isInPoint(e.getX(),e.getY(),conv(circledata.xc[i]),conv(circledata.yc[i]))){
        movingCenter=true;
        cent=i;
      }
    }
  }
 
  public void mouseDragged(MouseEvent e) {
    e.consume();
    if (movingCenter){
      blankx=conv(circledata.xc[cent]);
      blanky=conv(circledata.yc[cent]);
      blankr=conv(circledata.rad[cent]);
      circledata.shiftCent(cent,inv(e.getX()),inv(e.getY()));
      if (drawmode!=1)
	prevmode=drawmode;
      drawmode=1;
      repaint();
    }
  }

  public void mouseMoved(MouseEvent e) {
  }

  public void mouseReleased(MouseEvent e) {
    e.consume();
    if (movingCenter){
      circledata.shiftCent(cent,inv(e.getX()),inv(e.getY()));
      movingCenter=false;
      drawmode=prevmode;
      repaint();
    }
  }

  public void mouseEntered(MouseEvent e) {
  }

  public void mouseExited(MouseEvent e) {
  }

  public void mouseClicked(MouseEvent e) {
    e.consume();
    if (! menuClicked(e.getX(),e.getY())) {
	if (zoomin){
	    circledata.zoom(2,inv(e.getX()),inv(e.getY()));
	    repaint();
	}
	if (moveit){
	    circledata.move(inv(e.getX()),inv(e.getY()));
	    repaint();
	}
    }  
  }
}