結果

問題 No.152 貯金箱の消失
コンテスト
ユーザー komori3
提出日時 2017-10-20 05:58:13
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 13 ms / 5,000 ms
コード長 2,768 bytes
コンパイル時間 741 ms
コンパイル使用メモリ 78,676 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-11-21 08:15:59
合計ジャッジ時間 1,623 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
other AC * 12
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _CRT_SECURE_NO_WARNINGS
#include <iostream>
#include <climits>
#include <cmath>
#include <string>
#include <vector>
#include <algorithm>
#include <queue>
#include <map>
#include <functional>
#include <set>
#include <numeric>
#include <stack>
#include <utility>
#include <time.h>
#include <iterator>
//#include "util.h"

using namespace std;
typedef unsigned uint;
typedef long long lint;
typedef unsigned long long ulint;



//呪文
template <typename _KTy, typename _Ty> ostream& operator << (ostream& ostr, const pair<_KTy, _Ty>& m) { cout << "{" << m.first << ", " << m.second << "}"; return ostr; }
template <typename _KTy, typename _Ty> ostream& operator << (ostream& ostr, const map<_KTy, _Ty>& m) { if (m.empty()) { cout << "{ }"; return ostr; } cout << "{" << *m.begin(); for (auto itr = ++m.begin(); itr != m.end(); itr++) { cout << ", " << *itr; } cout << "}"; return ostr; }
template <typename _Ty> ostream& operator << (ostream& ostr, const vector<_Ty>& v) { if (v.empty()) { cout << "{ }"; return ostr; }	cout << "{" << v.front(); for (auto itr = ++v.begin(); itr != v.end(); itr++) { cout << ", " << *itr; }	cout << "}"; return ostr; }
template <typename _Ty> ostream& operator << (ostream& ostr, const set<_Ty>& s) { if (s.empty()) { cout << "{ }"; return ostr; } cout << "{" << *(s.begin()); for (auto itr = ++s.begin(); itr != s.end(); itr++) { cout << ", " << *itr; }	cout << "}"; return ostr; }
//template <typename T> void print(T* v, int N) { if (N == 0) cout << "{ }"; cout << "{" << v[0]; for (int i = 1; i < N; i++) cout << ", " << v[i]; cout << "}"; }
#define PI 3.14159265358979323846
#define EPS 1e-8
#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n)  FOR(i,0,n)
#define all(x) (x).begin(), (x).end()



// 互いに素
int coprime(int a, int b)
{
	if (a == 1 || b == 1) return true;
	int tmp = 0;
	while (tmp != 1) {
		tmp = a % b;
		if (tmp == 0) return false;
		a = b;
		b = tmp;
	}
	return true;
}

// 周長 L 以下の原始ピタゴラス数の総数
int count_primitive_Pythagorean_triple(int L)
{
	int res = 0;
	for (int m = 1; 2 * (m + 1) * (2 * m + 1) <= L; m++) {
		for (int n = m + 1; 2 * n * (m + n) <= L; n += 2) {
			if (coprime(m, n)) res++;
		}
	}
	return res;
}

// 周長 L 以下のピタゴラス数の総数
int count_Pythagorean_triple(int L)
{
	int res = 0;
	for (int m = 1; 2 * (m + 1) * (2 * m + 1) <= L; m++) {
		for (int n = m + 1; 2 * n * (m + n) <= L; n += 2) {
			if (coprime(m, n)) res += L / (2 * n * (m + n));
		}
	}
	return res;
}

int yuki0152()
{
	int L;
	cin >> L;
	
	cout << count_primitive_Pythagorean_triple(L / 4) << endl;

	return 0;
}

int main()
{
	//clock_t start, end;
	//start = clock();

	yuki0152();

	//end = clock();
	//printf("%d msec.\n", end - start);

	return 0;
}
0