Luogu P4562 [JXOI2018]游戏

题目
我们用埃氏筛从(l,r)筛一遍，每次把没有被筛掉的数的倍数筛掉。
易知最后剩下来的数（这个集合记为(S)）的个数就是我们需要选的数，设有(s)个，令(n=r-l+1)。
记(f_i)为第(i)次操作筛完的方案数，我们要求的就是(sumlimits_{i=s}^n if_i)
枚举最后一次操作，方案数为(s)。
后面的位置从任选不属于(S)的数的方案数为(n-schoose n-i)。
前面(i-1)个数是哪些就确定了，排列方案数为((i-1)!)。
后面的同理((n-i)!)。
所以(f_i=s{n-schoose n-i}(i-1)!(n-i)!)。
答案(ans=sumlimits_{i=s}^n if_i=ssumlimits_{i=s}^n{n-schoose n-i}i!(n-i)!=s(n-s)!sumlimits_{i=s}^nfrac{i!}{(i-s)!})。

#include<bits/stdc++.h>
using namespace std;
const int N=10000007,P=1000000007;
bitset<N>t;int fac[N],ifac[N];
int inc(int a,int b){a+=b;return a>=P? a-P:a;}
int mul(int a,int b){return 1ll*a*b%P;}
int power(int a,int k=P-2){int r=1;for(;k;k>>=1,a=mul(a,a))if(k&1)r=mul(a,r);return r;}
int main()
{
    int i,j,l,r,s=0,ans=0,n;scanf("%d%d",&l,&r),n=r-l+1;
    for(i=l;i<=r;++i) if(!t[i]) for(++s,j=i;j<=r;j+=i) t[j]=1;
    for(fac[0]=i=1;i<=r;++i) fac[i]=mul(i,fac[i-1]);
    for(ifac[r]=power(fac[r]),i=r;i;--i) ifac[i-1]=mul(i,ifac[i]);
    for(i=s;i<=n;++i) ans=inc(ans,mul(fac[i],ifac[i-s]));
    return !printf("%d",mul(mul(s,ans),fac[n-s]));
}