Description
给出一个长度为 n n 的序列,要求找满足 lcm(b1,b2,...,bn)=max(b1,b2,...,bn),1≤bi≤ai l c m ( b 1 , b 2 , . . . , b n ) = m a x ( b 1 , b 2 , . . . , b n ) , 1 ≤ b i ≤ a i 的 b b 序列个数
Input
第一行一整数表示序列长度,之后 n n 个正整数
Output
输出满足条件的 b b 序列个数,结果模
Sample Input
4
1 4 3 2
Sample Output
15
Solution
一个序列的最大公倍数等于其最大值说明这个序列所有元素均为最大值的因子,那么我们考虑枚举最大值为 m m ,那么每个必须是 m m 的因子,其某些必须等于 m m ,假设的因子从小到大拿出来是 c1,c2,...,cd c 1 , c 2 , . . . , c d ,对于第 i i 个数字,最多可以取到,假设 ck≤ai<ck+1 c k ≤ a i < c k + 1 ,那么从 1 1 ~中可以取到 k k 个不同的的因子,对 a a 序列排序,每次二分搜出介于之间的 ai a i 数量 cntk c n t k ,记大于等于 m m 的数量为 cnt c n t ,那么对于 ai<m a i < m 的位置,其取值方案就是 ∏k=1d−1kcntk ∏ k = 1 d − 1 k c n t k ,而对于这 cnt c n t 个不小于 m m 的位置,必须要保证其中至少有一个位置取到,故方案数为 mcnt−(m−1)cnt m c n t − ( m − 1 ) c n t ,即最大值为 m m 的合法序列有,枚举 m m 累加答案即可,时间复杂度
Code
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
typedef pair<int,int>P;
const int INF=0x3f3f3f3f,maxn=100005;
#define mod 1000000007
vector<int>f[maxn];
void init(int n=1e5)
{
for(int i=1;i<=n;i++)
for(int j=i;j<=n;j+=i)
f[j].push_back(i);
}
int Pow(int a,int b)
{
int ans=1;
while(b)
{
if(b&1)ans=(ll)ans*a%mod;
a=(ll)a*a%mod;
b>>=1;
}
return ans;
}
int n,a[maxn];
int main()
{
init();
while(~scanf("%d",&n))
{
int m=0;
for(int i=0;i<n;i++)scanf("%d",&a[i]),m=max(m,a[i]);
sort(a,a+n);
int ans=0;
for(int i=1;i<=m;i++)
{
int res=1,size=f[i].size();
for(int j=0;j<size-1;j++)
{
int t1=f[i][j],t2=f[i][j+1];
int num=lower_bound(a,a+n,t2)-lower_bound(a,a+n,t1);
res=(ll)res*Pow(j+1,num)%mod;
}
int num=lower_bound(a,a+n,m+1)-lower_bound(a,a+n,f[i][size-1]);
res=(ll)res*(Pow(size,num)-Pow(size-1,num)+mod)%mod;
ans+=res;
if(ans>=mod)ans-=mod;
}
printf("%d\n",ans);
}
return 0;
}