洛谷—数论/简单数学问题

SuperGCD

解题思路： 其实很简单c++选手可能会当做一道高精度去做，而高精度相对就会比较麻烦，处理细节
后续： 但是参考了其他大牛，有人用更相减损法，很简短的代码

链接： 代码很长，可以参考 传送门

计算系数

解题思路： 先理解多项式展开的原理，分析题型是一个dp,不断更新项的系数，之后就是确定dp公式
代码：

    #include<iostream>
    #include<cstdio>
    #include<algorithm>
    #include<cmath>
    #include<cstring>
    using namespace std;
    int a,b,k,n,m;
    long long f[1005][1005];
    int main()
    {
        scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);
        f[0][0]=1;
        for (int i=0;i<=n;i++)
            for (int j=0;j<=m;j++)
            {
                if (i==0 && j==0) continue;
                f[i][j]=0;
                if (i>0)
                    f[i][j]=(f[i][j]+f[i-1][j]*a)%10007;
                if (j>0)
                    f[i][j]=(f[i][j]+f[i][j-1]*b)%10007;
            }
        printf("%lld\n",f[n][m]);
        return 0;
    }

扫雷

解题思路： 看到网上大神都在用dp,但是一开始我也是想到深度搜索，将可能性遍及一次，
代码： 学习一下dp，原理大概是确定第一个然后再去确定后面的

#include<bits/stdc++.h>
using namespace std;
int i,j,k,m,n,ans=2,a[10001],b[10001];
void checkx()
{
    for(int i=2;i<=n+1;i++)
    {
    b[i]=a[i-1]-b[i-1]-b[i-2];
    if (!(b[i]==1||b[i]==0))
          {
          ans--;
          break;
          }
          if (i==n+1&&b[i]!=0)
                {
                 ans--;
                 break;
             }
    }
}
int main()
{
    ios::sync_with_stdio(false);
    cin>>n;
    b[0]=0;
    for (i=1;i<=n;i++)    
    cin>>a[i];
    b[1]=1;
    checkx();
    b[1]=0;
     checkx();
    cout<<ans;   
return 0;    
}

余数求和

解题思路： 简单看会以为是暴力，但是其实是有一定转化的思想，取模的公式不只是a%b,还可以表示为a-b*(a/b),利用这个去化简步骤，不会超时
代码：

#include<cstdio>
#include<algorithm>
using namespace std;
typedef long long ll;

int main() {
    ll n,k;
    scanf("%lld%lld",&n,&k);
    ll ans=n*k;
    for(ll l=1,r;l<=n;l=r+1) {
        if(k/l!=0) r=min(k/(k/l),n); 
        else r=n;
        ans-=(k/l)*(r-l+1)*(l+r)/2;
    }
    printf("%lld",ans);
    return 0;
}

2^k进制数

解题思路： 都是关于排列组合的问题，虽然知道题型，但是就是怎么也写不出公式，根据题意，从所有n个数中选m个数后，只有一种组合方式，去计算m的等式即可
**代码：**查看大佬言简意赅的代码

#include<iostream>
#include<cstdio>
#define rint register int
#define ini inline int 
using namespace std;
int a[30001][201],c[201],ans[201],k,w;
ini jia(int *a,int *b)
{
    int lenc=0,x=0;
    while (lenc<a[0] || lenc<b[0])
    {
        a[++lenc]=a[lenc]+b[lenc]+x;
        x=a[lenc]/10;
        a[lenc]%=10;    
    }
    if (x>0) a[++lenc]=x;
    a[0]=lenc;
}
int main()
{
    scanf("%d%d",&k,&w);
    int kk=w%k;int hh=w/k;
    int y=0;
    for (rint i=1;i<=kk;i++)
        y+=1<<(i-1);     //最高位最大是几
    int minn=(1<<k)-1;//第一位最大是多少
    if (hh==1 || (hh==2 && kk==0))
    {
        if (kk==0) y=minn;int tot=0;
        for (rint i=1;i<=y;i++)
            tot+=minn-i;
        cout<<tot;
        return 0;
    }//特判两位以内的情况
    for (rint i=1;i<=minn-1;i++)
    {
        a[i][1]=i;a[i][0]=1;
        jia(ans,a[i]);
    }//第二位是minn-1到1各有多少种可能
    for (rint i=3;i<=hh;i++)
        for (rint j=1;j<=minn-i+1;j++)
//第i位最大是minn-i+1
//第i位是1到minn-i+1各有多少种可能
        {
            jia(a[j],a[j-1]);
            jia(ans,a[j]);
        }
//以下是单独处理最高位，因为最高位有限制
    for (rint j=1;j<=minn-hh;j++)
        jia(a[j],a[j-1]);
    for (rint j=minn-hh;j>=max(minn-hh-y+1,1);j--)
        jia(ans,a[j]);
//我们的前缀和是倒着存的
    for (rint i=ans[0];i>=1;i--)
        printf("%d",ans[i]);
}