传送门唉我觉得这题数据范围1e5都能做啊。。。居然只出了2000
考完听zxyzxyzxy说我的贪心可以卡但过了?
可能今天本来是0+10+00+10+00+10+0只是运气好T1T1T1骗了100pts100pts100pts吧233.
下面讲讲贪心:
令wi=∏i=1kpiaiw_i=\prod_{i=1}^kp_i^{a_i}wi=∏i=1kpiai
令A=∑aipi是好质数,B=∑ai,pi是坏质数A=\sum_{a_i} p_i是好质数,B=\sum_{a_i},p_i是坏质数A=∑aipi是好质数,B=∑ai,pi是坏质数
那么fi=A−Bf_i=A-Bfi=A−B
令gi=gcd(w1,w2,...,wi)g_i=gcd(w_1,w_2,...,w_i)gi=gcd(w1,w2,...,wi)
考虑j<ij<ij<i
有gi∣gjg_i|g_jgi∣gj
令k=gjgik=\frac {g_j} {g_i}k=gigj
那么有fgj=fk+fgif_{g_j}=f_k+f_{g_i}fgj=fk+fgi
所以如果先在iii那个地方进行一次操作,无论怎么样都不会令结果变小。
于是我们从后向前算出所有有正的贡献的iii加到答案里面就行了。
代码:
#include<bits/stdc++.h>
#include<tr1/unordered_map>
using namespace std;
const int N=32005;
int n,m,pri[N],a[5005],b[5005],tot=0,sum=0,ans=0;
bool vis[N];
tr1::unordered_map<int,int>S;
inline int read(){
int ans=0;
char ch=getchar();
while(!isdigit(ch))ch=getchar();
while(isdigit(ch))ans=(ans<<3)+(ans<<1)+(ch^48),ch=getchar();
return ans;
}
inline void init(){
for(int i=2;i<=N-5;++i){
if(!vis[i])pri[++tot]=i;
for(int j=1;j<=tot&&i*pri[j]<=N-5;++j){
int k=i*pri[j];
vis[k]=1;
if(i%pri[j]==0)break;
}
}
}
inline int calc(int x){
int ret=0;
for(int i=1;i<=tot&&pri[i]*pri[i]<=x;++i){
if(x==1)break;
if(x!=x/pri[i]*pri[i])continue;
int cnt=0;
while(x==x/pri[i]*pri[i])++cnt,x/=pri[i];
if(S[pri[i]])ret-=cnt;
else ret+=cnt;
}
if(x^1){
if(S[x])--ret;
else ++ret;
}
return ret;
}
int main(){
init();
n=read(),m=read(),S.clear();
for(int i=1;i<=n;++i)a[i]=read();
for(int i=1;i<=m;++i)S[read()]=1;
for(int i=1;i<=n;++i){
sum+=calc(a[i]);
if(i^1)a[i]=__gcd(a[i],a[i-1]);
b[i]=-calc(a[i]);
}
int pos=n,delta=0;
while(pos){
if(b[pos]-delta>0)sum+=(b[pos]-delta)*pos,delta+=b[pos]-delta;
--pos;
}
cout<<sum;
return 0;
}