【jxoi2018】游戏 组合数学

首先令$n=r-l+1$。

令$k$表示区间$[l,r]$中存在多少个数$x$，使得$x$不存在小于$x$且在区间$[l,r]$中的因数，我们把包含这些数的数集称为$S$

我们来先想一个$O(nk)$的$min-max$容斥做法吧。。。。。

显然这一题我们可以转化为min-max容斥的模型（将这k个数选完期望需要选多少次）

$max_{S}=\sum_{T∈S}(-1)^{|T+1|}min_{T}$。

令$P_x=\sum_{T∈S\ and\ |T|=x} min_{T}$。

我们推一推式子就会发现$P_i=x!(n-x)!\sum_{i=1}{n-k+1}i\binom{n-i}{k-i}$。

然后我们发现这个式子是$O(n^2)$的，而且非常难以推出。

代码如下（这个代码可能有点假）

 #include<bits/stdc++.h>
 #define L long long
 #define MOD 1000000007
 #define M 10000005
 using namespace std;

 L pow_mod(L x,L k){L ans=; for(;k;k>>=,x=x*x%MOD) if(k&) ans=ans*x%MOD; return ans;}
 L fac[M]={},invfac[M]={};
 L C(int n,int m){return fac[n]*invfac[m]%MOD*invfac[n-m]%MOD;}

 int vis[M]={};
 int n,k=;

 L p[M]={};

 int main(){
     fac[]=; for(int i=;i<M;i++) fac[i]=fac[i-]*i%MOD;
     invfac[M-]=pow_mod(fac[M-],MOD-);
     for(int i=M-;~i;i--) invfac[i]=invfac[i+]*(i+)%MOD;

     int l,r; cin>>l>>r; n=r-l+;
     for(int i=l;i<=r;i++){
         if(vis[i]) continue;
         k++;
         for(int j=i;j<=r;j+=i) vis[j]=;
     }

     for(int x=;x<=k;x++){
         L now=;
         for(int i=;i<=n-x+;i++){
             L s=i;
             for(int j=;j<x;j++) s=s*(n-i-j+)%MOD;
             now=(now+s)%MOD;
         }
         p[x]=now*x%MOD*fac[n-x]%MOD;
     }
     L ans=;
     for(L x=,zf=;x<=k;x++,zf=-zf){
         ans=(ans+zf*p[x]*C(k,x)%MOD+MOD)%MOD;
     }
     cout<<ans<<endl;
 }

考虑一些简单的方法

我们考虑回题目中的枚举排列。令$F_i$表示 $t(p)=i$的排列个数，那么答案显然为$\sum_{i=k}^{n}F_i$

不难发现，一种$t(p)=i$的排列，其前$i-1$项中必包含有数集$S$中$k-1$个数，且第i个数必为数集$S$中的数。

那么不难求出$F_i=k(n-k)!\dfrac{i!}{(i-k)!}$

答案即为$k(n-k)!\sum_{i=k}^{n} \dfrac{i!}{(i-k)!}$

随便求一求就好了

 #include<bits/stdc++.h>
 #define L long long
 #define MOD 1000000007
 #define M 10000005
 using namespace std;

 L pow_mod(L x,L k){L ans=; for(;k;k>>=,x=x*x%MOD) if(k&) ans=ans*x%MOD; return ans;}
 L fac[M]={},invfac[M]={};
 L C(int n,int m){return fac[n]*invfac[m]%MOD*invfac[n-m]%MOD;}

 int vis[M]={};
 int n,k=;

 L p[M]={};

 int main(){
     fac[]=; for(int i=;i<M;i++) fac[i]=fac[i-]*i%MOD;
     invfac[M-]=pow_mod(fac[M-],MOD-);
     for(int i=M-;~i;i--) invfac[i]=invfac[i+]*(i+)%MOD;

     int l,r; cin>>l>>r; n=r-l+;
     for(int i=l;i<=r;i++){
         if(vis[i]) continue;
         k++;
         for(int j=i;j<=r;j+=i) vis[j]=;
     }
     L ans=k*fac[n-k]%MOD,sum=;
     for(int i=k;i<=n;i++)
     sum=(sum+fac[i]*invfac[i-k])%MOD;
     cout<<ans*sum%MOD<<endl;
 }