代码源每日一题Div.1 平方计数(枚举因子,埃氏筛)

 题目链接

题意

 

思路

     设a_{i}^{2} + a_{j} = x^{2},则a_{j} = x^{2}-a_{i}^{2}=\left ( x-a_{i} \right )*(x+a_{i})

     分别设f_1=\left ( x-a_{i} \right )f_2=\left ( x+a_{i} \right )

      则 f_1,f_2 均为 a_{j} 的的因数,且 f_2-f_1=2*a_i,可枚举每一个 a_{j} 的因数,确定 a_i 的值,若暴力枚举,复杂度 O(n\sqrt{n}),会T;因为对于每一个数 x,我们枚举了  1-\sqrt{x} 每一个数,但其中有很多不是 x 的因数,怎样才能只枚举因数,而不枚举其他数;

     我们对于每一个因数 x,枚举其能够作为谁的因数,答案显然是 x 的倍数都可以,至此问题转变成了枚举每一个数的倍数,再思考一下枚举倍数的过程,埃氏筛!!!事先保存每一个数出现的次数,通过埃氏筛枚举倍数,复杂度O(nlogn)

代码

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;

int a[1000010], tong[1000010];

int main()
{
	ios::sync_with_stdio(false); cin.tie(0); cout.tie(0);
	int n;
	cin >> n;
	for(int i = 1; i <= n; i++){
		cin >> a[i]; tong[a[i]]++;
	}
	
	ll ans = 0;
	for(int i = 1; i <= 1000000; i++){
		for(int j = i; j <= 1000000; j += i){
			int mx = max(j / i, i), mn = min(j / i, i); //表示两个因子 
			int tmp = mx - mn; 
			if(tmp % 2) continue; //判断是否存在这样的ai 
			ans += 1ll * tong[j] * tong[tmp / 2]; //j表示aj,tmp / 2表示ai 
		}
	}
	cout << ans / 2 << "\n"; //一对因子会计算两次,如2*5=10,5*2=10 
	return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值