結果

問題 No.2979 直角三角形の個数
ユーザー kotatsugamekotatsugame
提出日時 2024-12-25 06:27:34
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 1,144 bytes
コンパイル時間 987 ms
コンパイル使用メモリ 105,900 KB
実行使用メモリ 7,680 KB
最終ジャッジ日時 2024-12-25 06:27:51
合計ジャッジ時間 14,536 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 11 ms
7,680 KB
testcase_01 AC 11 ms
7,552 KB
testcase_02 AC 11 ms
7,424 KB
testcase_03 AC 13 ms
7,552 KB
testcase_04 AC 12 ms
7,564 KB
testcase_05 AC 11 ms
7,612 KB
testcase_06 AC 11 ms
7,560 KB
testcase_07 AC 11 ms
7,424 KB
testcase_08 AC 11 ms
7,680 KB
testcase_09 AC 11 ms
7,648 KB
testcase_10 AC 11 ms
7,680 KB
testcase_11 AC 11 ms
7,600 KB
testcase_12 AC 11 ms
7,484 KB
testcase_13 AC 11 ms
7,680 KB
testcase_14 AC 14 ms
7,552 KB
testcase_15 AC 15 ms
7,680 KB
testcase_16 AC 23 ms
7,608 KB
testcase_17 AC 46 ms
7,552 KB
testcase_18 AC 110 ms
7,552 KB
testcase_19 AC 149 ms
7,572 KB
testcase_20 AC 406 ms
7,552 KB
testcase_21 AC 622 ms
7,632 KB
testcase_22 AC 3,365 ms
7,680 KB
testcase_23 AC 3,318 ms
7,552 KB
testcase_24 AC 11 ms
7,552 KB
testcase_25 AC 11 ms
7,676 KB
testcase_26 AC 10 ms
7,540 KB
testcase_27 AC 12 ms
7,468 KB
testcase_28 TLE -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<cmath>
#include<vector>
#include<cassert>
using namespace std;
long N;
int pf[1<<20];
int ps[7],T[1<<7];
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(nullptr);
	for(int i=2;i<1<<20;i++)if(pf[i]==0)
	{
		for(int j=i;j<1<<20;j+=i)pf[j]=i;
	}
	cin>>N;
	long ans=0;
	for(int m=2;;m++)
	{
		long up=(double)N/(2*m);
		int l=m+1,r=(int)min(up,2*m-1L);
		if(l>r)break;
		int pssz=0;
		{
			int v=m;
			while(v%2==0)v/=2;
			ps[pssz++]=2;
			while(v>1)
			{
				int p=pf[v];
				ps[pssz++]=p;
				while(v%p==0)v/=p;
			}
		}
		assert(pssz<=7);
		long sum=0;
		T[0]=1;
		for(int i=0;i<1<<pssz;i++)
		{
			if(i)
			{
				int j=0;
				while(!(i>>j&1))j++;
				T[i]=T[i^1<<j]*ps[j];
			}
			int t=T[i];
			int nl=(double)(l+t-1)/t,nr=(double)r/t;
			long nup=(double)up/t;
			long cur=0;
			while(nl<=nr)
			{
				long q=(double)nup/nl;
				int nnl=nr;
				if(nup<q*nr)nnl=(double)nup/q;
				//int nnl=(int)min((long)nr,(long)((double)nup/q))+1;
				cur+=q*(nnl-nl+1);
				nl=nnl+1;
			}
			//for(int j=nl;j<=nr;j++)cur+=nup/j;
			if(__builtin_parity(i))sum-=cur;
			else sum+=cur;
		}
		ans+=sum;
	}
	cout<<ans<<endl;
}
0