题目描述
俗话说,有钱就是任性。我们的高富帅鱼丸同学打算去看电影。鱼丸到了电影院以后,发现座位的编号正好是 11 到 200200 。但是有一些座位号对应的座位坏掉了,没法坐,不妨假设还剩下 NN 个能坐的椅子。电影的老板告诉鱼丸,如果你要包下一个集合 SS 里的所有椅子,就要付出这些椅子的编号的最小公倍数的钱。鱼丸很任性地同意了。
来这里玩了很多天以后,鱼丸发现自己正好来了 2^N-12
N
−1 天,并且由于他非常任性,对于这 NN 个椅子的每一种可能的非空子集,他都包下过来看电影。鱼丸大少爷虽然不在乎花了多少钱,但你毕竟是他的助理,于是你想知道鱼丸一共花了多少钱。由于钱的数量实在太大,请对答案 mod~1e9+7mod 1e9+7之后输出。
输入格式
第一行输入一个数 N(1 \le N \le 200)N(1≤N≤200),表示能做的椅子的数量。
接下来一行 NN 个 11 到 200200 之间的整数,用空格隔开。
输出格式
一行输出答案。
样例
样例输入
2
2 4
样例输出
10
数据范围与提示
对于 20%20% 的数据,有 N \le 20N≤20
对于 100%100% 的数据,有 N \le 200N≤200
显示分类标签
题解:
最先对于开多维数组我是拒绝的,可它太香了
但是根据最后计算答案时的方法,可以发现对于大于13的质数新设一维和直接乘进方案数是没有区别的。。。。
对于大于13的质数单独考虑,小于的直接枚举,暴力转移的新方法。
#include<bits/stdc++.h>
using namespace std;
const long long mod=1e9+7;
// 2 3 5 7 11 13
long long f[8][5][4][3][3][3][2];
long long g[8][5][4][3][3][3][2];
long long n;
long long tot,prim[205],v[205];
inline void shai(long long x){
for(long long i=2;i<=x;++i){
if(!v[i]){
v[i]=i;prim[++tot]=i;
}
for(long long j=1;j<=tot;++j){
if(v[i]<prim[j]||i*prim[j]>x)break;
v[i*prim[j]]=prim[j];
}
}
//for(long long i=1;i<=tot;++i)cout<<prim[i]<<endl;
}
inline long long power(long long x,long long c){
long long now=1;
x=x%mod;
while(c){
if(c&1)now=1LL*now*x%mod;
x=1LL*x*x%mod;c>>=1;
}
return now;
}
struct node{
long long x,up;
}a[205];
bool cmp(node a,node b){return a.up<b.up;}
long long fe[10];
inline void fen(long long x){
for(long long i=1;i<=6;++i){
fe[i]=0;
while(x%prim[i]==0){
fe[i]++;x=x/prim[i];
}
}
for(long long i=7;i<=tot;++i){
if(x%prim[i]==0)fe[7]=i;
}
}
void test(long long i1,long long i2,long long i3,long long i4,long long i5,long long i6){
if(f[i1][i2][i3][i4][i5][i6][0])
cout<<i1<<" "<<i2<<" "<<i3<<" "<<i4<<" "<<i5<<" "<<i6<<" "<<f[i1][i2][i3][i4][i5][i6][0]<<endl;
}
void Test(){
for(long long i1=7;i1>=0;--i1)
for(long long i2=4;i2>=0;--i2)
for(long long i3=3;i3>=0;--i3)
for(long long i4=2;i4>=0;--i4)
for(long long i5=2;i5>=0;--i5)
for(long long i6=2;i6>=0;--i6)
test(i1,i2,i3,i4,i5,i6);
}
int main(){
scanf("%lld",&n);
shai(200);
for(long long i=1;i<=n;++i){
scanf("%lld",&a[i].x);
for(long long j=tot;j;--j){
if(a[i].x%prim[j]==0){
a[i].up=j;
break;
}
}
}
sort(a+1,a+1+n,cmp);
//cout<<1<<endl;
long long las=1;
f[0][0][0][0][0][0][0]=0;
for(long long i=1;i<=n;++i){
//cout<<a[i].up<<endl;
if(a[i].up>6)break;
las=i+1;
memset(g,0,sizeof(g));
//cout<<i<<endl;
fen(a[i].x);
for(long long i1=7;i1>=0;--i1)
for(long long i2=4;i2>=0;--i2)
for(long long i3=3;i3>=0;--i3)
for(long long i4=2;i4>=0;--i4)
for(long long i5=2;i5>=0;--i5)
for(long long i6=2;i6>=0;--i6)
f[max(fe[1],i1)][max(fe[2],i2)][max(fe[3],i3)][max(fe[4],i4)][max(fe[5],i5)][max(fe[6],i6)][0]+=f[i1][i2][i3][i4][i5][i6][0],
f[max(fe[1],i1)][max(fe[2],i2)][max(fe[3],i3)][max(fe[4],i4)][max(fe[5],i5)][max(fe[6],i6)][0]%=mod;
f[fe[1]][fe[2]][fe[3]][fe[4]][fe[5]][fe[6]][0]+=1;
f[fe[1]][fe[2]][fe[3]][fe[4]][fe[5]][fe[6]][0]%=mod;
//cout<<i<<endl;
}
//cout<<1<<endl;
//cout<<las<<endl;
// Test();
for(;las<=n;){
for(long long i=las;i<=n;++i){
las=i+1;
//memset(g,0,sizeof(g));
fen(a[i].x);
for(long long i1=7;i1>=0;--i1)
for(long long i2=4;i2>=0;--i2)
for(long long i3=3;i3>=0;--i3)
for(long long i4=2;i4>=0;--i4)
for(long long i5=2;i5>=0;--i5)
for(long long i6=2;i6>=0;--i6)
f[max(fe[1],i1)][max(fe[2],i2)][max(fe[3],i3)][max(fe[4],i4)][max(fe[5],i5)][max(fe[6],i6)][1]+=
1LL*f[i1][i2][i3][i4][i5][i6][0]*prim[a[i].up]%mod+f[i1][i2][i3][i4][i5][i6][1],
f[max(fe[1],i1)][max(fe[2],i2)][max(fe[3],i3)][max(fe[4],i4)][max(fe[5],i5)][max(fe[6],i6)][1]%=mod;
f[fe[1]][fe[2]][fe[3]][fe[4]][fe[5]][fe[6]][1]+=prim[a[i].up];
f[fe[1]][fe[2]][fe[3]][fe[4]][fe[5]][fe[6]][1]%=mod;
if(i<n&&a[i+1].up!=a[i].up)break;
}
for(long long i1=7;i1>=0;--i1)
for(long long i2=4;i2>=0;--i2)
for(long long i3=3;i3>=0;--i3)
for(long long i4=2;i4>=0;--i4)
for(long long i5=2;i5>=0;--i5)
for(long long i6=2;i6>=0;--i6)
f[i1][i2][i3][i4][i5][i6][0]=(f[i1][i2][i3][i4][i5][i6][1]+f[i1][i2][i3][i4][i5][i6][0])%mod,
f[i1][i2][i3][i4][i5][i6][1]=0;
//Test();
//cout<<endl;
}
long long ans=0;
for(long long i1=7;i1>=0;--i1)
for(long long i2=4;i2>=0;--i2)
for(long long i3=3;i3>=0;--i3)
for(long long i4=2;i4>=0;--i4)
for(long long i5=2;i5>=0;--i5)
for(long long i6=2;i6>=0;--i6)
ans=(ans+1LL*power(2,i1)*power(3,i2)%mod*power(5,i3)*power(7,i4)%mod*power(11,i5)*power(13,i6)%mod*f[i1][i2][i3][i4][i5][i6][0]%mod)%mod;
//Test();
printf("%lld\n",(ans+mod)%mod);
return 0;
}
/*
11
195 149 193 181 177 164 27 109 142 62 171
*/