2017-9-22 NOIP模拟赛[xxy][数论]

XXY 的 的 NOIP 模拟赛 4 4 —— 数学专场

A

Description
定义 f（x）表示 x 的约数和，例：f（12）=1+2+3+4+6+12=28
给出 x，y，求Σf（i），i∈[x,y]
Input
一行两个整数 x，y
Output
一行表示答案
Example
两组输入数据
2 4
123 321
对应输出
14
72543
Hint
对于 20%的数据，1<=x<=y<=1000
对于 40%的数据，1<=x<=y<=1e7
对于 100%的数据，1<=x<=y<=2e9

#include<iostream>
#include<cstdio>
using namespace std;
long long x,y;
long long count(int x){
    long long res=;
    long long l=,r;
    while(){
        r=x/(x/l);
        long long a=x/l;
        res+=((l+r)*(r-l+)/)*a;
        l=r+;
        if(x/l==)break;
    }
    return res;
}
int main(){
    freopen("A.in","r",stdin);
    freopen("A.out","w",stdout);
    cin>>x>>y;
    long long ansl=count(x-);
    long long ansr=count(y);
    cout<<ansr-ansl;
}

90分 出现整数被0除的错误

#include<iostream>
#include<cstdio>
using namespace std;
long long x,y;
long long count(int x){
    long long res=;
    long long l=,r;
    while(){
        if(x/l==)r=x;
        else r=x/(x/l);
        long long a=x/l;
        res+=((l+r)*(r-l+)/)*a;
        l=r+;
        if(x/l==)break;
    }
    return res;
}
int main(){
    freopen("A.in","r",stdin);
    freopen("A.out","w",stdout);
    cin>>x>>y;
    long long ansl=count(x-);
    long long ansr=count(y);
    cout<<ansr-ansl;
}

100分

B

Description
求满足以下条件的 x 的个数
① x∈[1,n!]
② 设 pi 为质数且 pi|x，那么 pi>m，对于所有的 pi 均成立
答案对 100000007 取模
Input
一行两个整数 n，m
Output
一行表示答案
Example
3 组输入样例
100 10
100 20
10000 9000
Output
对应 3 组输出
43274465
70342844
39714141
Hint
对于 20%的数据，n，m<=8
对于 60%的数据，n<=100000,m<=8
对于 100%的数据，n<=10000000，保证 n-m<=100000，m<=n

学到的一点东西

1.x的所有素因子大于m，则x与m!互素
2.已知phi[(i-1)!],递推地求phi[i!]
　　如果i是素数,phi[i!]=phi[(i-1)!]*(i-1)
　　如果i不是素数，phi[i!]=phi[(i-1)!]*i

3.phi[n!]=phi[m!]*(n!/m!)

#include<iostream>
#include<cstdio>
#include<algorithm>
using namespace std;
int p[],n,m,w=,ans,cnt;
bool vis[],q[];
void dfs(int pos,long long sum){
    if(sum>w)return;
    if(!q[sum]&&sum!=)ans++,q[sum]=;
    //cout<<sum<<' ';
    for(int i=pos;i<=cnt;i++){
        if(sum*p[i]>w)return;
        dfs(i,sum*p[i]);
    }
}
int main(){
    //freopen("Cola.txt","r",stdin);
    freopen("B.in","r",stdin);
    freopen("B.out","w",stdout);
    scanf("%d%d",&n,&m);
    for(int i=;i<=n;i++)w*=i;
    for(int i=;i<=w;i++){
        if(!vis[i])p[++cnt]=i;
        for(int j=;j<=cnt&&i*p[j]<=w;j++){
            vis[i*p[j]]=;
            if(i%p[j]==)break;
        }
    }
    int pos=upper_bound(p+,p+cnt+,m)-p;
    dfs(pos,);
    for(int i=pos;i<=cnt;i++){
        if(cnt>w)break;
        if(!q[p[i]])ans++;
    }
    cout<<ans;
}

20分 暴力

#include<cstdio>
#define N 10000001
#define mod 100000007
using namespace std;
int cnt,p[],phi[N];
long long phifac[N];
bool v[N];
int main()
{
    freopen("B.in","r",stdin);
    freopen("B.out","w",stdout);
    phi[]=;
    for(int i=;i<N;i++)
    {
        if(!v[i])
        {
            p[++cnt]=i;
            phi[i]=i-;
        }
        for(int j=;j<=cnt;j++)
        {
            if(i*p[j]>=N) break;
            v[i*p[j]]=true;
            if(i%p[j]) phi[i*p[j]]=phi[i]*(p[j]-);
            else
            {
                phi[i*p[j]]=phi[i]*p[j];
                break;
            }
        }
    }
    int n,m; long long ans;
    scanf("%d%d",&n,&m);
        ans=;
        phifac[]=;
        for(int i=;i<=m;i++)
        if(!v[i]) phifac[i]=phifac[i-]*(i-)%mod;
        else phifac[i]=phifac[i-]*i%mod;
        ans=phifac[m];
        for(int i=m+;i<=n;i++) ans=ans*i%mod;
        printf("%I64d\n",ans-);
}

100分

C

Description
令 f（k）=n 表示 有 n 种方式，可以把正整数 k 表示成几个素数的乘积的形式。
例 10=2*5=5*2，所以 f（10）=2
给出 n，求最小的 k
Input
一行一个整数 n
Output
一行表示答案
Example
四组输入样例
1
2
3
105
对应输出：
2
6
12
720
Hint
20%的数据，k<=20
另外 20%的数据，n<=20
另外 20%的数据，k<=1e6
100%的数据，k,n<=2^60

#include<iostream>
#include<cstdio>
using namespace std;
long long a[];
int k,n,bin[],p[],cnt;
bool vis[];
bool check(int x){
    int l=,pos=;
    while(x!=){
        while(x%p[pos]!=)pos++;
        l=l+;
        bin[l]=;
        while(x%p[pos]==){
            bin[l]++;
            x=x/p[pos];
        }
    }
    long long res=;
    int sum=;
    for(int i=;i<=l;i++)sum+=bin[i];
    res=a[sum];
    for(int i=;i<=l;i++)res/=a[bin[i]];
    if(res==n)return ;
    return ;
}
int main(){
    //freopen("Cola.txt","r",stdin);
    freopen("C.in","r",stdin);
    freopen("C.out","w",stdout);
    a[]=;
    for(int i=;i<=;i++)a[i]=a[i-]*i;
    scanf("%d",&n);
    for(int i=;i<=;i++){
        if(!vis[i])p[++cnt]=i;
        for(int j=;j<=cnt&&i*p[j]<=;j++){
            vis[i*p[j]]=;
            if(i%p[j]==)break;
        }
    }
    for(int i=;i<=;i++){//k的所有可能值
        if(check(i)){
            printf("%d",i);
            return ;
        }
    }
}

30分 暴力枚举答案

#include<cstdio>
#include<iostream>
using namespace std;
typedef  long long LL;
int p[]={,,,,,,,,,,,,,,,,,,,,};
LL ans,n;
LL C[][];
void solve(int num,int lim,LL tot,LL now,int last)
{
    if(now>ans) return;
    if(tot==n) { ans=now; return ; }
    if(tot>n || num>) return;
    LL t=;
    for(int i=;i<=lim;i++)
    {
        t*=p[num];
        if(now>=ans/t) return;
        solve(num+,i,tot*C[last+i][i],now*t,last+i);
    }
}
int main()
{
    freopen("C.in","r",stdin);
    freopen("C.out","w",stdout);
    C[][]=;
    for(int i=;i<;i++)
    {
        C[i][]=;
        for(int j=;j<=i;j++)
            C[i][j]=C[i-][j-]+C[i-][j];
    }
    scanf("%I64d",&n);
    if(n==)
    {
           printf("1 2\n");
           return ;
    }
    ans=(LL)<<;
    solve(,,,,);
    printf("%I64d\n",ans);
}

100分 搜索+剪枝