結果

問題 No.1274 楽しい格子点
コンテスト
ユーザー ok
提出日時 2020-10-31 00:16:42
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,309 bytes
コンパイル時間 6,372 ms
コンパイル使用メモリ 281,876 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-07-22 03:58:11
合計ジャッジ時間 6,218 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 38 WA * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<string>
#include<iomanip>
#include<cmath>
#include<vector>
#include<algorithm>
#include<utility>
#include <boost/multiprecision/cpp_dec_float.hpp>

using namespace std;
namespace mp = boost::multiprecision;

#define int long long
#define endl "\n"

constexpr long long INF = (long long)1e18;
constexpr long long MOD = 1'000'000'007; 

struct fast_io {
	fast_io(){
		std::cin.tie(nullptr);
		std::ios::sync_with_stdio(false);
	};
} fio;


long long gcd(long long a, long long b){return b?gcd(b,a%b):a;}

signed main(){
	cout<<fixed<<setprecision(15);
	
	int A, B;
	// long double ans = 0;
	mp::cpp_dec_float_100 ans = 0;
    int GCD = 0;
    
    cin>>A>>B;
    
    if(A == 0 && B == 0) {
        cout<<0.5<<endl;
        return 0;
    }
    
    A = abs(A);
    B = abs(B);
    
    GCD = gcd(A, B);
    
    if(GCD == 1) {
        cout<<0.337189<<endl;
        return 0;
    }
    
    // for(int i = 1; i <= 500; i += GCD){
        // for(int j = 1; j <= 500; j += GCD){
            // mp::cpp_dec_float_100 x = i + j;
            // ans += mp::pow(x,-x);
        // }
    // }
    
    for(int i = 1; i <= 1000; i += GCD){
         mp::cpp_dec_float_100 x = i+1, y = (i + GCD - 1)/GCD;
         
         ans += mp::pow(x,-x) * y;
    }
    
    cout<<ans<<endl;
	
	return 0;
}
0