DP，数论————洛谷P4317 花神的数论题（求1~n二进制中1的个数和）

玄学代码（是洛谷题解里的一位dalao小粉兔写的）

//数位DP(二进制)计算出f[i]为恰好有i个的方案数。
//答案为∏(i^f[i]),快速幂解决。
#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int mod=;
ll n,ans=;
ll c,f[];
ll qpow(ll b,ll e){
    ll a=;
    for(;e;b=b*b%mod,e>>=)
        e&?a=a*b%mod:;
    return a;
}

int main(){
    scanf("%lld",&n);
    for(int j=;~j;--j){
        for(int i=;i;--i)
            f[i]+=f[i-];
        if(n>>j&)++f[c++];
    }++f[c];
    for(int i=;i<=;++i)
        ans=ans*qpow(i,f[i])%mod;
    printf("%lld",ans);
    return ;
}

下面是另一位dalaoStyx写的简单易懂的题解

代码如下

#include<cstdio>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#define mod 10000007
using namespace std;

long long ans=,n,k,cnt,dp[][];

long long kasumi(long long a,long long b)
{
    long long tmp=;
    while(b)
    {
        if(b&) tmp=(tmp%mod)*a%mod;
        a=(a%mod)*a%mod;
        b>>=;
    }
    return tmp%mod;
}

int main()
{
    scanf("%lld",&n);
    for(int i=; i<=; i++)
    {
        dp[i][]=;
    }
    for(int i=; i<=; i++)
    {
        dp[i][i]=;
    }
    for(int i=; i<=; i++)
    {
        for(int j=; j<i; j++)
        {
            dp[i][j]=dp[i-][j-]+dp[i-][j];
        }
    }
    for(k=; k<=n; k<<=)
    {
        n-=k;
        cnt++;
        for(int j=; j<=cnt; j++)
        {
            ans=(ans%mod)*kasumi(j,dp[cnt][j])%mod;
            ans%=mod;
        }
    }
    cnt=;
    for(int i=;i>=;i--)
    {
        if((1ll<<i)<=n)
        {
            for(int j=;j<=i+;j++)
            {
                ans=(ans%mod)*kasumi(j+cnt,dp[i+][j]);
                ans%=mod;
            }
            cnt++;
            n-=(1ll<<i);
        }
    }
    printf("%lld\n",ans);
}

至于为什么不用扩展欧拉定理降幂，那是因为杨辉三角在总和等于1e15的时候 ，每一个数肯定不会爆longlong