結果
問題 | No.306 さいたま2008 |
ユーザー |
![]() |
提出日時 | 2020-03-18 01:20:41 |
言語 | C#(csc) (csc 3.9.0) |
結果 |
AC
|
実行時間 | 31 ms / 2,000 ms |
コード長 | 3,533 bytes |
コンパイル時間 | 981 ms |
コンパイル使用メモリ | 116,480 KB |
実行使用メモリ | 26,908 KB |
最終ジャッジ日時 | 2025-02-14 12:16:14 |
合計ジャッジ時間 | 2,625 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
other | AC * 26 |
コンパイルメッセージ
Microsoft (R) Visual C# Compiler version 3.9.0-6.21124.20 (db94f4cc) Copyright (C) Microsoft Corporation. All rights reserved.
ソースコード
using System;using System.Collections.Generic;using System.Collections;using System.Collections.Specialized;using System.Linq;using System.Text;using System.IO;using System.Reflection;using static System.Math;using System.Numerics;static class Program{const int mod=(int)1e9+7;const decimal eps=(decimal)1e-11;static int[] a,b;static void Main(){Sc sc=new Sc();a=sc.Ia;b=sc.Ia;if(a[1]>b[1]){a[0]^=b[0];b[0]^=a[0];a[0]^=b[0];a[1]^=b[1];b[1]^=a[1];a[1]^=b[1];}Console.WriteLine(Gss(a[1]-1,b[1]+1));}static readonly decimal phi=(1+Sqrt(5))*(decimal)0.5;static decimal Gss(decimal a,decimal b){decimal lb=a,ub=b,c1=0,c2=0,d1=0,d2=0;c1=(lb*phi+ub)/(1+phi);c2=(lb+phi*ub)/(1+phi);bool bo=true;d1=Fu(c1);while(ub-lb>eps){if(bo){d2=Fu(c2);}else{d1=Fu(c1);}if(d1>d2){lb=c1;c1=c2;c2=(lb+phi*ub)/(1+phi);d1=d2;bo=true;}else{ub=c2;c2=c1;c1=(lb*phi+ub)/(1+phi);d2=d1;bo=false;}}return c1;}static decimal Fu(decimal c){return Sqrt(a[0]*a[0]+(a[1]-c)*(a[1]-c))+Sqrt(b[0]*b[0]+(b[1]-c)*(b[1]-c));}static decimal Sqrt(decimal x,decimal? guess=null){var ourGuess=guess.GetValueOrDefault(x/2m);var result=x/ourGuess;var average=(ourGuess+result)/2m;if(average==ourGuess){return average;}else{return Sqrt(x,average);}}}public class Sc{public int I{get{return int.Parse(Console.ReadLine());}}public long L{get{return long.Parse(Console.ReadLine());}}public double D{get{return double.Parse(Console.ReadLine());}}public string S{get{return Console.ReadLine();}}public int[] Ia{get{return Array.ConvertAll(Console.ReadLine().Split(),int.Parse);}}public long[] La{get{return Array.ConvertAll(Console.ReadLine().Split(),long.Parse);}}public double[] Da{get{return Array.ConvertAll(Console.ReadLine().Split(),double.Parse);}}public string[] Sa{get{return Console.ReadLine().Split();}}public object[] Oa{get{return Console.ReadLine().Split();}}public int[] Ia2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),int.Parse);}}public int[] Ia3(int a){return Array.ConvertAll((a.ToString()+" "+Console.ReadLine()).Split(),int.Parse);}public int[] Ia3(bool a,int b,bool c,int d){return Array.ConvertAll(((a?b.ToString()+" ":"")+Console.ReadLine()+(c?" "+d.ToString():"")).Split(),int.Parse);}public long[] La2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),long.Parse);}}public long[] La3(int a){return Array.ConvertAll((a.ToString()+" "+Console.ReadLine()).Split(),long.Parse);}public long[] La3(bool a,int b,bool c,int d){return Array.ConvertAll(((a?b.ToString()+" ":"")+Console.ReadLine()+(c?" "+d.ToString():"")).Split(),long.Parse);}public double[] Da2{get{return Array.ConvertAll(("0 "+Console.ReadLine()+" 0").Split(),double.Parse);}}public double[] Da3(int a){return Array.ConvertAll((a.ToString()+" "+Console.ReadLine()).Split(),double.Parse);}public double[] Da3(bool a,int b,bool c,int d){return Array.ConvertAll(((a?b.ToString()+" ":"")+Console.ReadLine()+(c?" "+d.ToString():"")).Split(),double.Parse);}public T[] Arr<T>(int n,Func<T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f();}return a;}public T[] Arr<T>(int n,Func<int,T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(i);}return a;}public T[] Arr<T>(int n,Func<string[],T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(Console.ReadLine().Split());}return a;}public T[] Arr<T>(int n,Func<int,string[],T> f){var a=new T[n];for(int i=0;i<n;i++){a[i]=f(i,Console.ReadLine().Split());}return a;}}