acwing数论

目录

1291 轻拍牛头

1294 樱花

1307 牡牛和牝牛

1308 方程的解

1309 车的放置

1310 数三角形

1312 序列统计


1291 轻拍牛头

今天是贝茜的生日,为了庆祝自己的生日,贝茜邀你来玩一个游戏.

贝茜让 N 头奶牛(编号 1 到 N)坐成一个圈。

除了 1 号与 N 号奶牛外,i 号奶牛与 i−1 号和 i+1 号奶牛相邻,N号奶牛与 1 号奶牛相邻。

农夫约翰用很多纸条装满了一个桶,每一张纸条中包含一个 1 到 1000000 之间的数字。

接着每一头奶牛 i从桶中取出一张纸条,纸条上的数字用 Ai 表示。

所有奶牛都选取完毕后,每头奶牛轮流走上一圈,当走到一头奶牛身旁时,如果自己手中的数字能够被该奶牛手中的数字整除,则拍打该牛的头。

牛们希望你帮助他们确定,每一头奶牛需要拍打的牛的数量。即共有 N 个整数 A1,A2,…,AN,对于每一个数 Ai,求其他的数中有多少个是它的约数。

输入格式

第一行包含整数 N

接下来 N 行,每行包含一个整数 Ai

输出格式

共 N 行,第 i 行的数字为第 i 头牛需要拍打的牛的数量。

数据范围

1≤N≤105
1≤Ai≤106

输入样例:

5
2
1
2
3
4

输出样例:

2
0
2
1
3

#include<bits/stdc++.h>
using namespace std;
int a[100010],cnt[1000010],s[1000010];
int main()
{
    int n;
    scanf("%d",&n);
    for(int i=0;i<n;i++){
        scanf("%d",&a[i]);
        cnt[a[i]]++;
    }
    for(int i=1;i<=1000000;i++){
        for(int j=i;j<=1000000;j+=i){
            s[j]+=cnt[i];
        }
    }
    for(int i=0;i<n;i++)printf("%d\n",s[a[i]]-1);
    return 0;
}

1294 樱花

给定一个整数 n,求有多少正整数数对 (x,y) 满足 1/x+1/y=1/n!。

输入格式

一个整数 n。

输出格式

一个整数,表示满足条件的数对数量。

答案对 109+7取模。

数据范围

1≤n≤106

输入样例:

2

输出样例:

3

样例解释

共有三个数对 (x,y) 满足条件,分别是 (3,6),(4,4),(6,3)

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e6+10,mod=1e9+7;
int p[N],cnt;
bool st[N];
void get(int n)
{
    for(int i=2;i<=n;i++){
        if(!st[i])p[cnt++]=i;
        for(int j=0;p[j]*i<=n;j++){
            st[p[j]*i]=true;
            if(i%p[j]==0)break;
        }
    }
}
int main()
{
    int n;
    cin>>n;
    get(n);
    LL res=1;
    for(int i=0;i<cnt;i++){
        int s=p[i];
        int g=0;
        for(int j=n;j;j/=s)g+=j/s;
        res=res*(2*g+1)%mod;
    }
    cout<<res;
    return 0;
}

1307 牡牛和牝牛

约翰要带 N 只牛去参加集会里的展示活动,这些牛可以是牡牛,也可以是牝牛。

牛们要站成一排,但是牡牛是好斗的,为了避免牡牛闹出乱子,约翰决定任意两只牡牛之间至少要有 K 只牝牛。

请计算一共有多少种排队的方法,所有牡牛可以看成是相同的,所有牝牛也一样,答案对 5000011 取模。

输入格式

一行,输入两个整数 N 和 K。

输出格式

一个整数,表示排队的方法数。

数据范围

1≤N≤105
0≤K<N

输入样例:

4 2

输出样例:

6

样例解释

6 种方法分别是:牝牝牝牝,牡牝牝牝,牝牡牝牝,牝牝牡牝,牝牝牝牡,牡牝牝牡。

#include<bits/stdc++.h>
using namespace std;
const int N=100010,mod=5000011;
int f[N],s[N];
int n,m;
int main()
{
    cin>>n>>m;
    f[0]=s[0]=1;
    for(int i=1;i<=n;i++){
        f[i]=s[max(0,i-m-1)];
        s[i]=(s[i-1]+f[i])%mod;
    }
    cout<<s[n];
    return 0;
}

1308 方程的解

佳佳碰到了一个难题,请你来帮忙解决。
对于不定方程a1+a2+…+ak-1+ak=g(x),其中k≥2且k∈N,x是正整数,g(x)=x^x mod 1000(即x^x除以1000的余数),x,k是给定的数。我们要求的是这个不定方程的正整数解组数。
举例来说,当k=3,x=2时,分别为(a1,a2,a3)=(2,1,1),(1,2,1),(1,1,2)。

输入格式

有且只有一行,为用空格隔开的两个正整数,依次为 k,x

有且只有一行,为方程的正整数解组数。

数据范围

1≤k≤100
1≤x<231
k≤g(x)

输入样例:

3 2

输出样例:

3

#include<bits/stdc++.h>
using namespace std;
const int N=150;
int k,x;
int f[1000][100][N];
int qmi(int a,int k,int p)
{
    int res=1;
    while(k)
    {
        if(k&1)res=res*a%p;
        a=a*a%p;
        k>>=1;
    }
    return res;
}
void add(int c[],int a[],int b[])
{
    for(int i=0,t=0;i<N;i++){
        t+=a[i]+b[i];
        c[i]=t%10;
        t/=10;
    }
}
int main()
{
    cin>>k>>x;
    int n=qmi(x%1000,x,1000);
    for(int i=0;i<n;i++){
        for(int j=0;j<=i&&j<k;j++){
            if(!j)f[i][j][0]=1;
            else add(f[i][j],f[i-1][j-1],f[i-1][j]);
        }
    }
    int *g=f[n-1][k-1];
    int i=N-1;
    while(!g[i])i--;
    while(i>=0)cout<<g[i--];
}

1309 车的放置

 

问题 AM: 车的放置

时间限制: 1.000 Sec  内存限制: 128 MB
提交 状态

题目描述

有下面这样的一个网格棋盘,a,b,c,d表示了对应边长度,也就是对应格子数。


当a=b=c=d=2时,对应下面这样一个棋盘,

要在这个棋盘上放K个相互不攻击的车,也就是这K个车没有两个车在同一行,也没有两个车在同一列,问有多少种方案。同样只需要输出答案mod 100003后的结果。

输入

5个非负整数a, b, c, d和k。

输出

1个正整数,为答案mod 100003后的结果。

样例输入 Copy

2 2 2 2 2

样例输出 Copy

38

提示

对于部分数据,有b = 0;
对于部分数据,有a,b,c,d≤4。
对于100%的数据,a,b,c,d,k≤1000,且保证了至少有一种可行方案。

 

 

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=2010,mod=100003;
int f[N],inf[N];
int qmi(int a,int k)
{
    int res=1;
    while(k)
    {
        if(k&1)res=(LL)res*a%mod;
        a=(LL)a*a%mod;
        k>>=1;
    }
    return res;
}
int C(int a,int b)
{
    if(a<b)return 0;
    return (LL)f[a]*inf[a-b]%mod*inf[b]%mod;
}
int p(int a,int b)
{
    if(a<b)return 0;
    return (LL)f[a]*inf[a-b]%mod;
}
int main()
{
    f[0]=inf[0]=1;
    for(int i=1;i<N;i++){
        f[i]=(LL)f[i-1]*i%mod;
        inf[i]=(LL)qmi(f[i],mod-2);
    }
    int a,b,c,d,k;
    cin>>a>>b>>c>>d>>k;
    int ans=0;
    for(int i=0;i<=k;i++){
        ans=(ans+(LL)C(b,i)*p(a,i)%mod*C(d,k-i)%mod*p(a+c-i,k-i))%mod;
    }
    cout<<ans;
    return 0;
}

1310 数三角形

问题 AN: 数三角形

时间限制: 1.000 Sec  内存限制: 128 MB
提交 状态

题目描述

给定一个nxm的网格,请计算三点都在格点上的三角形共有多少个。下图为4x4的网格上的一个三角形。

注意:三角形的三点不能共线。n×m的网格共有(n+1)×(m+1)个格点。

输入

输入一行,包含两个空格分隔的正整数m和n(1<=m,n<=1000)。

输出

输出一个正整数,为所求三角形数量。

样例输入 Copy

2 2

样例输出 Copy

76

 

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;

int gcd(int a,int b)
{
    return b?gcd(b,a%b):a;
}
LL C(int n)
{
    return (LL)n*(n-1)*(n-2)/6;
}
int main()
{
    LL n,m;
    cin>>n>>m;
    n++,m++;
    LL ans=C(n*m)-n*C(m)-m*C(n);
    for(int i=1;i<=n;i++){
        for(int j=1;j<=m;j++){
            ans-=2ll*(gcd(i,j)-1)*(n-i)*(m-j);
        }
    }
    cout<<ans;
    return 0;
}

1312 序列统计

问题 AP: 序列统计

时间限制: 1.000 Sec  内存限制: 128 MB
提交 状态

题目描述

给定三个正整数N、L和R,统计长度在1到N之间,元素大小都在L到R之间的单调不降序列的数量。输出答案对10^6+3取模的结果。

输入

输入第一行包含一个整数T,表示数据组数。
第2到第T+1行每行包含三个整数N、L和R,N、L和R的意义如题所述。
1≤N,L,R≤10^9,1≤T≤100,输入数据保证L≤R。

输出

输出包含T行,每行有一个数字,表示你所求出的答案对10^6+3取模的结果。

样例输入 Copy

2
1 4 5
2 4 5

样例输出 Copy

2
5

 

#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int p=1e6+3;
int qmi(int a,int k)
{
    int res=1;
    while(k)
    {
        if(k&1)res=(LL)res*a%p;
        a=(LL)a*a%p;
        k>>=1;
    }
    return res;
}
int c(int a,int b)
{
    if(a<b)return 0;
    int c1=1,c2=1;
    for(int i=a,j=1;j<=b;i--,j++){
        c1=(LL)c1*i%p;
        c2=(LL)c2*j%p;
    }
    return (LL)c1*qmi(c2,p-2)%p;
}
int L(int a,int b)
{
    if(a<p&&b<p)return c(a,b);
    return (LL)L(a/p,b/p)*c(a%p,b%p)%p;
}
int main()
{
    int t;
    cin>>t;
    while(t--)
    {
        int n,l,r;
        cin>>n>>l>>r;
        cout<<(L(r-l+n+1,r-l+1)+p-1)%p<<endl;
    }
    return 0;
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值