#2402. 任性（willful）

题目描述
俗话说，有钱就是任性。我们的高富帅鱼丸同学打算去看电影。鱼丸到了电影院以后，发现座位的编号正好是 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

*/