四川大学第二届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;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值