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

 链接：https://ac.nowcoder.com/acm/contest/1838/E
来源：牛客网
 

时间限制：C/C++ 2秒，其他语言4秒
空间限制：C/C++ 131072K，其他语言262144K
64bit IO Format: %lld

题目描述

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 

我醉了 

((a%p)*(b%p))%p就TM会TLE

(a*b)%p就A了

我再也不 （k）

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N=1e6+30;
ll pri[N];
bool vis[N];
ll tot;
ll a[N];
ll n,p,tmp;

void prime()
{
    tot = 0;
    memset(vis, 0, sizeof(vis));
    for (ll i = 2; i <N; i++)
    {
        if (!vis[i])
        {
            pri[++tot] = i;
        }
        for (ll j = 1; j<=tot &&i*pri[j]<N; j++)
        {
            vis[i*pri[j]] = 1;
            if (i % pri[j] == 0)
                break;
        }
    }
}

void solve()
{
    memset(a,0,sizeof(a));
    a[0]=1;
    for(ll i=1; i<=tmp; i++)
    {
        a[i]=(a[i-1]*i)%p;
    }
}

int main()
{
    prime();
    solve();
    ll t;
    scanf("%lld",&t);
    while(t--)
    {
        scanf("%lld%lld",&n,&p);
        if(n>p)
            n=p-1;
        ll id=lower_bound(pri,pri+tot,n)-pri;
        tmp=pri[id];
        if(tmp>n)
            tmp=pri[id-1];
        solve();
        ll ans=0;
        for(ll i=1; i<=id; i++)
        {
            if(pri[i]>tmp)
                break;
            ans=(ans+a[pri[i]])%p;
        }
        cout<<ans<<'\n';
    }
    return 0;
}