#1033176 (C++14) No.2979 直角三角形の個数

• 問題: No.2979 直角三角形の個数
• ユーザー: kotatsugame
• 提出日時: 2024-12-25 06:44:00
• 言語: C++14 (gcc 13.3.0 + boost 1.87.0)
• 結果: AC
• 実行時間: 1,972 ms / 4,000 ms
• コード長: 1,043 bytes
• コンパイル時間: 442 ms
• コンパイル使用メモリ: 66,232 KB
• 実行使用メモリ: 7,552 KB
• 最終ジャッジ日時: 2024-12-25 06:44:10
• 合計ジャッジ時間: 7,455 ms
• ジャッジサーバーID （参考情報）: judge3 / judge5

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	10 ms; 7,424 KB
testcase_01	AC	10 ms; 7,472 KB
testcase_02	AC	12 ms; 7,404 KB
testcase_03	AC	12 ms; 7,424 KB
testcase_04	AC	12 ms; 7,500 KB
testcase_05	AC	11 ms; 7,448 KB
testcase_06	AC	10 ms; 7,424 KB
testcase_07	AC	10 ms; 7,552 KB
testcase_08	AC	12 ms; 7,424 KB
testcase_09	AC	11 ms; 7,504 KB
testcase_10	AC	11 ms; 7,552 KB
testcase_11	AC	10 ms; 7,424 KB
testcase_12	AC	10 ms; 7,424 KB
testcase_13	AC	11 ms; 7,424 KB
testcase_14	AC	12 ms; 7,552 KB
testcase_15	AC	12 ms; 7,424 KB
testcase_16	AC	17 ms; 7,424 KB
testcase_17	AC	28 ms; 7,520 KB
testcase_18	AC	82 ms; 7,424 KB
testcase_19	AC	77 ms; 7,464 KB
testcase_20	AC	199 ms; 7,424 KB
testcase_21	AC	315 ms; 7,552 KB
testcase_22	AC	1,537 ms; 7,424 KB
testcase_23	AC	1,575 ms; 7,552 KB
testcase_24	AC	10 ms; 7,552 KB
testcase_25	AC	11 ms; 7,444 KB
testcase_26	AC	12 ms; 7,548 KB
testcase_27	AC	12 ms; 7,424 KB
testcase_28	AC	1,972 ms; 7,508 KB

ソースコード

#include<iostream>
#include<cassert>
using namespace std;
long N;
int pf[1<<20];
int ps[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;
			while(v>1)
			{
				int p=pf[v];
				ps[pssz++]=p;
				while(v%p==0)v/=p;
			}
		}
		assert(pssz<=6);
		long sum=0;
		for(int i=0;i<1<<pssz;i++)
		{
			int t=1;
			for(int j=0;j<pssz;j++)if(i>>j&1)t*=ps[j];
			int nl=(double)(l+t-1)/t,nr=(double)r/t;
			if(nl%2==0)nl++;
			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;
				if(nnl%2==0)nnl--;
				//int nnl=(int)min((long)nr,(long)((double)nup/q))+1;
				cur+=q*((nnl-nl)/2+1);
				nl=nnl+2;
			}
			if(__builtin_parity(i))sum-=cur;
			else sum+=cur;
		}
		ans+=sum;
	}
	cout<<ans<<endl;
}