結果

問題 No.3344 Common Tangent Line
コンテスト
ユーザー tobisatis
提出日時 2025-11-14 00:38:10
言語 C#
(.NET 8.0.404)
結果
WA  
実行時間 -
コード長 2,698 bytes
コンパイル時間 10,701 ms
コンパイル使用メモリ 171,508 KB
実行使用メモリ 197,056 KB
最終ジャッジ日時 2025-11-14 00:38:58
合計ジャッジ時間 42,988 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample RE * 1
other WA * 40
権限があれば一括ダウンロードができます
コンパイルメッセージ
  復元対象のプロジェクトを決定しています...
  /home/judge/data/code/main.csproj を復元しました (112 ミリ秒)。
  main -> /home/judge/data/code/bin/Release/net8.0/main.dll
  main -> /home/judge/data/code/bin/Release/net8.0/publish/

ソースコード

diff # 

#nullable enable

using System.Numerics;

#region
var (_input, _iter) = (Array.Empty<string>(), 0);
T I<T>() where T : IParsable<T>
{
    while (_iter >= _input.Length) (_input, _iter) = (Console.ReadLine()!.Split(' '), 0);
    return T.Parse(_input[_iter++], null);
}
#endregion

static double V(Complex p1, Complex p2, double scale, Complex rotation, double sx, double sy)
{
    p1 = p1 * scale * rotation;
    p2 = p2 * scale * rotation;
    var x1 = p1.Real + sx;
    var y1 = p1.Imaginary + sy;
    var x2 = p2.Real + sx;
    var y2 = p2.Imaginary + sy;
    var m = (y2 - y1) / (x2 - x1);
    var c = y1 - x1 * m;
    var s = Math.Max(1, Math.Max(Math.Abs(m), Math.Abs(c)));
    return Math.Abs(m + c + 1) / s;
}

static double F((double, double, double) c1, (double, double, double) c2)
{
    if (c1.Item3 > c2.Item3) (c1, c2) = (c2, c1);
    var (x1, y1, r1) = c1;
    var (x2, y2, r2) = c2;
    x2 -= x1;
    y2 -= y1;
    var arg2 = new Complex(x2, y2);
    arg2 /= arg2.Magnitude;
    var center2 = new Complex(x2, y2) / arg2;
    (x2, y2) = (center2.Real, center2.Imaginary);
    x2 /= r1;
    y2 /= r1;
    r2 /= r1;
    var h = Math.PI / 2;
    var argI = new Complex(-arg2.Real, arg2.Imaginary);

    var res = 0.0;
    {
        var p1 = new Complex();
        var p2 = p1;
        var pass = 0.0;
        var fail = h;
        while (Math.Abs(pass - fail) > 0.00000001)
        {
            var mid = (pass + fail) / 2;
            p1 = new Complex(Math.Cos(mid), Math.Sin(mid));
            p2 = new Complex(x2, y2) - p1 * r2;
            var v = p1 - p2;
            var r = (v / p1).Phase;
            if (r > h) pass = mid;
            else fail = mid;
        }
        res += V(p1, p2, r1, arg2, x1, y1);
        res += V(p1, p2, r1, argI, x1, y1);
    }
    {
        var p3 = new Complex();
        var p4 = p3;
        var pass = h;
        var fail = Math.PI;
        while (Math.Abs(pass - fail) > 0.00000001)
        {
            var mid = (pass + fail) / 2;
            p3 = new Complex(Math.Cos(mid), Math.Sin(mid));
            p4 = new Complex(x2, y2) + p3 * r2;
            var v = p3 - p4;
            var r = (v / p3).Phase;
            if (r > h) pass = mid;
            else fail = mid;
        }
        res += V(p3, p4, r1, arg2, x1, y1);
        res += V(p3, p4, r1, argI, x1, y1);
    }
    return res;
}

var t = I<int>();
var ans = new List<string>();
for (var i = 0; i < t; i++)
{
    var x1 = I<int>();
    var y1 = I<int>();
    var r1 = I<int>();
    var x2 = I<int>();
    var y2 = I<int>();
    var r2 = I<int>();
    var z = F((x1, y1, r1), (x2, y2, r2));
    ans.Add(z.ToString("F9"));
}
Console.WriteLine(string.Join(Environment.NewLine, ans));
0