#include <stdio.h>
#include <projects.h>

#define WKBUFSIZE 1024

static char* JGD2000_ZAHYOKEI[]
= {"+init=epsg:2443", "+init=epsg:2444", "+init=epsg:2445", "+init=epsg:2446",
   "+init=epsg:2447", "+init=epsg:2448", "+init=epsg:2449", "+init=epsg:2450",
   "+init=epsg:2451", "+init=epsg:2452", "+init=epsg:2453", "+init=epsg:2454",
   "+init=epsg:2455", "+init=epsg:2456", "+init=epsg:2457", "+init=epsg:2458",
   "+init=epsg:2459", "+init=epsg:2460", "+init=epsg:2461"};

int main()
{
  char* str_param[1];
  char  wkbuf[WKBUFSIZE];
  float x,y;
  int a;

  PJ* proj_param;
  XY idata;
  LP odata;

  while(fgets(wkbuf, WKBUFSIZE, stdin) != NULL){
    /* パラメータが Y,X の順なのに注意！ */
    if(sscanf(wkbuf, "%f %f %d", &y, &x, &a) != 3){
        fprintf(stderr, "format error: %s\n", wkbuf);
        continue;
    }
    if(a < 1 || a > (sizeof JGD2000_ZAHYOKEI / sizeof JGD2000_ZAHYOKEI[0])){
      fprintf(stderr, "illegal number: %d\n", a);
      continue;
    }
    str_param[0] = JGD2000_ZAHYOKEI[a-1];
    idata.u = x;
    idata.v = y;
    proj_param = pj_init(sizeof(str_param)/sizeof(char *), str_param);
    odata = pj_inv(idata, proj_param);
    pj_free(proj_param);
    odata.u /= DEG_TO_RAD;
    odata.v /= DEG_TO_RAD;
    printf("%.6f %.6f\n", odata.u, odata.v);
  }
  return 0;
}