結果
問題 | No.245 貫け! |
ユーザー |
![]() |
提出日時 | 2015-07-17 23:31:37 |
言語 | C++11 (gcc 13.3.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 4,567 bytes |
コンパイル時間 | 905 ms |
コンパイル使用メモリ | 96,724 KB |
実行使用メモリ | 5,376 KB |
最終ジャッジ日時 | 2024-07-08 09:37:13 |
合計ジャッジ時間 | 1,854 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 4 |
other | AC * 4 WA * 12 |
ソースコード
#include <iostream>#include <iomanip>#include <sstream>#include <vector>#include <string>#include <set>#include <unordered_set>#include <map>#include <unordered_map>#include <stack>#include <queue>#include <deque>#include <algorithm>#include <functional>#include <iterator>#include <limits>#include <numeric>#include <utility>#include <cmath>#include <cassert>#include <cstdio>using namespace std; using namespace placeholders;using LL = long long;using ULL = unsigned long long;using VI = vector< int >;using VVI = vector< vector< int > >;using VS = vector< string >;using SS = stringstream;using PII = pair< int, int >;using VPII = vector< pair< int, int > >;template < typename T = int > using VT = vector< T >;template < typename T = int > using VVT = vector< vector< T > >;template < typename T = int > using LIM = numeric_limits< T >;template < typename T > inline istream& operator>>( istream &s, vector< T > &v ){ for ( T &t : v ) { s >> t; } return s; }template < typename T > inline ostream& operator<<( ostream &s, const vector< T > &v ){ for ( int i = 0; i < int( v.size() ); ++i ){ s << ( " " + !i) << v[i]; } return s; }template < typename T > inline T fromString( const string &s ) { T res; istringstream iss( s ); iss >> res; return res; };template < typename T > inline string toString( const T &a ) { ostringstream oss; oss << a; return oss.str(); };#define REP2( i, n ) REP3( i, 0, n )#define REP3( i, m, n ) for ( int i = ( int )( m ); i < ( int )( n ); ++i )#define GET_REP( a, b, c, F, ... ) F#define REP( ... ) GET_REP( __VA_ARGS__, REP3, REP2 )( __VA_ARGS__ )#define FOR( e, c ) for ( auto &e : c )#define ALL( c ) ( c ).begin(), ( c ).end()#define AALL( a, t ) ( t* )a, ( t* )a + sizeof( a ) / sizeof( t )#define DRANGE( c, p ) ( c ).begin(), ( c ).begin() + ( p ), ( c ).end()#define SZ( v ) ( (int)( v ).size() )#define PB push_back#define EM emplace#define EB emplace_back#define BI back_inserter#define EXIST( c, e ) ( ( c ).find( e ) != ( c ).end() )#define MP make_pair#define fst first#define snd second#define DUMP( x ) cerr << #x << " = " << ( x ) << endl#include <complex>using Point = complex< double >;constexpr double EPS = 1e-8;constexpr double PI = acos( -1 );constexpr Point O( 0, 0 );// 入力ストリームから実数二つをとって Point へistream& operator>> ( istream &s, Point &a ){double r, i;s >> r >> i;a = Point( r, i );return s;}// 内積(ドット積)double dot( const Point &a, const Point &b ){return a.real() * b.real() + a.imag() * b.imag();}// 外積(クロス積)double cross( const Point &a, const Point &b ){return a.real() * b.imag() - a.imag() * b.real();}// 点 p が線分 ( q1, q2 ) に乗っているか// include : dot, crossbool intersectPS( const Point &p, const Point &q1, const Point &q2 ){// 端点からの p へのベクトルの外積が 0 ⇔ 端点からのベクトルが並行 ⇔ 直線 ( q1, q2 ) に点 p が乗る// p から端点へのベクトルの内積が 0 以下 ⇔ 点 p が q1, q2 の内側にあるreturn abs( cross( q1 - p, q2 - p ) ) <= EPS && dot( q1 - p, q2 - p ) <= EPS;}// 直線 ( p1, p2 ) と直線 ( q1, q2 ) の交点// include : crossPoint intersectionLL( const Point &p1, const Point &p2, const Point &q1, const Point &q2 ){return p1 + ( p2 - p1 ) * ( cross( q2 - q1, q1 - p1 ) / cross( q2 - q1, p2 - p1 ) );}// 線分 ( p1, p2 ) と線分 ( q1, q2 ) が交差するか// include : cross, intersectPS, intersectionLLbool intersectSS( const Point &p1, const Point &p2, const Point &q1, const Point &q2 ){if ( abs( cross( p1 - p2, q1 - q2 ) ) <= EPS ) // 二直線が並行な場合{// 片方の線分に他方の端点が乗るかreturn intersectPS( q1, p1, p2 ) ||intersectPS( q2, p1, p2 ) ||intersectPS( p1, q1, q2 ) ||intersectPS( p2, q1, q2 );}else{// 二直線の交点が両方の線分に乗るかPoint r( intersectionLL( p1, p2, q1, q2 ) );return intersectPS( r, p1, p2 ) && intersectPS( r, q1, q2 );}}int main(){cin.tie( 0 );ios::sync_with_stdio( false );int N;cin >> N;VT< Point > ps( N * 2 );cin >> ps;int res = 0;REP( i, N * 2 ){REP( j, i ){const Point a = ( ps[j] - ps[i] ) * Point( 1000000 );const Point p = ps[i] - a, q = ps[i] + a;int tmp = 0;for ( int k = 0; k < 2 * N; k += 2 ){tmp += intersectSS( p, q, ps[k], ps[ k + 1 ] );}res = max( res, tmp );}}cout << res << endl;return 0;}