#include <stdio.h>
#include <string.h>
#include <math.h>
#include <float.h>
 
#define FLG_CONTACT_LESS (0)
#define FLG_A_TOUCH_ON_B (1)
#define FLG_A_IS_IN_B    (2)
 
#define SIMULATION_SECONDS (20.0)
 
typedef struct {
    /* 球の物理的な値 */
    double r; /* 球の半径(m) */
    double m; /* 球の質量(kg) */
 
    /* 球の初期状態を表す変数 */
    double x0; /* 球のx軸上の初期位置 (m) */
    double y0; /* 球のy軸上の初期位置 (m) */
    double vx0; /* 球のx軸上の初期速度 (m/s) */
    double vy0; /* 球のy軸上の初期速度 (m/s) */
 
    /* 球の現在の状態を表す変数 */
    double x; /* 球の中心座標のx軸上の現在地 (m) */
    double y; /* 球の中心座標のy軸上の現在地 (m) */
    double vx; /* x軸方向の現在の速度 (m/s) */
    double vy; /* y軸方向の現在の速度 (m/s) */
 
} sphere_object_t;
 
/* 正負は数学の平面を表すxy軸と一致する。
   実験1:
   物体a: 半径 1cm, 100g, 初期位置右に 10cm , 上に 0cm の状態で、
          左に 1 cm/sec の速度で進む。
   物体b: 半径 1cm, 100g, 初期位置左に 10cm , 上に 0cm の状態で、
          右に 1 cm/sec の速度で進む。
   10 秒後に物体a,b が衝突する。
 
   実験2: 現在は適当。
   物体a  半径 2cm, 100g, 初期位置右に 10cm , 上に 0cm の状態で、
          左に 1 cm/sec の速度で進む。
   物体b: 半径 10cm, 100g, 初期位置左に 10cm , 上に 0cm の状態で、
          右に 1 cm/sec の速度で進む。
 */
 
int so_init(sphere_object_t *so,
            double r, double m,
            double x0, double y0,
            double vx0, double vy0)
{
    /* 球の物理的な値 */
    so->r = r; /* 球の半径(m) */
    so->m = m; /* 球の質量(kg) */
 
    /* 球の初期状態を表す変数 */
    so->x0 = x0; /* 球のx軸上の初期位置 (m) */
    so->y0 = y0; /* 球のy軸上の初期位置 (m) */
    so->vx0 = vx0; /* 球のx軸上の初期速度 (m/s) */
    so->vy0 = vy0; /* 球のy軸上の初期速度 (m/s) */
 
    /* 球の現在の状態を表す変数 */
    so->x = x0; /* 球の中心座標のx軸上の現在地 (m) */
    so->y = y0; /* 球の中心座標のy軸上の現在地 (m) */
    so->vx = vx0; /* x軸方向の現在の速度 (m/s) */
    so->vy = vy0; /* y軸方向の現在の速度 (m/s) */
 
    return 0;
}
 
int so_view_current(sphere_object_t *so)
{
    fprintf(stdout, "                       so = %p\n", so);
 
    /* 球の現在の状態 */
    fprintf(stdout, " so->(r, x, y, m, vx, vy) = (%.3f, %.3f, %.3f, %.3f, %.3f, %.3f)\n", so->r, so->x, so->y, so->m, so->vx, so->vy);
 
    return 0;
}
 
int so_view(sphere_object_t *so)
{
    fprintf(stdout, "                    so = %p\n", so);
    /* 球の物理的な値 */
    fprintf(stdout, "            so->(r, m) = (%.3f, %.3f)\n", so->r, so->m);
 
    /* 球の初期状態を表す変数 */
    fprintf(stdout, "so->(x0, y0, vx0, vy0) = (%.3f, %.3f, %.3f, %.3f)\n", so->x0, so->y0, so->vx0, so->vy0);
 
    /* 球の現在の状態を表す変数 */
    fprintf(stdout, "    so->(x, y, vx, vy) = (%.3f, %.3f, %.3f, %.3f)\n", so->x, so->y, so->vx, so->vy);
 
 
    return 0;
}
 
int so_copy(sphere_object_t *dst, sphere_object_t *src)
{
    memcpy(dst, src, sizeof(sphere_object_t));
 
    return 0;
}
 
int so_collision_detect(sphere_object_t *a, sphere_object_t *b)
{
    double distance_of_a_and_b, sum_of_a_r_and_b_r;
    double dx, dx2, dy, dy2;
    int hit_flg = FLG_CONTACT_LESS;
 
    sum_of_a_r_and_b_r = a->r + b->r;
 
    dx = a->x - b->x;
    dx2 = dx * dx;
    dy = a->y - b->y;
    dy2 = dy * dy;
    distance_of_a_and_b = sqrt(dx2 + dy2);
 
    fprintf(stdout, " sum_of_a_r_and_b_r = %f\n", sum_of_a_r_and_b_r);
    fprintf(stdout, "distance_of_a_and_b = %f\n", distance_of_a_and_b);
    if (fabs(distance_of_a_and_b - sum_of_a_r_and_b_r) < DBL_EPSILON){
        fprintf(stdout, "/* a touch on b */\n");
        hit_flg |= FLG_A_TOUCH_ON_B;
    }
    if (distance_of_a_and_b < sum_of_a_r_and_b_r) {
        fprintf(stdout, "/* a is in b */\n");
        hit_flg |= FLG_A_IS_IN_B;
    }
 
    return hit_flg;
}
 
/* sphere_object の何秒後の状態を知りたいか、sec引数にて指定する。*/
int so_update(sphere_object_t *dst, sphere_object_t *src, double sec)
{
    sphere_object_t tmp_, *tmp = &tmp_;
 
    so_copy(tmp, src);
    tmp->x += tmp->vx * sec;
    tmp->y += tmp->vy * sec;
    so_copy(dst, tmp);
 
    return 0;
}
 
int so_clash(double *v1_, double *v2_,
             double m1, double v1,
             double m2, double v2,
             double e)
{
    /*
    1.   m1 * v1 + m2 + v2 = m1 * v1' + m2 * v2'
    2.   v1' - v2' = -e * (v1 - v2)
                 e = - ((v1' - v2') / (v1 - v2))
    2-1. v2' = v1' + e * (v1 - v2)
    2-2. v1' = v2' - e * (v1 - v2)
    1. に 2-1. を代入して、
         m1 * v1 + m2 + v2 = m1 * v1' + m2 * (v1' + e * (v1 - v2))
                           = v1' * (m1 + m2) + e * m2 * (v1 - v2)
         v1' * (m1 + m2) = m1 * v1 + m2 * v2 - (e * m2 * (v1 - v2))
 
               m1 * v1 + m2 * v2 - m2 * (e * (v1 - v2))
         v1' = ----------------------------------------
               m1 + m2
 
    1. に 2-2. を代入して、
               m1 * v1 + m2 * v2 + m1 * (e * (v1 - v2))
         v2' = ----------------------------------------
               m1 + m2
 
    */
    *v1_ = (m1 * v1 + m2 * v2 - m2 * (e * (v1 - v2))) / (m1 + m2);
    *v2_ = (m1 * v1 + m2 * v2 + m1 * (e * (v1 - v2))) / (m1 + m2);
 
    return 0;
}
 
int so_change_mv(sphere_object_t *a, sphere_object_t *b, double e)
{
    /* v1 as a->vx, v2 as b->vx */
    so_clash(&a->vx, &b->vx, a->m, a->vx, b->m, b->vx, e);
    so_clash(&a->vy, &b->vy, a->m, a->vy, b->m, b->vy, e);
 
    return 0;
}
 
/* 完全弾性衝突の実験 */
void experiment_1(double e)
{
    sphere_object_t a_, b_, *a = &a_, *b = &b_;
    double sec, delta;
    int collision_flg;
 
    /*             r,    m,     x0,    y0,    vx0,   vy0 */
    so_init(a, 0.010, 0.10,  0.100, 0.000, -0.010, 0.000);
    so_init(b, 0.010, 0.10, -0.100, 0.000,  0.010, 0.000);
 
    delta = 1.0;
    for (sec=0.0;sec<SIMULATION_SECONDS + 1.0;sec+=delta) {
        fprintf(stdout, "現在は実験開始から %.1f 秒後の状態です。\n", sec);
        so_view_current(a);
        so_view_current(b);
        collision_flg = so_collision_detect(a, b);
        if (collision_flg) {
            fprintf(stdout, "a detect collision of b.(collision_flg=%d)\n", collision_flg);
            so_change_mv(a, b, e);
        }
        fprintf(stdout, "\n");
 
        so_update(a, a, delta);
        so_update(b, b, delta);
    }
}
 
int main(void)
{
    /* 完全弾性衝突の実験 */
    experiment_1(1.0);
 
    return 0;
}

I2luY2x1ZGUgPHN0ZGlvLmg+CiNpbmNsdWRlIDxzdHJpbmcuaD4KI2luY2x1ZGUgPG1hdGguaD4KI2luY2x1ZGUgPGZsb2F0Lmg+CgojZGVmaW5lIEZMR19DT05UQUNUX0xFU1MgKDApCiNkZWZpbmUgRkxHX0FfVE9VQ0hfT05fQiAoMSkKI2RlZmluZSBGTEdfQV9JU19JTl9CICAgICgyKQoKI2RlZmluZSBTSU1VTEFUSU9OX1NFQ09ORFMgKDIwLjApCgp0eXBlZGVmIHN0cnVjdCB7CiAgICAvKiDnkIPjga7niannkIbnmoTjgarlgKQgKi8KICAgIGRvdWJsZSByOyAvKiDnkIPjga7ljYrlvoQobSkgKi8KICAgIGRvdWJsZSBtOyAvKiDnkIPjga7os6rph48oa2cpICovCgogICAgLyog55CD44Gu5Yid5pyf54q25oWL44KS6KGo44GZ5aSJ5pWwICovCiAgICBkb3VibGUgeDA7IC8qIOeQg+OBrnjou7jkuIrjga7liJ3mnJ/kvY3nva4gKG0pICovCiAgICBkb3VibGUgeTA7IC8qIOeQg+OBrnnou7jkuIrjga7liJ3mnJ/kvY3nva4gKG0pICovCiAgICBkb3VibGUgdngwOyAvKiDnkIPjga546Lu45LiK44Gu5Yid5pyf6YCf5bqmIChtL3MpICovCiAgICBkb3VibGUgdnkwOyAvKiDnkIPjga556Lu45LiK44Gu5Yid5pyf6YCf5bqmIChtL3MpICovCgogICAgLyog55CD44Gu54++5Zyo44Gu54q25oWL44KS6KGo44GZ5aSJ5pWwICovCiAgICBkb3VibGUgeDsgLyog55CD44Gu5Lit5b+D5bqn5qiZ44GueOi7uOS4iuOBruePvuWcqOWcsCAobSkgKi8KICAgIGRvdWJsZSB5OyAvKiDnkIPjga7kuK3lv4PluqfmqJnjga556Lu45LiK44Gu54++5Zyo5ZywIChtKSAqLwogICAgZG91YmxlIHZ4OyAvKiB46Lu45pa55ZCR44Gu54++5Zyo44Gu6YCf5bqmIChtL3MpICovCiAgICBkb3VibGUgdnk7IC8qIHnou7jmlrnlkJHjga7nj77lnKjjga7pgJ/luqYgKG0vcykgKi8KCn0gc3BoZXJlX29iamVjdF90OwoKLyog5q2j6LKg44Gv5pWw5a2m44Gu5bmz6Z2i44KS6KGo44GZeHnou7jjgajkuIDoh7TjgZnjgovjgIIKICAg5a6f6aiTMToKICAg54mp5L2TYTog5Y2K5b6EIDFjbSwgMTAwZywg5Yid5pyf5L2N572u5Y+z44GrIDEwY20gLCDkuIrjgasgMGNtIOOBrueKtuaFi+OBp+OAgQogICAgICAgICAg5bem44GrIDEgY20vc2VjIOOBrumAn+W6puOBp+mAsuOCgOOAggogICDniankvZNiOiDljYrlvoQgMWNtLCAxMDBnLCDliJ3mnJ/kvY3nva7lt6bjgasgMTBjbSAsIOS4iuOBqyAwY20g44Gu54q25oWL44Gn44CBCiAgICAgICAgICDlj7PjgasgMSBjbS9zZWMg44Gu6YCf5bqm44Gn6YCy44KA44CCCiAgIDEwIOenkuW+jOOBq+eJqeS9k2EsYiDjgYzooZ3nqoHjgZnjgovjgIIKCiAgIOWun+mokzI6IOePvuWcqOOBr+mBqeW9k+OAggogICDniankvZNhICDljYrlvoQgMmNtLCAxMDBnLCDliJ3mnJ/kvY3nva7lj7PjgasgMTBjbSAsIOS4iuOBqyAwY20g44Gu54q25oWL44Gn44CBCiAgICAgICAgICDlt6bjgasgMSBjbS9zZWMg44Gu6YCf5bqm44Gn6YCy44KA44CCCiAgIOeJqeS9k2I6IOWNiuW+hCAxMGNtLCAxMDBnLCDliJ3mnJ/kvY3nva7lt6bjgasgMTBjbSAsIOS4iuOBqyAwY20g44Gu54q25oWL44Gn44CBCiAgICAgICAgICDlj7PjgasgMSBjbS9zZWMg44Gu6YCf5bqm44Gn6YCy44KA44CCCiAqLwoKaW50IHNvX2luaXQoc3BoZXJlX29iamVjdF90ICpzbywKICAgICAgICAgICAgZG91YmxlIHIsIGRvdWJsZSBtLAogICAgICAgICAgICBkb3VibGUgeDAsIGRvdWJsZSB5MCwKICAgICAgICAgICAgZG91YmxlIHZ4MCwgZG91YmxlIHZ5MCkKewogICAgLyog55CD44Gu54mp55CG55qE44Gq5YCkICovCiAgICBzby0+ciA9IHI7IC8qIOeQg+OBruWNiuW+hChtKSAqLwogICAgc28tPm0gPSBtOyAvKiDnkIPjga7os6rph48oa2cpICovCgogICAgLyog55CD44Gu5Yid5pyf54q25oWL44KS6KGo44GZ5aSJ5pWwICovCiAgICBzby0+eDAgPSB4MDsgLyog55CD44GueOi7uOS4iuOBruWIneacn+S9jee9riAobSkgKi8KICAgIHNvLT55MCA9IHkwOyAvKiDnkIPjga556Lu45LiK44Gu5Yid5pyf5L2N572uIChtKSAqLwogICAgc28tPnZ4MCA9IHZ4MDsgLyog55CD44GueOi7uOS4iuOBruWIneacn+mAn+W6piAobS9zKSAqLwogICAgc28tPnZ5MCA9IHZ5MDsgLyog55CD44Gueei7uOS4iuOBruWIneacn+mAn+W6piAobS9zKSAqLwoKICAgIC8qIOeQg+OBruePvuWcqOOBrueKtuaFi+OCkuihqOOBmeWkieaVsCAqLwogICAgc28tPnggPSB4MDsgLyog55CD44Gu5Lit5b+D5bqn5qiZ44GueOi7uOS4iuOBruePvuWcqOWcsCAobSkgKi8KICAgIHNvLT55ID0geTA7IC8qIOeQg+OBruS4reW/g+W6p+aomeOBrnnou7jkuIrjga7nj77lnKjlnLAgKG0pICovCiAgICBzby0+dnggPSB2eDA7IC8qIHjou7jmlrnlkJHjga7nj77lnKjjga7pgJ/luqYgKG0vcykgKi8KICAgIHNvLT52eSA9IHZ5MDsgLyogeei7uOaWueWQkeOBruePvuWcqOOBrumAn+W6piAobS9zKSAqLwoKICAgIHJldHVybiAwOwp9CgppbnQgc29fdmlld19jdXJyZW50KHNwaGVyZV9vYmplY3RfdCAqc28pCnsKICAgIGZwcmludGYoc3Rkb3V0LCAiICAgICAgICAgICAgICAgICAgICAgICBzbyA9ICVwXG4iLCBzbyk7CgogICAgLyog55CD44Gu54++5Zyo44Gu54q25oWLICovCiAgICBmcHJpbnRmKHN0ZG91dCwgIiBzby0+KHIsIHgsIHksIG0sIHZ4LCB2eSkgPSAoJS4zZiwgJS4zZiwgJS4zZiwgJS4zZiwgJS4zZiwgJS4zZilcbiIsIHNvLT5yLCBzby0+eCwgc28tPnksIHNvLT5tLCBzby0+dngsIHNvLT52eSk7CgogICAgcmV0dXJuIDA7Cn0KCmludCBzb192aWV3KHNwaGVyZV9vYmplY3RfdCAqc28pCnsKICAgIGZwcmludGYoc3Rkb3V0LCAiICAgICAgICAgICAgICAgICAgICBzbyA9ICVwXG4iLCBzbyk7CiAgICAvKiDnkIPjga7niannkIbnmoTjgarlgKQgKi8KICAgIGZwcmludGYoc3Rkb3V0LCAiICAgICAgICAgICAgc28tPihyLCBtKSA9ICglLjNmLCAlLjNmKVxuIiwgc28tPnIsIHNvLT5tKTsKCiAgICAvKiDnkIPjga7liJ3mnJ/nirbmhYvjgpLooajjgZnlpInmlbAgKi8KICAgIGZwcmludGYoc3Rkb3V0LCAic28tPih4MCwgeTAsIHZ4MCwgdnkwKSA9ICglLjNmLCAlLjNmLCAlLjNmLCAlLjNmKVxuIiwgc28tPngwLCBzby0+eTAsIHNvLT52eDAsIHNvLT52eTApOwoKICAgIC8qIOeQg+OBruePvuWcqOOBrueKtuaFi+OCkuihqOOBmeWkieaVsCAqLwogICAgZnByaW50ZihzdGRvdXQsICIgICAgc28tPih4LCB5LCB2eCwgdnkpID0gKCUuM2YsICUuM2YsICUuM2YsICUuM2YpXG4iLCBzby0+eCwgc28tPnksIHNvLT52eCwgc28tPnZ5KTsKCgogICAgcmV0dXJuIDA7Cn0KCmludCBzb19jb3B5KHNwaGVyZV9vYmplY3RfdCAqZHN0LCBzcGhlcmVfb2JqZWN0X3QgKnNyYykKewogICAgbWVtY3B5KGRzdCwgc3JjLCBzaXplb2Yoc3BoZXJlX29iamVjdF90KSk7CgogICAgcmV0dXJuIDA7Cn0KCmludCBzb19jb2xsaXNpb25fZGV0ZWN0KHNwaGVyZV9vYmplY3RfdCAqYSwgc3BoZXJlX29iamVjdF90ICpiKQp7CiAgICBkb3VibGUgZGlzdGFuY2Vfb2ZfYV9hbmRfYiwgc3VtX29mX2Ffcl9hbmRfYl9yOwogICAgZG91YmxlIGR4LCBkeDIsIGR5LCBkeTI7CiAgICBpbnQgaGl0X2ZsZyA9IEZMR19DT05UQUNUX0xFU1M7CgogICAgc3VtX29mX2Ffcl9hbmRfYl9yID0gYS0+ciArIGItPnI7CgogICAgZHggPSBhLT54IC0gYi0+eDsKICAgIGR4MiA9IGR4ICogZHg7CiAgICBkeSA9IGEtPnkgLSBiLT55OwogICAgZHkyID0gZHkgKiBkeTsKICAgIGRpc3RhbmNlX29mX2FfYW5kX2IgPSBzcXJ0KGR4MiArIGR5Mik7CgogICAgZnByaW50ZihzdGRvdXQsICIgc3VtX29mX2Ffcl9hbmRfYl9yID0gJWZcbiIsIHN1bV9vZl9hX3JfYW5kX2Jfcik7CiAgICBmcHJpbnRmKHN0ZG91dCwgImRpc3RhbmNlX29mX2FfYW5kX2IgPSAlZlxuIiwgZGlzdGFuY2Vfb2ZfYV9hbmRfYik7CiAgICBpZiAoZmFicyhkaXN0YW5jZV9vZl9hX2FuZF9iIC0gc3VtX29mX2Ffcl9hbmRfYl9yKSA8IERCTF9FUFNJTE9OKXsKICAgICAgICBmcHJpbnRmKHN0ZG91dCwgIi8qIGEgdG91Y2ggb24gYiAqL1xuIik7CiAgICAgICAgaGl0X2ZsZyB8PSBGTEdfQV9UT1VDSF9PTl9COwogICAgfQogICAgaWYgKGRpc3RhbmNlX29mX2FfYW5kX2IgPCBzdW1fb2ZfYV9yX2FuZF9iX3IpIHsKICAgICAgICBmcHJpbnRmKHN0ZG91dCwgIi8qIGEgaXMgaW4gYiAqL1xuIik7CiAgICAgICAgaGl0X2ZsZyB8PSBGTEdfQV9JU19JTl9COwogICAgfQoKICAgIHJldHVybiBoaXRfZmxnOwp9CgovKiBzcGhlcmVfb2JqZWN0IOOBruS9leenkuW+jOOBrueKtuaFi+OCkuefpeOCiuOBn+OBhOOBi+OAgXNlY+W8leaVsOOBq+OBpuaMh+WumuOBmeOCi+OAgiovCmludCBzb191cGRhdGUoc3BoZXJlX29iamVjdF90ICpkc3QsIHNwaGVyZV9vYmplY3RfdCAqc3JjLCBkb3VibGUgc2VjKQp7CiAgICBzcGhlcmVfb2JqZWN0X3QgdG1wXywgKnRtcCA9ICZ0bXBfOwoKICAgIHNvX2NvcHkodG1wLCBzcmMpOwogICAgdG1wLT54ICs9IHRtcC0+dnggKiBzZWM7CiAgICB0bXAtPnkgKz0gdG1wLT52eSAqIHNlYzsKICAgIHNvX2NvcHkoZHN0LCB0bXApOwoKICAgIHJldHVybiAwOwp9CgppbnQgc29fY2xhc2goZG91YmxlICp2MV8sIGRvdWJsZSAqdjJfLAogICAgICAgICAgICAgZG91YmxlIG0xLCBkb3VibGUgdjEsCiAgICAgICAgICAgICBkb3VibGUgbTIsIGRvdWJsZSB2MiwKICAgICAgICAgICAgIGRvdWJsZSBlKQp7CiAgICAvKgogICAgMS4gICBtMSAqIHYxICsgbTIgKyB2MiA9IG0xICogdjEnICsgbTIgKiB2MicKICAgIDIuICAgdjEnIC0gdjInID0gLWUgKiAodjEgLSB2MikKICAgICAgICAgICAgICAgICBlID0gLSAoKHYxJyAtIHYyJykgLyAodjEgLSB2MikpCiAgICAyLTEuIHYyJyA9IHYxJyArIGUgKiAodjEgLSB2MikKICAgIDItMi4gdjEnID0gdjInIC0gZSAqICh2MSAtIHYyKQogICAgMS4g44GrIDItMS4g44KS5Luj5YWl44GX44Gm44CBCiAgICAgICAgIG0xICogdjEgKyBtMiArIHYyID0gbTEgKiB2MScgKyBtMiAqICh2MScgKyBlICogKHYxIC0gdjIpKQogICAgICAgICAgICAgICAgICAgICAgICAgICA9IHYxJyAqIChtMSArIG0yKSArIGUgKiBtMiAqICh2MSAtIHYyKQogICAgICAgICB2MScgKiAobTEgKyBtMikgPSBtMSAqIHYxICsgbTIgKiB2MiAtIChlICogbTIgKiAodjEgLSB2MikpCgogICAgICAgICAgICAgICBtMSAqIHYxICsgbTIgKiB2MiAtIG0yICogKGUgKiAodjEgLSB2MikpCiAgICAgICAgIHYxJyA9IC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KICAgICAgICAgICAgICAgbTEgKyBtMgoKICAgIDEuIOOBqyAyLTIuIOOCkuS7o+WFpeOBl+OBpuOAgQogICAgICAgICAgICAgICBtMSAqIHYxICsgbTIgKiB2MiArIG0xICogKGUgKiAodjEgLSB2MikpCiAgICAgICAgIHYyJyA9IC0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0tLS0KICAgICAgICAgICAgICAgbTEgKyBtMgoKICAgICovCiAgICAqdjFfID0gKG0xICogdjEgKyBtMiAqIHYyIC0gbTIgKiAoZSAqICh2MSAtIHYyKSkpIC8gKG0xICsgbTIpOwogICAgKnYyXyA9IChtMSAqIHYxICsgbTIgKiB2MiArIG0xICogKGUgKiAodjEgLSB2MikpKSAvIChtMSArIG0yKTsKCiAgICByZXR1cm4gMDsKfQoKaW50IHNvX2NoYW5nZV9tdihzcGhlcmVfb2JqZWN0X3QgKmEsIHNwaGVyZV9vYmplY3RfdCAqYiwgZG91YmxlIGUpCnsKICAgIC8qIHYxIGFzIGEtPnZ4LCB2MiBhcyBiLT52eCAqLwogICAgc29fY2xhc2goJmEtPnZ4LCAmYi0+dngsIGEtPm0sIGEtPnZ4LCBiLT5tLCBiLT52eCwgZSk7CiAgICBzb19jbGFzaCgmYS0+dnksICZiLT52eSwgYS0+bSwgYS0+dnksIGItPm0sIGItPnZ5LCBlKTsKCiAgICByZXR1cm4gMDsKfQoKLyog5a6M5YWo5by+5oCn6KGd56qB44Gu5a6f6aiTICovCnZvaWQgZXhwZXJpbWVudF8xKGRvdWJsZSBlKQp7CiAgICBzcGhlcmVfb2JqZWN0X3QgYV8sIGJfLCAqYSA9ICZhXywgKmIgPSAmYl87CiAgICBkb3VibGUgc2VjLCBkZWx0YTsKICAgIGludCBjb2xsaXNpb25fZmxnOwoKICAgIC8qICAgICAgICAgICAgIHIsICAgIG0sICAgICB4MCwgICAgeTAsICAgIHZ4MCwgICB2eTAgKi8KICAgIHNvX2luaXQoYSwgMC4wMTAsIDAuMTAsICAwLjEwMCwgMC4wMDAsIC0wLjAxMCwgMC4wMDApOwogICAgc29faW5pdChiLCAwLjAxMCwgMC4xMCwgLTAuMTAwLCAwLjAwMCwgIDAuMDEwLCAwLjAwMCk7CgogICAgZGVsdGEgPSAxLjA7CiAgICBmb3IgKHNlYz0wLjA7c2VjPFNJTVVMQVRJT05fU0VDT05EUyArIDEuMDtzZWMrPWRlbHRhKSB7CiAgICAgICAgZnByaW50ZihzdGRvdXQsICLnj77lnKjjga/lrp/pqJPplovlp4vjgYvjgokgJS4xZiDnp5Llvozjga7nirbmhYvjgafjgZnjgIJcbiIsIHNlYyk7CiAgICAgICAgc29fdmlld19jdXJyZW50KGEpOwogICAgICAgIHNvX3ZpZXdfY3VycmVudChiKTsKICAgICAgICBjb2xsaXNpb25fZmxnID0gc29fY29sbGlzaW9uX2RldGVjdChhLCBiKTsKICAgICAgICBpZiAoY29sbGlzaW9uX2ZsZykgewogICAgICAgICAgICBmcHJpbnRmKHN0ZG91dCwgImEgZGV0ZWN0IGNvbGxpc2lvbiBvZiBiLihjb2xsaXNpb25fZmxnPSVkKVxuIiwgY29sbGlzaW9uX2ZsZyk7CiAgICAgICAgICAgIHNvX2NoYW5nZV9tdihhLCBiLCBlKTsKICAgICAgICB9CiAgICAgICAgZnByaW50ZihzdGRvdXQsICJcbiIpOwoKICAgICAgICBzb191cGRhdGUoYSwgYSwgZGVsdGEpOwogICAgICAgIHNvX3VwZGF0ZShiLCBiLCBkZWx0YSk7CiAgICB9Cn0KCmludCBtYWluKHZvaWQpCnsKICAgIC8qIOWujOWFqOW8vuaAp+ihneeqgeOBruWun+mokyAqLwogICAgZXhwZXJpbWVudF8xKDEuMCk7CgogICAgcmV0dXJuIDA7Cn0=