标签：cot int daimayuan aj 因子 计数 ai 4000 一题
思路
a
i
2
+
a
j
=
x
2
a_i^2+a_j=x^2
ai2+aj=x2
移项得：
a
j
=
x
2
−
a
i
2
a_j=x^2-a_i^2
aj=x2−ai2
平方差公式得:
a
j
=
(
x
+
a
i
)
(
x
−
a
i
)
a_j=(x+a_i)(x-a_i)
aj=(x+ai)(x−ai)
可以发现
a
j
a_j
aj 为其两个因子的乘积且满足
(
x
+
a
i
)
−
(
x
−
a
i
)
=
2
∗
a
i
(x+a_i) - (x-a_i) = 2*a_i
(x+ai)−(x−ai)=2∗ai
1
e
6
1e6
1e6范围内数因子个数最多量级为
4000
4000
4000
那么可以枚举因子
利用短除法求因子
O
(
n
∗
n
)
O(n*\sqrt{n})
O(n∗n
)可能会T，这里利用倍增
O
(
n
∗
l
o
g
(
n
)
)
O(n*log(n))
O(n∗log(n))求因子
for(int i = 1;i <= 1000000;i++)
for(int j = i;j <= 1000000;j+= i)
// 第一层枚举一个因子 i , j/i为第二个因子，乘积为j
code
int n;
int a[1000005];
int cot[1000005];
signed main()
{
_orz;
cin>>n;
forr(i,1,n) cin>>a[i],cot[a[i]]++;
int res = 0;
// 倍数取因子;
// 会出现两次 例如: i = 1 j = 3 // i = 3 j = 3
//
for(int i = 1;i <= 1000000;i++){
for(int j = i;j <= 1000000;j+= i){
int mi = min(j/i,i),ma = max(j/i,i);
int d = ma-mi;
if(d%2==0){
res += cot[j]*cot[d/2];
// if(cot[j] && cot[d/2])
// cout << j <<" "<< d/2 << endl;
}
}
}
cout << res/2 << '\n';
return 0;
}
标签：cot,int,daimayuan,aj,因子,计数,ai,4000,一题 来源： https://blog.csdn.net/qq_51687628/article/details/123649508
本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享； 2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关； 3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关； 4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除； 5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。