关闭

[中等题] Project Euler 608 Divisor Sums

97人阅读 评论(0) 收藏 举报
分类:

这个题怎么Difficulty rating 80% 啊,送经验的感觉啊

D(m,n)=======d|mk=1nσ0(kd)d|mk=1na|kb|d[(a,b)=1]a=1nnad|mb|d[(a,b)=1]a=1nnad|mb|di|a,i|bμ(i)i|mμ(i)×(σ01)(mi)idnnidi|mμ(i)×(σ01)(mi)d=1ninidi|mμ(i)×(σ01)(mi)×f(ni)

这个 σ01 就是 111 ,就是把一个数分成三个有序的数的乘积的方案数。
满足条件且 μ(i)0i 不会太多,大概是 4.4×107 个,直接搜索。
f 这个东西全部预处理好,总计算量我算算了也就 109 级别,还是能接受的。

于是跑了 2.5min 答案就出来了

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

const ll n=1e12,m=200;
const int P=1e9+7;

int p[205],num;
int vst[205];

int q[205];

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

inline ll C(ll x){
  return (x+2)*(x+1)/2;
}

ll Map1[1000005],Map2[1000005];

inline ll Calc(ll x){
  ll ret=0;
  for (ll i=1,j,t;i<=x;i=j+1){
    j=x/(t=x/i);
    ret+=t*(j-i+1);
  }
  return (x<=1e6?Map1[x]:Map2[n/x])=ret%P;
}
inline ll Map(ll x){
  return x<=1e6?Map1[x]:Map2[n/x];
}

ll ans=0;

inline void dfs(ll t,ll cur,ll mu,ll phi){
  static int Cnt=0;
  if (t==num+1){
    if ((++Cnt)%1000000==0)
      printf("%d\n",Cnt);
    if (mu==1)
      ans+=phi*Map(n/cur)%P;
    else
      ans+=P-phi*Map(n/cur)%P;
    return;
  }
  if (cur*p[t]<=n)
    dfs(t+1,cur*p[t],-mu,phi*C(q[t]-1)%P);
  dfs(t+1,cur,mu,phi*C(q[t])%P);
}

inline ll Brute(){
  ll mm=1;
  for (int i=1;i<=m;i++) mm*=i;
  int ret=0;
  for (int i=1;i<=mm;i++)
    if (mm%i==0)
      for (int j=1;j<=n;j++)
    for (int k=1;k<=i*j;k++)
      if ((i*j)%k==0)
        ret++;
  printf("%d\n",ret);
  return ret;
}
inline void GG(){
  //ll ret=0;
  ll Cnt=0;
  for (ll i=1,j;i<=n;i=j+1){
    j=n/(n/i);
    Calc(n/i);
    if ((++Cnt)%10000==0)
      printf("%lld %lld\n",i,Cnt);
  }
  //printf("%lld\n",ret);
  //return ret;
}

int main(){
  Pre(m);
  for (int i=1;i<=num;i++){
    q[i]=0;
    for (ll t=p[i];t<=m;t*=p[i])
      q[i]+=m/t;
  }
  GG();
  dfs(1,1,1,1);
  printf("%lld\n",ans%P);
  //Brute();
  return 0;
}
0
0
查看评论

最小费用最大流

最小费用最大流就是指网络最大流的集合中找出最小费用的那个最大流。 其基本思想是从某个可行流F出发,找到关于这个流的一个可改进路经P,然后沿着P调整F,对新的可行流试图寻找关于他的可改进路经,如此反复直至求得最大流。现在要找最小费用的最大流,可以证明,若F是流量为V(F)的流中费用最小者,而P是关于...
  • sunmenggmail
  • sunmenggmail
  • 2012-03-31 00:01
  • 564

ProjectEuler题解(更新到100题)

欧拉项目题解,源代码从Github上下载:https://github.com/cloudzfy/euler Multiples of 3 and 5Even Fibonacci numbersLargest prime factorLargest palindro...
  • cloudzfy1
  • cloudzfy1
  • 2016-12-14 12:10
  • 536

Project Euler 1-5题

第1题 题目来源ProjectEuler这个题求的是严格小于1000的数中,是3或5的倍数的数的和。(刚开始理解错below的意思了,把1000算进去了,尴尬)int main(){ int ans=0; for (int i=1;i<1000;i++){ ...
  • pfccWang
  • pfccWang
  • 2017-09-10 21:25
  • 190

Project Euler一句话题解(此贴持续更新)

好好搞一波数学,从PE第三页开刷101:Lagrange插值多项式 Ans:37076114526
  • MrBird_to_fly
  • MrBird_to_fly
  • 2016-10-20 20:45
  • 305

Project-Euler problem 1-50

最近闲的做了下Project Euler 上的题目,前面50题都比较简单,简单总结下,一下代码一般是Python和C/C++的 用Python 做这些题目简直是酸爽啊 一下代码可能不一定是我的,因为不知道论坛里面的回复不是永久的,所以我的代码有的丢了,可能找个和我的意思相近的代码。题目翻译是从 欧...
  • u011401504
  • u011401504
  • 2015-01-31 00:09
  • 2028

【Project Euler】3 第三题

 //The prime factors of 13195 are 5, 7, 13 and 29.         //What is the largest prime factor of the number 600...
  • NoMasp
  • NoMasp
  • 2015-02-07 12:21
  • 1107

Project Euler 10

本文章来自我的个人网站,如感兴趣,欢迎访问我的个人网站:http://www.qingshuimonk.com/ 今天终于搞完了PE的第十题,说起来这道题做了得有小半个月了。 算法类似于PE的第十题,还是用之前的质数去除,以此来判断是否是质数,所以代码就在第七题的基础上稍稍改动一下就可以了。...
  • michaelxi007
  • michaelxi007
  • 2013-11-17 23:16
  • 786

Project Euler 31-35题

第31题 题目来源ProjectEuler这一题求200可以由1,2,5,10,20,50,100,200以多少种不同的方式相加而成。 将硬币分别标记为f[1]=1;f[2]=2;f[3]=5;f[4]=10;f[5]=20;f[6]=50;f[7]=100];f[8]=200; 使用dp[i...
  • pfccWang
  • pfccWang
  • 2017-09-20 13:02
  • 125

[反演] Project Euler 608. Divisor Sums

D(m,n)=&#x2211;d|m&#x2211;k=1n&#x03C3;0(kd)" role="presentation">D(m,n)=∑d|m∑k=1nσ0(kd)D(m,n)=∑d|m∑k=1nσ0(kd)D(m,n)=\dis...
  • Coldef
  • Coldef
  • 2018-01-29 07:47
  • 55

Project Euler Problem 51-70

题目翻译是从 欧拉计划 | Project Euler 中文翻译站上面Copy 的 Problem 51 Prime digit replacements 通过置换*3的第一位得到的9个数中,有六个是质数:13,23,43,53,73和83。 通过用同样的数字...
  • u011401504
  • u011401504
  • 2015-02-08 02:22
  • 1685
    个人资料
    • 访问:364253次
    • 积分:13016
    • 等级:
    • 排名:第1238名
    • 原创:975篇
    • 转载:3篇
    • 译文:0篇
    • 评论:56条
    文章分类
    最新评论