結果

問題 No.306 さいたま2008
ユーザー Goryudyuma
提出日時 2015-12-08 23:30:56
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
WA  
実行時間 -
コード長 9,898 bytes
コンパイル時間 1,502 ms
コンパイル使用メモリ 169,648 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-14 20:26:30
合計ジャッジ時間 2,331 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 22 WA * 1
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <bits/stdc++.h>
using namespace std;
//提出時に必ず0とし、確認実行すること！！
#define DEBUG 1
class D3
{
public:
    long double X, Y, Z, EPS;
    D3( long double, long double, long double );
    bool operator== ( D3 Partner );
    bool operator==( long double );
    D3 operator*( long double );
    D3 operator/( long double );
    friend ostream& operator<<( ostream& os, const D3& );
};
class Point :private D3
{
public:
    Point( long double, long double, long double );
    long double getX();
    long double getY();
    long double getZ();
    bool operator== ( Point );
    bool operator<( Point );
    bool operator>( Point );
    friend bool operator<( const Point&, const Point& );
    friend bool operator>( const Point&, const Point& );
    long double S_point( Point, Point );
    long double V_point( Point, Point, Point );
    friend ostream& operator<<( ostream& os, const Point& );
};
class Vec
{
private:
    Point SP;//Starting Point 始点
    D3 D;//Direction 方向
    Vec( D3 Direction_, Point SP_ );
public:
    Vec( long double x_, long double y_, long double z_, Point );
    Vec operator+( Vec Partner );
    Vec operator*( Vec Partner );
    Vec operator*( long double ld );
    Vec operator-( Vec Partner );
    Vec operator/( Vec Partner );
    Vec operator/( long double );
    bool operator== ( Vec );
    bool operator==( long double );
    Vec Cross_product( Vec );
    double Inner_product( Vec );
    Vec( Point A, Point B );
    Point getSP();
    Point getGP();
    D3 getD();
    long double length();
    Vec Unit_vec();
    Vec Inverse_vec();
    Vec Reverse_vec();
    bool Vertical( Vec );
    bool Parallel( Vec );
    bool isIntersection( Vec );
    Point IntersectionPoint( Vec );
    bool onVecPoint( Point );
    long double S_vec( Vec );
    long double V_vec( Vec, Vec );
    friend ostream& operator<<( ostream& os, const Vec& );
};
//D3
D3::D3( long double x_ = 0.0L, long double y_ = 0.0L, long double z_ = 0.0L ) {
    X = x_;
    Y = y_;
    Z = z_;
    EPS = 1e-9L;
}
bool D3::operator== ( D3 Partner ) {
    return abs( X - Partner.X ) < EPS&&abs( Y - Partner.Y ) < EPS&&abs( Z - Partner.Z ) < EPS;
}
bool D3::operator==( long double ld ) {
    return (*this) == D3();
}
D3 D3::operator*( long double ld ) {
    return D3( X*ld, Y*ld, Z*ld );
}
D3 D3::operator/( long double ld ) {
    return  *this*(ld*-1.L);
}
ostream& operator<<( ostream& os, const D3& d ) {
    if( DEBUG ) {
        os << "D3" << endl;
        os << "X:" << d.X << endl;
        os << "Y:" << d.Y << endl;
        os << "Z:" << d.Z << endl;
    }
    return os;
}
//Point
Point::Point( long double x_ = 0.0L, long double y_ = 0.0L, long double z_ = 0.0L ) {
    X = x_;
    Y = y_;
    Z = z_;
}
//x座標を返す
long double Point::getX() {
    return X;
}
//y座標を返す
long double Point::getY() {
    return Y;
}
//z座標を返す
long double Point::getZ() {
    return Z;
}
bool Point::operator== ( Point Partner ) {
    return D3( *this ) == D3( Partner );
}
bool Point::operator<( Point Partner ) {
    if( (*this).getX() != Partner.getX() ) {
        return (*this).getX() < Partner.getX();
    }
    if( (*this).getY() != Partner.getY() ) {
        return (*this).getY() < Partner.getY();
    }
    return (*this).getZ() < Partner.getZ();
}
bool Point::operator>( Point Partner ) {
    return Partner < (*this);
}
//三点から面積を求める
long double Point::S_point( Point B, Point C ) {
    Vec AB( *this, B ), AC( *this, C );
    return AB.S_vec( AC );
}
//四点から体積を求める
long double Point::V_point( Point B, Point C, Point D ) {
    return Vec( (*this), B ).V_vec( Vec( (*this), C ), Vec( (*this), D ) );
}
ostream & operator<<( ostream & os, const Point &p ) {
    if( DEBUG ) {
        os << "Point" << endl;
        os << "X:" << p.X << endl;
        os << "Y:" << p.Y << endl;
        os << "Z:" << p.Z << endl;
    }
    return os;
}
bool operator<( const Point& p1, const Point& p2 ) {
    if( p1.X == p2.X ) {
        if( p1.Y == p2.Y ) {
            return p1.Z < p2.Z;
        }
        return p1.Y < p2.Y;
    }
    return p1.X < p2.X;
}
bool operator>( const Point& p1, const Point& p2 ) {
    return p2 < p1;
}
//Vec
//ベクトルと始点からベクトルを作る
Vec::Vec( D3 Direction_, Point SP_ ) {
    D = Direction_;
    SP = SP_;
}
//始点と各ベクトルの長さからベクトルを作る
Vec::Vec( long double x_ = 0.0L, long double y_ = 0.0L, long double z_ = 0.0L, Point SP_ = Point() ) {
    *this = Vec( D3( x_, y_, z_ ), SP_ );
}
Vec Vec::operator+( Vec Partner ) {
    return Vec( D.X + Partner.D.X, D.Y + Partner.D.Y, D.Z + Partner.D.Z, SP );
}
Vec Vec::operator*( Vec Partner ) {
    return Vec( D.X * Partner.D.X, D.Y * Partner.D.Y, D.Z * Partner.D.Z, SP );
}
Vec Vec::operator*( long double ld ) {
    return Vec( D*ld, SP );
}
Vec Vec::operator-( Vec Partner ) {
    return *this + (Partner*(-1.L));
}
Vec Vec::operator/( Vec Partner ) {
    return Vec( D.X / Partner.D.X, D.Y / Partner.D.Y, D.Z / Partner.D.Z, SP );
}
Vec Vec::operator/( long double ld ) {
    return *this*(1.L / ld);
}
bool Vec::operator== ( Vec Partner ) {
    return D == Partner.D&&SP == Partner.SP;
}
bool Vec::operator==( long double ) {
    return (*this) == Vec();
}
//外積
Vec Vec::Cross_product( Vec Partner ) {
    return Vec( D.Y*Partner.getD().Z - D.Z*Partner.getD().Y, D.Z*Partner.getD().X - D.X*Partner.getD().Z, D.X*Partner.getD().Y - D.Y*Partner.getD().X 
        );
}
//内積
double Vec::Inner_product( Vec Partner ) {
    return D.X*Partner.getD().X + D.Y*Partner.getD().Y + D.Z*Partner.getD().Z;
}
//２つの点からベクトルを作る
Vec::Vec( Point A, Point B ) {
    *this = Vec( B.getX() - A.getX(), B.getY() - A.getY(), B.getZ() - A.getZ(), A );
}
//始点を返す
Point Vec::getSP() {
    return SP;
}
//終点を返す
Point Vec::getGP() {
    return Point( SP.getX() + D.X, SP.getY() + D.Y, SP.getZ() + D.Z );
}
//方向を返す
D3 Vec::getD() {
    return D;
}
//長さを返す
long double Vec::length() {
    return sqrtl( D.X*D.X + D.Y*D.Y + D.Z*D.Z );
}
//単位ベクトル(長さ1)を返す
Vec Vec::Unit_vec() {
    return Vec( *this ) / (*this).length();
}
//逆ベクトルを返す(始点は変わらず。)
Vec Vec::Inverse_vec() {
    return (*this)*-1.L;
}
//逆ベクトルを返す(始点と終点の入れ替え。)
Vec Vec::Reverse_vec() {
    return Vec( (*this).Inverse_vec().getD(), (*this).getGP() );
}
//垂直かどうか
bool Vec::Vertical( Vec Partner ) {
    return (*this).Inner_product( Partner ) == 0;
}
//平行かどうか
bool Vec::Parallel( Vec Partner ) {
    return (*this).Unit_vec().getD() == Partner.Unit_vec().getD() || (*this).Unit_vec().Inverse_vec().getD() == Partner.Unit_vec().getD();
}
//２つのベクトルが交差しているかどうか
bool Vec::isIntersection( Vec Partner ) {
    int count = 0;
    Vec O = Partner.Cross_product( Vec( Partner.getSP(), (*this).getSP() ) );
    Vec P = Partner.Cross_product( Vec( Partner.getSP(), (*this).getGP() ) );
    Vec Q = (*this).Cross_product( Vec( (*this).getSP(), Partner.getSP() ) );
    Vec R = (*this).Cross_product( Vec( (*this).getSP(), Partner.getGP() ) );
    return (!(*this).Parallel( Partner )) && (
        O.getD().X*P.getD().X < D3().EPS&& O.getD().Y*P.getD().Y < D3().EPS&& O.getD().Z*P.getD().Z < D3().EPS) && (
            Q.getD().X*R.getD().X < D3().EPS&& Q.getD().Y*R.getD().Y < D3().EPS&& Q.getD().Z*R.getD().Z < D3().EPS);
}
//ベクトル同士の交点を返す
//必ず交点があるかどうかを確認すること！
//TODO:未実装
Point Vec::IntersectionPoint( Vec ) {
    return Point();
}
//ベクトル上にポイントが有るかどうかを返す
bool Vec::onVecPoint( Point A ) {
    Vec Test( (*this).getSP(), A );
    if( A == ((*this).getSP()) || (Test.Parallel( *this ) && (*this).getD().X / Test.getD().X >= 1.0L) ) {
        return true;
    }
    return false;
}
//同一始点2ベクトルから面積を求める
long double Vec::S_vec( Vec B ) {
    if( ((*this).getSP() == B.getSP()) ) {
        Vec G = (*this).Cross_product( B );
        return sqrtl( G.Inner_product( G ) ) / 2.L;
    }
    if( (*this).getGP() == B.getSP() ) {
        return (*this).Reverse_vec().S_vec( B );
    }
    if( (*this).getSP() == B.getGP() ) {
        return (*this).S_vec( B.Reverse_vec() );
    }
    if( (*this).getGP() == B.getGP() ) {
        return (*this).Reverse_vec().S_vec( B.Reverse_vec() );
    }
    return -1.L;
}
//同一始点かもしくは終点と始点が繋がっていなくてはいけない
//しかし今のところその判定はなし。
//TODO 始点と終点の同一判定
long double Vec::V_vec( Vec B, Vec C ) {
    Vec G = (*this).Cross_product( B );
    return G.Inner_product( C ) / 6.L;
    return 0;
}
ostream & operator<<( ostream & os, const Vec &v ) {
    if( DEBUG ) {
        os << "Vec" << endl;
        os << "start point" << endl;
        os << v.SP << endl;
        os << "vec" << endl;
        os << v.D << endl;
    }
    return os;
}
class Points
{
private:
    vector<Point>VP;
public:
    void push( Point p ) {
        VP.push_back( p );
    }
    bool erase( Point p ) {
        auto now = find( VP.begin(), VP.end(), p );
        if( now != VP.end() ) {
            VP.erase( now );
            return 0;
        } else {
            return 1;
        }
    }
    void print() {
        for( auto x : VP ) {
            cout << x << endl;
        }
    }
    //sortします。Xの大小->Yの大小->Zの大小。
    void sort() {
        std::sort( VP.begin(), VP.end() );
    }
};
class PointDistance
{
public:
    vector <int> findPoint( int x1, int y1, int x2, int y2 ) {
        Point A( x1, y1 ), B( x2, y2 );
        for( int i = -100; i <= 100; i++ ) {
            for( int j = -100; j <= 100; j++ ) {
                if( Vec( A, Point( i, j ) ).length() > Vec( B, Point( i, j ) ).length() && Vec( B, Point( i, j ) ).length() >= 1 ) {
                    return vector<int>{i, j };
                }
            }
        }
    }
};
int main() {
    long double x, y;
    cin >> x >> y;
    Point A( x, y );
    cin >> x >> y;
    Point B( x, y );
    long double mi = 0.0L, ma = 1000.0L;
    for( size_t i = 0; i < 10000; i++ ) {
        long double mid = (mi + ma) / 2.L;
        long double midadd = mid + 1e-7;
        if( Vec( A, Point( 0, mid ) ).length() + Vec( Point( 0, mid ), B ).length() < Vec( A, Point( 0, midadd ) ).length() + Vec( Point( 0, midadd 
            ), B ).length() ) {
            ma = mid;
        } else {
            mi = mid;
        }
    }
    cout << fixed << setprecision( 20 ) << mi << endl;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0