Description
对任意正整数$n$，$f(n)$为满足$|x^2-y^2|=n$的有序整数对$(x,y)$的个数。
现在给出$m$个正整数$a_1,a_2,\cdots,a_m$
令$sum=f(a_1)+f(a_2)+\cdots +f(a_m)-m$（保证$sum\gt 1$）
若$sum$为合数，求出在$[1,sum]$区间内所有满足以下条件的合数$x$的个数：
$$x|(2^x-1)$$
若$sum$为质数，求出在$[1,sum^2]$区间内所有满足以下条件的整数x的个数：
$$sum|(2^x-1)$$
Input Format
第一行给出了正整数$m$
接下来$m$行为$m$个正整数$a_1,a_2,\cdots ,a_m$
Output Format
第一行输出$sum$。
第二行输出满足条件的$x$的个数。
Sample Input
1 4
Sample Output
3 4
Hint
$(x,y)=(2,0),(-2,0),(0,2),(0,-2)$
故$sum=4-1=3$
$3$为质数
满足条件的$x$分别为：$2,4,6,8$，共$4$个
【数据规模与约定】
$40\%$：$m\le 1000$
$40\%$：$sum$为合数
$100\%$：$m\le 10^4，a_i\le 10^{10}$
解题报告
$f(x)$的值
根据题意，$f(n)$为满足$|x^2-y^2|=n$的有序数对$(x,y)$的个数。
将式子稍微转变一下，变成$|(x+y)(x-y)|=n$，所以设$a=|x+y|$，$b=|x-y|$。所以我们用$[a,b]$来代替$(x,y)$，则$[a,b]$需要满足$a \cdot b = n$，且$a$、$b$的奇偶性相同。
当$x、y \neq 0$时，由$(x,y)$衍生出的$(x,y)$、$(x,-y)$、$(-x,y)$、$(-x,-y)$、$(y,x)$、$(y,-x)$、$(-y,x)$、$(-y,-x)$，都是满足题意的，将它们写成$[a,b]$形式，只能写出两个，即$[a,b]$、$[b,a]$。
当$x=0$或$y=0$时，在$(x,0)$与$(0,y)$中，必然有$x=y$，所以由$(x,y)$可以衍生出$(0,x)$、$(x,0)$、$(0,-x)$、$(-x,0)$，而写成$[a,b]$的形式只能写出一个。
综上所述，存在满足条件的$(x,y)$序对的个数为$[a,b]$序对个数的$4$倍。
对$n$的取值进行分类讨论，当$n$为奇数时，$n$只能被拆分乘两个奇数的乘积。设$n = \prod_{i=1}^k P_i^{t_i}$（$P_i$为互不相等的奇质数），则通过乘法原理易得将$n$拆分成两个奇数的乘积的方案数为$\prod_{i=1}^k (t_i+1)$。
当$n$为偶数时，$n$只能被拆分成两个偶数的乘积。设$n = 2^a \cdot \prod_{i=1}^k P_i^{t_i}$，则得到将$n$拆分成两个偶数的乘积的方案数为$(a - 1) \cdot \prod_{i=1}^k (t_i+1)$。
最后将结果乘以$4$即可。
$Sum$为合数
当$Sum$为合数时，要求求出区间$[1,Sum]$内满足$x|(2^x-1)$的合数$x$的个数。
若$x$为偶数，由于$2^x-1$必为奇数，则不存在$x$满足$x|(2^x-1)$。
若$x$为奇数，则假设$x|(2^x-1)$。这个式子等价于：$2^x \equiv 1\pmod{x}$。
假设已经找出一个最小的$r$满足：$2^r \equiv 1\pmod{x}$，则必然满足$r|x$。
设$P$为$x$的一个质因子，则必然有$2^r \equiv 1\pmod{P}$。且$2^{P-1} \equiv 1\pmod{P}$，由此可得$r|(P-1)$。
当$P$为$x$的最小质因子时（由于$x$为奇数，所以$P \gt 2$），若存在$r$满足$r|(P-1)$，则$P$不为$x$的最小质因子，矛盾。所以不存在奇合数$x$满足$x|(2^x-1)$。
综上所述，当Sum为合数是，区间$[1,Sum]$内满足$x|(2^x-1)$的合数$x$的个数为0。
$Sum$为质数
当$Sum$为质数时，要求求出区间$[1,sum^2]$内满足$Sum|(2^x-1)$的整数$x$的个数。
同样的，这个式子等价于：$2^x \equiv 1\pmod{Sum}$。
找出一个最小的$r$满足：$2^r \equiv 1\pmod{Sum}$，则必然满足$r|x$。即任何一个满足条件的$x$都是$r$的倍数。
因为$2^r\equiv1\pmod{Sum}$，所以$2^{kr}=(2^r)^k\equiv1^k=1\pmod{Sum}$。（$k\in Z$）。即所以任何一个$r$的倍数都满足条件。
综上所述，在区间$[1,Sum^2]$内满足条件的$x$的个数即在该区间内$r$的倍数的个数。
代码：
#include <cstdio>
int m, f[100010], p[10000], pn;
long long a, sum;
long long F(long long x)
{
long long ret = 1, t = 0;
if (x % 4 == 2) return 0;
if (x % 4 == 0) x /= 4;
for (int i = 0, j, k; i < pn && x > 1; i++)
{
j = p[i];
k = 1;
while (x % j == 0)
{
k++;
x /= j;
}
ret *= k;
}
if (x != 1) ret *= 2;
return ret;
}
bool isPrime(long long x)
{
if (x < 100000)
{
return f[x] == 0;
}
for (int i = 0; i < pn; i++)
{
if (x % p[i] == 0)
{
return false;
}
}
return true;
}
long long pow(long long x, long long p, long long mod)
{
if (p == 0) return 1;
if (p == 1) return x % mod;
if (p == 2) return x * x % mod;
return pow(pow(x, p / 2, mod), 2, mod) * pow(x, p % 2, mod) % mod;
}
long long Find_r(long long x)
{
long long ret = -1;
x--;
for (long long i = 2; i * i <= x; i++)
{
if (x % i) continue;
if (pow(2, i, x + 1) == 1) return i;
if (pow(2, x / i, x + 1) == 1)
{
if (ret == -1 || x / i < ret)
{
ret = x / i;
}
}
}
return ret;
}
int main()
{
f[0] = f[1] = 1;
for (int i = 2; i <= 100000; i++)
{
if (f[i] == 0)
{
for (int j = i + i; j <= 100000; j += i)
{
f[j] = i;
}
p[pn++] = i;
}
}
scanf("%d", &m);
for (int i = 0; i < m; i++)
{
scanf("%lld", &a);
sum += F(a) * 4 - 1;
}
printf("%lld\n", sum);
if (isPrime(sum))
{
long long r = Find_r(sum);
if (r == -1) //找不到r时即r=phi(sum)=sum-1
{
r = sum - 1;
}
sum = sum * sum / r;
printf("%lld", sum);
}
else
{
printf("0");
}
}