洛谷P4562 [JXOI2018]游戏数论

最新推荐文章于 2025-02-05 17:22:58 发布

「已注销」

最新推荐文章于 2025-02-05 17:22:58 发布

阅读量173

点赞数

CC 4.0 BY-SA版权

本文链接：https://blog.csdn.net/qq_42692055/article/details/108634313

本文探讨了如何通过筛选法找出特定区间内无法由该区间数的因数表示的所有数，并进一步计算这些数在特定组合中的贡献。使用了埃氏筛法进行筛选，详细介绍了筛选过程及如何计算最终答案。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

正解:数论

解题报告:

传送门!

首先考虑怎么样的数可能出现在$t(i)$那个位置上?显然是$[l,r]$中所有无法被表示出来的数(就约数不在$[l,r]$内的数嘛$QwQ$

所以可以先把这些数筛出来

具体怎么筛的话,最原始的方法就埃氏筛?

然后显然可以线性筛,但我$jio$得大概快不到哪儿去而且麻烦一些懒得打了所以直接用的埃氏筛

然后现在就筛出来,有$sum$个满足条件的数了,考虑怎么算贡献?(先注明下,,,可能有些±1的细节下面都直接略过了$QAQ$

就先枚最后一个这样的数出现的位置$x$(也就是$t(i)$的取值

首先它自己会有$x$的贡献

然后在它后面可以从$n-sum$个数中选出$n-x$个数,而且这些数可以随便排列,所以有$C(n-x,n-sum)\cdot (n-x)!$

然后在它之前的数可以随便排列,所以有$(x-1)!$

最后这个数是从$sum$个数中选出来的所以还有$C(sum,1)$

综上,可以列出式子,$ans=\sum i\cdot sum\cdot (i-1)!\cdot (n-x)!\cdot C(n-x,n-sum)$

然后这题好像还有优化,,,但是我懒得落实了(bushi

反正一时半会儿应该不会再写优化了,,,如果很久很久之后考古出了这篇博客也许会填下QAQ?(四舍五入就是咕了:D

#include<bits/stdc++.h>
using namespace std;
#define il inline
#define ll long long
#define gc getchar()
#define ri register int
#define rc register char
#define rb register bool
#define rp(i,x,y) for(ri i=x;i<=y;++i)
#define my(i,x,y) for(ri i=x;i>=y;--i)
 
const int mod=1000000007,N=10000000+10;
int l,r,sum,fac[N],ifac[N];
bool vis[N];
 
 
il int read()
{
    rc ch=gc;ri x=0;rb y=1;
    while(ch!='-' && (ch>'9' || ch<'0'))ch=gc;
    if(ch=='-')ch=gc,y=0;
    while(ch>='0' && ch<='9')x=(x<<1)+(x<<3)+(ch^'0'),ch=gc;
    return y?x:-x;
}
il int power(ri x,ri y){int ret=1;while(y){if(y&1)ret=1ll*ret*x%mod;x=1ll*x*x%mod;y>>=1;}return ret;}
il void pre(ri n)
{
    fac[0]=1;rp(i,1,n)fac[i]=1ll*fac[i-1]*i%mod;
    ifac[n]=power(fac[n],mod-2);my(i,n-1,0)ifac[i]=1ll*ifac[i+1]*(i+1)%mod;
}
inline int C(ri n,ri m){return 1ll*fac[n]*ifac[m]%mod*ifac[n-m]%mod;}
 
int main()
{
//  freopen("4562.in","r",stdin);//freopen("4562.out","w",stdout);
    l=read();r=read();pre(r);
    rp(i,l,r)if(!vis[i]){++sum;for(ri j=i;j<=r;j+=i)vis[j]=1;}
    ri as=0,n=r-l+1;
    rp(i,sum,n)as=(as+1ll*i*sum%mod*C(n-sum,n-i)%mod*fac[n-i]%mod*fac[i-1]%mod)%mod;
    printf("%d\n",as);
    return 0;
}

最后放下代码鸭!