话不多说,直接看题:
1.
我们可以得到大致一个思路,就是先枚举1-1e6的质数,然后看看有几个即可。
我们怎么知道个数呢?
首先我们知道1---n中有n/p的下取整个为p的倍数。
因此,p的个数至少是n/p的下取整个,当然有些数有不止1个p的倍数,于是我们得到n/p^2+n/p^3+...直到p^i>n.
下面是AC代码:
#include<bits/stdc++.h>
using namespace std;
const int N=1000010;
int pri[N],cnt;
bool st[N];
void getpri(int n){
for(int i=2;i<=n;i++){
if(!st[i]) pri[cnt++]=i;
for(int j=0;pri[j]*i<=n;j++){
st[pri[j]*i]=1;
if(i%pri[j]==0) break;
}
}
}
int main(){
int n;
cin>>n;
getpri(n);
for(int i=0;i<cnt;i++){
int p=pri[i];
int s=0,t=n;
while(t){
s+=t/p;
t/=p;
}
printf("%d %d\n",p,s);
}
}
2.
我们把数按照基本算数定理拆分,我们对于每一个质数的奇数,答案就乘这个,若为偶数,那么就不用管,下面是AC代码:
#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
int main(){
LL n;
cin>>n;
LL res=1;
for(LL i=2;i*i<=n;i++){
if(n%i==0){
int s=0;
while(n%i==0){
s++;
n/=i;
}
if(s%2) res*=i;
}
}
if(n>1) res*=n;
cout<<res;
}
3.
下面是分析:
下面是AC代码:
#include<bits/stdc++.h>
using namespace std;
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
int T,a0,a1,b0,b1;
int ans;
int main(){
cin>>T;
while(T--){
ans=0;
scanf("%d%d%d%d",&a0,&a1,&b0,&b1);
int p1=a0/a1,p2=b1/b0;
for(int i=1;i*i<=b1;i++){
if(b1%i) continue;
if(i%a1==0&&gcd(i/a1,p1)==1&&gcd(b1/i,p2)==1) ans++;
if(i==b1/i) continue;
int ck=b1/i;
if(ck%a1==0&&gcd(ck/a1,p1)==1&&gcd(b1/ck,p2)==1) ans++;
}
printf("%d\n",ans);
}
}
4.
这里主要介绍一个定理:求两个数的约数等价于最大公约数的因子,于是我们求出gcd,然后求约数即可,下面是AC代码:
#include<bits/stdc++.h>
using namespace std;
int q[5000],cnt;
int gcd(int a,int b){
return b?gcd(b,a%b):a;
}
int main(){
int a,b;
cin>>a>>b;
int d=gcd(a,b);
for(int i=1;i<=d/i;i++){
if(d%i==0){
q[cnt++]=i;
if(i!=d/i) q[cnt++]=d/i;
}
}
sort(q,q+cnt);
int n;
cin>>n;
while(n--){
int L,R;
scanf("%d%d",&L,&R);
int f=0;
for(int i=cnt-1;i>=0;i--){
if(q[i]>=L&&q[i]<=R){
printf("%d\n",q[i]);
f=1;
break;
}
}
if(f==0) cout<<-1<<endl;
}
}