結果

問題 No.2979 直角三角形の個数
ユーザー V_MelvilleV_Melville
提出日時 2024-12-03 16:43:06
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,488 bytes
コンパイル時間 805 ms
コンパイル使用メモリ 88,276 KB
実行使用メモリ 31,616 KB
最終ジャッジ日時 2024-12-03 16:43:52
合計ジャッジ時間 46,019 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 89 ms
31,616 KB
testcase_01 AC 89 ms
21,120 KB
testcase_02 AC 90 ms
31,488 KB
testcase_03 AC 89 ms
21,120 KB
testcase_04 AC 88 ms
21,120 KB
testcase_05 AC 91 ms
31,488 KB
testcase_06 AC 92 ms
21,248 KB
testcase_07 AC 93 ms
24,360 KB
testcase_08 AC 130 ms
31,616 KB
testcase_09 AC 103 ms
17,668 KB
testcase_10 AC 759 ms
26,544 KB
testcase_11 AC 174 ms
15,872 KB
testcase_12 TLE -
testcase_13 AC 1,780 ms
15,872 KB
testcase_14 TLE -
testcase_15 TLE -
testcase_16 TLE -
testcase_17 TLE -
testcase_18 TLE -
testcase_19 WA -
testcase_20 WA -
testcase_21 TLE -
testcase_22 WA -
testcase_23 WA -
testcase_24 AC 83 ms
16,000 KB
testcase_25 AC 86 ms
16,000 KB
testcase_26 AC 83 ms
16,000 KB
testcase_27 TLE -
testcase_28 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <string.h>
#include <stdio.h>
#include <math.h>
 
using namespace std;
typedef long long LL;
 
const int N=10000005;
const int M=1005;
 
bool prime[N];
int p[N];
int pri[M];
int num[M];
int k,cnt,tmp,tmp1,count;
 
void isprime()
{
    k=0;
    int i,j;
    memset(prime,true,sizeof(prime));
    for(i=2;i<N;i++)
    {
        if(prime[i])
        {
            p[k++]=i;
            for(j=i+i;j<N;j+=i)
            {
                prime[j]=false;
            }
        }
    }
}
 
void Find(int n)
{
    cnt=0;
    int t=(int)sqrt(n*1.0),i,a;
    for(i=0;p[i]<=t;i++)
    {
        if(n%p[i]==0)
        {
            a=0;
            pri[cnt]=p[i];
            while(n%p[i]==0)
            {
                a++;
                n/=p[i];
            }
            num[cnt]=2*a;
            cnt++;
        }
    }
    if(n>1)
    {
        pri[cnt]=n;
        num[cnt]=2;
        cnt++;
    }
}
 
void dfs(int dept, LL product=1)
{
    if(dept==cnt)
    {
        if(product%2==0&&product>tmp&&product<tmp1)
           count++;
        return;
    }
    for(int i=0;i<=num[dept];i++)
    {
        dfs(dept+1,product);
        product*=pri[dept];
    }
}
 
int main()
{
    int n;
    scanf("%d",&n);
    
    isprime();
    
    int ans = 0;
    for (int i = 1; i <= n; ++i) {
        count=0;
        tmp=(LL)(sqrt(2.0)*i);
        tmp1=2*i;
        Find(i);
        dfs(0,1);
        ans += count;
    }
    
    cout << ans << '\n';
    
    return 0;
}
0