[数位DP 莫比乌斯反演] 2017 计蒜之道 复赛 A. 阿里云秘钥池

大概可以推出来转移是
fx=d|xμ(d)P1di=1gid=d|xμ(d)G(d)
转移是 O(nlnn)
然后就直接数位dp咯

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

#define read(x) scanf("%d",&(x))

const ll PP=2e18;

inline void add(ll &x,ll y){
  x+=y; if (x>=PP) x-=PP;
}

const int maxn=100005;
int prime[maxn],vst[maxn];
int num,mu[maxn];

inline void Pre(int n){
  mu[1]=1;
  for (int i=2;i<=n;i++){
    if (!vst[i]) prime[++num]=i,mu[i]=-1;
    for (int j=1;j<=num && (ll)i*prime[j]<=n;j++){
      vst[i*prime[j]]=1;
      if (i%prime[j]==0){
    mu[i*prime[j]]=0; break;
      }else{
    mu[i*prime[j]]=-mu[i];
      }
    }
  }
}

int P;

int a[70],m;
ll f[70][maxn],g[70][maxn];

inline void calc(ll *f,ll *g){
  for (int i=0;i<P;i++) g[i]=0;
  for (int i=1;i<P;i++)
    for (int j=i;j<P;j+=i)
      add(g[i],f[j]);
}

inline ll Solve(ll n){
  m=0;
  while (n) a[++m]=n%P,n/=P;

  int flag=1;
  for (int i=0;i<P;i++) f[m][i]=0;
  for (int i=1;i<a[m];i++) add(f[m][i],1);
  calc(f[m],g[m]);

  for (int t=m-1;t;t--){
    for (int i=0;i<P;i++) f[t][i]=0;
    for (int i=1;i<P;i++)
      if (mu[i])
    for (int j=i;j<P;j+=i)
      add(f[t][j],mu[i]==1?g[t+1][i]:(PP-g[t+1][i])%PP);
    for (int i=1;i<P;i++)
      add(f[t][i],1);
    if (flag){
      for (int i=1;i<a[t];i++)
    if (__gcd(i,a[t+1])==1)
      add(f[t][i],1);
      if (a[t]==0 || __gcd(a[t],a[t+1])!=1)
    flag=0;
    }
    calc(f[t],g[t]);
  }

  ll ret=0;
  for (int i=1;i<P;i++)
    add(ret,f[1][i]);
  return ret;
}

int main(){
  int T; ll L,R;
  freopen("t.in","r",stdin);
  freopen("t.out","w",stdout);
  read(T); Pre(100000);
  while (T--){
    scanf("%lld%lld%d",&L,&R,&P);
    printf("%lld\n",Solve(R+1)-Solve(L));
  }
  return 0;
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值