結果

問題 No.1143 面積Nの三角形
ユーザー xyz2606xyz2606
提出日時 2020-08-01 00:30:09
言語 C++14
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 188 ms / 800 ms
コード長 3,090 bytes
コンパイル時間 1,968 ms
コンパイル使用メモリ 115,440 KB
実行使用メモリ 13,012 KB
最終ジャッジ日時 2023-09-21 04:14:41
合計ジャッジ時間 3,499 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 1 ms
4,376 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 AC 2 ms
4,376 KB
testcase_07 AC 11 ms
4,380 KB
testcase_08 AC 14 ms
4,376 KB
testcase_09 AC 22 ms
4,472 KB
testcase_10 AC 11 ms
4,380 KB
testcase_11 AC 82 ms
8,144 KB
testcase_12 AC 102 ms
9,140 KB
testcase_13 AC 106 ms
9,564 KB
testcase_14 AC 12 ms
4,380 KB
testcase_15 AC 2 ms
4,376 KB
testcase_16 AC 185 ms
13,008 KB
testcase_17 AC 188 ms
13,012 KB
testcase_18 AC 156 ms
11,284 KB
testcase_19 AC 148 ms
11,160 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<vector>
#include<set>
#include<map>
#include<queue>
#include<string>
#include<algorithm>
#include<iostream>
#include<bitset>
#include<functional>
#include<numeric>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<cassert>
#include<cmath>
#include<iomanip>
#include<random>
#include<ctime>
#include<complex>
using namespace std;
typedef long long LL;
typedef double D;
#define all(v) (v).begin(), (v).end()
mt19937 gene(233);
typedef complex<double> Complex;
#define fi first
#define se second
#define ins insert
#define pb push_back
inline char GET_CHAR(){
	const int maxn = 131072;
	static char buf[maxn],*p1=buf,*p2=buf;
	return p1==p2&&(p2=(p1=buf)+fread(buf,1,maxn,stdin),p1==p2)?EOF:*p1++;
}
inline int getInt() {
	int res(0);
	char c = getchar();
	while(c < '0') c = getchar();
	while(c >= '0') {
		res = res * 10 + (c - '0');
		c = getchar();
	}
	return res;
}

inline LL fastpo(LL x, LL n, LL mod) {
	LL res(1);
	while(n) {
		if(n & 1) {
			res = res * (LL)x % mod;
		}
		x = x * (LL) x % mod;
		n /= 2;
	}
	return res;
}
LL gcd(LL a, LL b) { return b ? gcd(b, a % b) : a; }

inline string itoa(int x, int width = 0) {
  string res;
  if(x == 0) res.push_back('0');
  while(x) {
    res.push_back('0' + x % 10);
    x /= 10;
  }
  while((int)res.size() < width) res.push_back('0');
  reverse(res.begin(), res.end());
  return res;
}
const int _B = 131072;
char buf[_B];
int _bl = 0;
inline void flush() {
	fwrite(buf, 1, _bl, stdout);
	_bl = 0;
}
__inline void _putchar(char c) {
	if(_bl == _B) flush();
	buf[_bl++] = c;
}
inline void print(LL x, char c) {
	static char tmp[20];
	int l = 0;
	if(!x) tmp[l++] = '0';
	else {
		while(x) {
			tmp[l++] = x % 10 + '0';
			x /= 10;
		}
	}
	for(int i = l - 1; i >= 0; i--) _putchar(tmp[i]);
	_putchar(c);
}
struct P {
	D x, y;
};
const int N = 300033;
const int LOG = 20;
const int mod = 1e9 + 7;
const int inf = 1e9 + 7;
int n, m;
int dx[4] = {1, 0, -1, 0};
int dy[4] = {0, 1, 0, -1};
int main() {
	LL n;
	vector<LL> v, u;
	cin >> n;
	for(int p = 1; p <= n; p++) {
		if(n % p == 0) v.push_back(p);
	}
	for(int i = 0; i < (int)v.size(); i++) {
		for(int j = 0; j <= i; j++) {
			u.push_back(v[i] * v[j]);
		}
	}
	sort(u.begin(), u.end());
	u.resize(unique(u.begin(), u.end()) - u.begin());
	int ans = 0;
	//cout << u.size() << endl;
	vector<vector<LL>> w(u.size());
	for(int i = 0; i < (int)u.size(); i++) {
		for(int  j = 0; j <= i; j++) {
			if(u[i] % u[j] == 0) {
				w[i].pb(u[j]);
			}
		}
	}
	for(int id = 0; id < (int)u.size(); id++) {
		LL p = u[id];
		LL nnp = lower_bound(u.begin(), u.end(), n * n / p) - u.begin();
		for(LL pa : w[nnp]) {
			if(pa >= p) break;
			LL a = p - pa;
			//cout << p << ' ' << a << endl;
			LL nb = lower_bound(u.begin(), u.end(), n * n / p / (p - a)) - u.begin();
			for(LL pb : w[nb]) {
				
				LL b = p - pb;
				if(b < max(a, p - a + 1) || b >= p) continue;
				//cout << p << ' ' << a << ' ' << b << endl;
				LL c = p * 2 - a - b;
				if(c < a || c < b) continue;
				if(n * n / p / (p - a) / (p - b) != p - c) continue;
				ans++;
			}
		}
	}
	cout << ans << endl;
}
0