#137721 (C++14) No.461 三角形はいくつ？

提出ソース

結果

問題	No.461 三角形はいくつ？
ユーザー	chaemon
提出日時	2016-12-12 00:51:55
言語	C++14 (gcc 12.3.0 + boost 1.83.0)
結果	AC
実行時間	2,108 ms / 5,000 ms
コード長	3,627 bytes
コンパイル時間	993 ms
コンパイル使用メモリ	98,780 KB
実行使用メモリ	5,248 KB
最終ジャッジ日時	2024-11-29 04:34:01
合計ジャッジ時間	34,481 ms
ジャッジサーバーID （参考情報）	judge4 / judge3

このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示

入力	結果	実行時間実行使用メモリ
testcase_00	AC	2 ms 5,248 KB
testcase_01	AC	2 ms 5,248 KB
testcase_02	AC	2 ms 5,248 KB
testcase_03	AC	2 ms 5,248 KB
testcase_04	AC	5 ms 5,248 KB
testcase_05	AC	945 ms 5,248 KB
testcase_06	AC	514 ms 5,248 KB
testcase_07	AC	629 ms 5,248 KB
testcase_08	AC	426 ms 5,248 KB
testcase_09	AC	4 ms 5,248 KB
testcase_10	AC	133 ms 5,248 KB
testcase_11	AC	719 ms 5,248 KB
testcase_12	AC	319 ms 5,248 KB
testcase_13	AC	1,061 ms 5,248 KB
testcase_14	AC	1,017 ms 5,248 KB
testcase_15	AC	1,076 ms 5,248 KB
testcase_16	AC	729 ms 5,248 KB
testcase_17	AC	766 ms 5,248 KB
testcase_18	AC	993 ms 5,248 KB
testcase_19	AC	1,016 ms 5,248 KB
testcase_20	AC	1,116 ms 5,248 KB
testcase_21	AC	1,067 ms 5,248 KB
testcase_22	AC	1,056 ms 5,248 KB
testcase_23	AC	570 ms 5,248 KB
testcase_24	AC	592 ms 5,248 KB
testcase_25	AC	2,101 ms 5,248 KB
testcase_26	AC	2,103 ms 5,248 KB
testcase_27	AC	6 ms 5,248 KB
testcase_28	AC	7 ms 5,248 KB
testcase_29	AC	2,007 ms 5,248 KB
testcase_30	AC	2,108 ms 5,248 KB
testcase_31	AC	1,797 ms 5,248 KB
testcase_32	AC	1,665 ms 5,248 KB
testcase_33	AC	7 ms 5,248 KB
testcase_34	AC	7 ms 5,248 KB
testcase_35	AC	6 ms 5,248 KB
testcase_36	AC	597 ms 5,248 KB
testcase_37	AC	580 ms 5,248 KB
testcase_38	AC	593 ms 5,248 KB
testcase_39	AC	583 ms 5,248 KB
testcase_40	AC	597 ms 5,248 KB
testcase_41	AC	604 ms 5,248 KB
testcase_42	AC	688 ms 5,248 KB
testcase_43	AC	668 ms 5,248 KB
testcase_44	AC	711 ms 5,248 KB

権限があれば一括ダウンロードができます

ソースコード

raw source code

// #includes {{{
#include <algorithm>
#include <numeric>
#include <iostream>
#include <string>
#include <vector>
#include <queue>
#include <list>
#include <deque>
#include <stack>
#include <set>
#include <map>
#include <cstdio>
#include <cstdlib>
#include <cassert>
#include <cstring>
#include <cmath>
using namespace std;
// }}}
// pre-written code {{{
#define REP(i,n) for(int i=0;i<(int)(n);++i)
#define RREP(i,a,b) for(int i=(int)(a);i<(int)(b);++i)
#define FOR(i,c) for(__typeof((c).begin()) i=(c).begin();i!=(c).end();++i)
#define LET(x,a) __typeof(a) x(a)
//#define IFOR(i,it,c) for(__typeof((c).begin())it=(c).begin();it!=(c).end();++it,++i)
#define ALL(c) (c).begin(), (c).end()
#define MP make_pair

#define EXIST(e,s) ((s).find(e)!=(s).end())

#define RESET(a) memset((a),0,sizeof(a))
#define SET(a) memset((a),-1,sizeof(a))
#define PB push_back
#define DEC(it,command) __typeof(command) it=command

//debug
#define dump(x)  cerr << #x << " = " << (x) << endl;
#define debug(x) cerr << #x << " = " << (x) << " (L" << __LINE__ << ")" << " " << __FILE__ << endl;
#define debug2(x) cerr << #x << " = [";REP(__ind,(x).size()){cerr << (x)[__ind] << ", ";}cerr << "] (L" << __LINE__ << ")" << endl;

const int INF=0x3f3f3f3f;

typedef long long Int;
typedef unsigned long long uInt;
#ifdef __MINGW32__
typedef double rn;
#else
typedef long double rn;
#endif

typedef pair<int,int> pii;

/*
#ifdef MYDEBUG
#include"debug.h"
#include"print.h"
#endif
*/
// }}}

//{{{ quatient class
Int gcd(Int a, Int b) { return a > 0 ? gcd(b % a, a) : b; }
struct F {
	Int p, q;
	void normalize() { // keep q positive
		if (q < 0) p *= -1, q *= -1;
		Int d = gcd(p < 0 ? -p : p, q);
		if (d == 0) p = 0,  q = 1;
		else        p /= d, q /= d;
	}
	F(Int p, Int q = 1) : p(p), q(q) {
		normalize();
	}
	F &operator += (const F &a) {
		p = a.q * p + a.p * q; q = a.q * q;
		return *this;
	}
	F &operator -= (const F &a) {
		p = a.q * p - a.p * q; q = a.q * q;
		return *this;
	}
	F &operator *= (const F &a) {
		p *= a.p; q *= a.q;
		return *this;
	}
	F &operator /= (const F &a) {
		p *= a.q; q *= a.p;
		return *this;
	}
	F &operator - () {
		p *= -1;
		return *this;
	}
};
F operator + (const F &a, const F &b) {
	return F(a) += b;
}
F operator * (const F &a, const F &b) {
	return F(a) *= b;
}
F operator - (const F &a, const F &b) {
	return F(a) -= b;
}
F operator / (const F &a, const F &b) {
	return F(a) /= b;
}
bool operator < (const F &a, const F &b) { // avoid overflow
	return a.p * b.q < a.q * b.p;
}
bool operator <= (const F &a, const F &b) {
	return !(b < a);
}
bool operator > (const F &a, const F &b) {
	return b < a;
}
bool operator >= (const F &a, const F &b) {
	return !(a < b);
}
bool operator == (const F &a, const F &b) {
	return !(a < b) && !(b < a);
}
bool operator != (const F &a, const F &b) {
	return (a < b) || (b < a);
}
//}}}

int N;
vector<F> t[3];
vector<F> listed;

int main(){
	cin>>N;
	REP(i,3)t[i].push_back(1.0L);
	listed.push_back(1.0L);
	REP(i,N){
		int p,a,b;
		cin>>p>>a>>b;
		t[p].push_back(F(a,a+b));
		if(p==2)listed.push_back(t[p].back());
	}
	sort(ALL(listed));
	Int d = 0;
	REP(i,t[0].size())REP(j,t[1].size()){
		F u = F(2)-(t[0][i] + t[1][j]);
		u.normalize();
		d+=distance(lower_bound(ALL(listed),u),upper_bound(ALL(listed),u));
	}
	Int ans = 0;
	REP(i,t[0].size())REP(j,t[1].size()){
		if(t[0][i]+t[1][j]<F(1))continue;
		F u = max(F(1)-t[0][i],F(1)-t[1][j]);
		u.normalize();
		
		Int k = distance(lower_bound(ALL(listed),u),listed.end());
//		cout<<t[0][i]<<" "<<t[1][j]<<" "<<u<<" "<<k<<endl;
		ans+=k;
	}
//	cout<<ans<<" "<<d<<endl;
//	cout<<ans<<endl;
	cout<<ans-d<<endl;
	return 0;
}

yukicoder

結果

テストケース

ソースコード