题目
问题来自于一个精彩的故事:
有三个数学家,A,B与C。A选了两个正整数x与y满足x<=y。然后,A将x+y的值告诉了B,A又将x*y的值告诉了C。B与C都不知道x与y分别是什么,也不知道对方得到的值是什么。但B和C知道A告诉B的值是某两个正整数的“和”而告诉C的值是这两个数的“积”。而且这三个数学家的数学功底足够好。下面是B与C进行的对话:
B:“我确定你一定没有百分百的把握猜中我得到的数。”
C:“谢谢你的提示。现在我能确定你获得的数是 S。”
故事结束,回到问题。
这个故事中一共涉及3个未知参数x,y与S,其实由于S=x+y,所以实际一共只有两个未知参数而已。你可以带入一些正整数让这个故事没有逻辑漏洞。现在问题来了,在区间[L,R]上存在多少个数值t,使S=t时能找到对应的x与y,并让这个故事成立。输出这些t的和(1<=L<=R<=5,000,000)。
Input
多组测试数据,第一行一个整数T,表示数据个数,其中1<=T<=10.
之后有T行,每行两个数L与R,表示一组询问,其中1<=L<=R<=5,000,000.
Output
输出T行,每行一个整数,即[L,R]区间中所有符合条件的S的和。
Input示例
3
30 33
8 11
40 43
Output示例
33
19
0
分析
首先,怎么知道对方无法100%猜中?
考虑对方什么时候100%猜中。显然是拿了个质数。
因此和S一定不可以被表示成质数+1,不然对方100%猜中。
B说的那句话也就是告诉了C这样的信息“我的数不是质数+1”
那么为什么C就猜出来了呢?肯定有唯一解!
假设C拿到的是P, P为合数。
对于1*P=P这组拆分,P+1显然不是质数+1(P是合数嘛)
唯一解就肯定是1和P,那么其余所有解xy=P,都会有x+y-1是质数。
因此知道了P=S-1。
枚举x和y用n log n筛就行了。
代码
#include <bits/stdc++.h>
#define N 5000010
long long prime[N / 5];
long long sum[N];
bool notPrime[N];
void getSum()
{
int cnt = 0;
for (int i = 2; i <= N - 10; i++)
{
if (!notPrime[i])
prime[++cnt] = i;
for (int j = 1; j <= cnt && i <= (N - 10) / prime[j]; j++)
{
notPrime[i * prime[j]] = 1;
if (i % prime[j] == 0)
break;
}
}
for (int i = 1; i <= N - 10; i++)
sum[i] = i;
for (int i = 2; i <= N - 10; i++)
for (int j = i; j <= (N - 10) / i; j++)
if (notPrime[i + j - 1])
sum[i * j + 1] = 0;
for (int i = 1; i <= N - 10; i++)
if (!notPrime[i])
sum[i + 1] = 0;
sum[1] = sum[2] = 0;
for (int i = 1; i <= N - 10; i++)
sum[i] += sum[i - 1];
}
int main()
{
getSum();
int T;
scanf("%d",&T);
while (T--)
{
int l,r;
scanf("%d%d",&l,&r);
printf("%lld\n",sum[r] - sum[l - 1]);
}
}