四川大学第二届SCUACM新生赛（同步赛）E题 [模板]欧拉筛

题目链接： [模板]欧拉筛
题目描述：
lglg最近爱上了炒股，他想在股票市场内买一些股票。

现在有编号为1，2，3，…，100000000的100000000支股票，每个股票的价值等于其编号的阶乘(例如编号为5的股票的价值就是120)。
lglg是一个很挑剔的人，他只喜欢编号为质数的股票，但他很有钱，因此他希望买下所有编号小于等于N(1<=N<=1e8)并且编号为质数的股票，请你帮他算一算这些股票的价值之和。

由于价值和很大，他希望你能告诉他价值之和对P(1e3<=P<=1e5)取模后的值
输入描述:

输入包含多组数据 第一行是一个正整数 T(1<=T<=1000)，代表数据的组数。
每组数据只有一行，包含两个正整数 N 和 P 数字之间用空格隔开，
其代表的含义及范围已在题面中给出

输出描述:

每组数据输出一行，表示lglg希望买下的所有股票的价值之和对P取模后的值

示例1
输入

2
5 1001 
20 1001  

输出

128
86 

代码：

#include<iostream>
#include<math.h>
using namespace std;
typedef long long ll;
const ll MAX = 1e6 + 5;//100000000;
int prime[MAX];
ll v[MAX];
void init_prime()
{
	int i, j;
	for(i = 2;i <= sqrt(1000002.0); ++i)
	{
		if(!prime[i])
			for(j = i * i; j < 1000002; j += i)
				prime[j] = 1;
	}
	j = 0;
	for(i = 2;i <= 1000002; ++i)
		if(!prime[i]) 
			prime[j++] = i;
}

int main () {
    int t;
    ll n, p;
    scanf ("%d", &t);
    init_prime();
    while (t--) {
        scanf ("%lld %lld", &n, &p);
        n = min (n, p);
        v[1] = 1;
        ll sum = 0;
        for (ll i = 2; i <= n; i++) {
            v[i] = (i * v[i - 1]) % p;
        }
        for (ll i = 0; prime[i] <= n; i++) {
            sum = (sum  + v[prime[i]]) % p;
        }
        printf("%lld\n",sum);
    }
    return 0;
}