欧拉工程第70题:Totient permutation
和上面几题差不多的
Euler's Totient function, φ(n) [sometimes called the phi function]:小于等于n的数并且和n是互质的数的个数
存在这样的数:N的欧拉数φ(n),是N的一个排列
例如:φ(87109)=79180
求在1---10^7中n/φ(n) 取到最小的 n 是多少?
这里的是p是n的素因子,当素因子有相同的时候只取一个
任意一个正整数都能分解成若干个素数乘积的形式
直接利用上题的phi函数就可以求解
这个是跑的最快的函数了
long phi2(long n){
long res = 0;
if(n==1) return 0;
int pi=2;
res = n;
while(pi*pi <=n){
if(n%pi==0){
res/=pi;
res*=(pi-1);
while(n%pi==0){
n/=pi;
}
}
pi++;
}
if(n>1){
res/=n;
res*=(n-1);
}
return res;
}
// 8319823
// running time=89s43ms
运行结果也90s了
然后考虑到
这里的pi都是素数
所以我们可以只考虑素数的情况
当n%pi==0的时候这个pi就符合条件,并去重
void run2(){
long Max_n = 10000000;
int Prime_len=400;
int[] Prime = new int[Prime_len];
Prime[0] = 2;
int k=1;
int p=1;
while(k<Prime_len){
if(isPrime(p)){
Prime[k++] = p;
}
p+=2;
} long result = 0;
double minvalue=10000;
long euler = 0 ;
for(long n = 2;n<Max_n;n++){
euler = phi3(n,Prime);
if(isPerm(euler,n)){
double temp = n/(euler*1.0);
if(temp<minvalue){
minvalue = temp;
result = n;
}
}
}
System.out.println(result); }
结果:
// 8319823
// running time=58s426ms
时间少了30s
这里缺点是要取多少个素数?
我从100,200,300,400,在400的时候结果是正确的
全部程序:
package project61; public class P70{
void run2(){
long Max_n = 10000000;
int Prime_len=400;
int[] Prime = new int[Prime_len];
Prime[0] = 2;
int k=1;
int p=1;
while(k<Prime_len){
if(isPrime(p)){
Prime[k++] = p;
}
p+=2;
} long result = 0;
double minvalue=10000;
long euler = 0 ;
for(long n = 2;n<Max_n;n++){
euler = phi3(n,Prime);
if(isPerm(euler,n)){
double temp = n/(euler*1.0);
if(temp<minvalue){
minvalue = temp;
result = n;
}
}
}
System.out.println(result); }
// 8319823
// running time=58s426ms long phi3(long n,int[] Prime){
long res = 0;
if(n==1) return 0;
int len = Prime.length;
res = n;
int i=0; //&& n>=Prime[i]
while(i <len && n>=Prime[i] ){
if(n%Prime[i]==0){
res/=Prime[i];
res*=(Prime[i]-1);
while(n%Prime[i]==0){
n/=Prime[i];
}
}
i++;
}
if(n>1){
res/=n;
res*=(n-1);
}
return res;
}
void run(){
long Max_n = 10000000;
long result = 0;
double minvalue=10000;
long euler = 0 ;
for(long n = 2;n<Max_n;n++){
euler = phi2(n);
if(isPerm(euler,n)){
double temp = n/(euler*1.0);
if(temp<minvalue){
minvalue = temp;
result = n;
}
}
}
System.out.println(result);
}
boolean isPerm(long a,long b){
int[] label = new int[10];
while(a!=0&&b!=0){
label[(int) (a%10)]+=1;
label[(int) (b%10)]-=1;
a/=10;
b/=10; }
if(a!=0) return false;
if(b!=0) return false;
for(int i=0;i<10;i++)
if(label[i]!=0)
return false;
return true;
} boolean isPrime(int n){
if(n==2||n==3||n==5||n==7||n==11) return true;
if(n<2||n%2==0||n%3==0) return false;
for(int i=5;i<=Math.sqrt(n)+1;i++){
if(n%i==0) return false;
}
return true;
} long phi2(long n){
long res = 0;
if(n==1) return 0;
int pi=2;
res = n;
while(pi*pi <=n){
if(n%pi==0){
res/=pi;
res*=(pi-1);
while(n%pi==0){
n/=pi;
}
}
pi++;
}
if(n>1){
res/=n;
res*=(n-1);
}
return res;
}
// 8319823
// running time=89s43ms
// public static void main(String[] args){ long start = System.currentTimeMillis();
new P70().run2();
long end = System.currentTimeMillis();
long time = end - start;
System.out.println("running time="+time/1000+"s"+time%1000+"ms"); }
}
一直感觉这个程序还可以优化
是否还有好的算法
在解题论坛上,有个解答程序看不懂,然后我就敲成java,用时8s。。。快了好多的
long[] cal_phi(){
long[] phi=new long[Max_n+1];
for(int i=1;i<Max_n;i++){
phi[i]+=i;
for(int j=2*i;j<Max_n;j+=i)
phi[j]-=phi[i];
}
return phi;
}
8s
这个是产生Max_n的所以的欧拉函数,这个虽然是两层循环,但是效率很高的,变量完就得到欧拉函数值
下面一种形式:
long[] cal_phi2(){
long[] phi = new long[Max_n+1];
for(int i=1;i<Max_n;i++)
phi[i] = i;
for(int i=2;i<Max_n;i++){
if(phi[i]==i)
phi[i] = i - 1;
for(int j=2*i;j<Max_n;j+=i)
phi[j]*=(1-1.0/i);
}
return phi;
}
15s
// 8319823
// running time=14s577ms
至于上面是何总介绍?我还不明白
在解题报告中,还有一种是利用到素数,奇数,偶数还不明白为什么的。
全部程序:
package project61; public class P70_1{
int Max_n = 10000000;
void run(){
long[] phi=cal_phi();
int result = 0;
double minvalue = 10000.0;
for(int n=2;n<Max_n;n++){
long euler = phi[n];
if(isPerm(euler,n)){
double temp = n/(euler*1.0);
if(temp<minvalue){
minvalue = temp;
result = n;
}
} }
System.out.println(result); } long[] cal_phi(){
long[] phi=new long[Max_n+1];
for(int i=1;i<Max_n;i++){
phi[i]+=i;
for(int j=2*i;j<Max_n;j+=i)
phi[j]-=phi[i];
}
return phi;
}
// 8319823
// running time=7s855ms
long[] cal_phi2(){
long[] phi = new long[Max_n+1];
for(int i=1;i<Max_n;i++)
phi[i] = i;
for(int i=2;i<Max_n;i++){
if(phi[i]==i)
phi[i] = i - 1;
for(int j=2*i;j<Max_n;j+=i)
phi[j]*=(1-1.0/i);
}
return phi;
}
// 8319823
// running time=14s577ms
boolean isPerm(long a,long b){
int[] label = new int[10];
while(a!=0&&b!=0){
label[(int) (a%10)]+=1;
label[(int) (b%10)]-=1;
a/=10;
b/=10; }
if(a!=0) return false;
if(b!=0) return false;
for(int i=0;i<10;i++)
if(label[i]!=0)
return false;
return true;
} boolean isPrime(long n){
if(n==2||n==3||n==5||n==7) return true;
if(n<2||n%2==0||n%3==0) return false;
for(int i=5;i<=Math.sqrt(n);i++)
if(n%i==0) return false;
return true;
}
public static void main(String[] args){
long start = System.currentTimeMillis();
new P70_1().run();
long end = System.currentTimeMillis();
long time = end - start;
System.out.println("running time="+time/1000+"s"+time%1000+"ms");
}
}
上面的程序用Python实习时间好长的
欧拉工程第70题:Totient permutation的更多相关文章
- 欧拉工程第69题:Totient maximum
题目链接 欧拉函数φ(n)(有时也叫做phi函数)可以用来计算小于n 的数字中与n互质的数字的个数. 当n小于1,000,000时候,n/φ(n)最大值时候的n. 欧拉函数维基百科链接 这里的是p是n ...
- 欧拉工程第52题:Permuted multiples
题目链接 题目: 125874和它的二倍,251748, 包含着同样的数字,只是顺序不同. 找出最小的正整数x,使得 2x, 3x, 4x, 5x, 和6x都包含同样的数字. 这个题目相对比较简单 暴 ...
- 欧拉工程第68题:Magic 5-gon ring
题目链接 任意一条线上的三个数的和都等于9,顺时针,从最小的外圈开始,得到的序列是:432621213 和 序列 9位的字符串:三角环所能形成的最大字符串为432621513. ...
- 欧拉工程第67题:Maximum path sum II
By starting at the top of the triangle below and moving to adjacent numbers on the row below, the ma ...
- 欧拉工程第66题:Diophantine equation
题目链接 脑补知识:佩尔方差 上面说的貌似很明白,最小的i,对应最小的解 然而我理解成,一个循环的解了,然后就是搞不对,后来,仔细看+手工推导发现了问题.i从0开始变量,知道第一个满足等式的解就是最小 ...
- 欧拉工程第65题:Convergents of e
题目链接 现在做这个题目真是千万只草泥马在心中路过 这个与上面一题差不多 这个题目是求e的第100个分数表达式中分子的各位数之和 What is most surprising is that the ...
- 欧拉工程第56题:Powerful digit sum
题目链接 Java程序 package projecteuler51to60; import java.math.BigInteger; import java.util.Iterator; im ...
- 欧拉工程第55题:Lychrel numbers
package projecteuler51to60; import java.math.BigInteger; import java.util.Iterator; import java.util ...
- 欧拉工程第54题:Poker hands
package projecteuler51to60; import java.awt.peer.SystemTrayPeer; import java.io.BufferedReader; impo ...
随机推荐
- js设计模式(10)---观察者模式
0.前言 最近好多烦心事,由于自己的拖延懒惰造成事情堆积如山,看来最近得勤快些了,不然真的会死的很惨. 1.观察者模式是什么 又叫做发布者订阅者模式(publish/Subscribe),用来确定对象 ...
- CSS练习
看到一个CSS参考手册:http://css.doyoe.com/ 感谢感谢!
- PySide 简易教程<三>-------动手写起来
到目前为止,已经接触的Pyside的界面元素有如下几个:QWidget.QPushButton.QLabel.本次再介绍两个tooltip和messagebox.tooltip是一个鼠标悬浮提示信息, ...
- 《零成本实现Web自动化测试--基于Selenium》 第四章 Selenium 命令
Selenium 命令,通常被称为Selenese,由一系列运行测试案例所需要的命令构成.按顺序排列这些命令就构成了测试脚本. 一. 验证颜面元素 1.Assertion或者Verification ...
- Sales_item
#ifndef SALESITEM_H #define SALESITEM_H // Definition of Sales_item class and related functions goes ...
- boost muti-thread
背景 • 今天互联网应用服务程序普遍使用多线程来提高与多客户链接时的效率:为了达到最大的吞吐量,事务服务器在单独的线程上运行服务程序: GUI应用程序将那些费时,复杂的处理以线程的形式单独 ...
- 微软职位内部推荐-Senior Software Development Engineer
微软近期Open的职位: Job posting title: Senior Software Development Engineer Location: China, Beijing Divisi ...
- “Microsoft Visual Studio遇到了问题,需要关闭”解决办法
运行VS2008,打开项目,弹出错误界面 . 解决办法:将项目中的所有设计窗体关闭并保存,重新打开就OK~
- OC面向对象继承关系和组合关系笔记
继承关系是描述类和类之间的关系,两个类分别称为子类和父类,子类继承了父类,子类就拥有了父类的属性和方法: 继承的关系特点描述出来就是:** “是” ** (例如:学生类 是 人类) 组合关系描述的语 ...
- Careercup - Facebook面试题 - 4907555595747328
2014-05-02 07:49 题目链接 原题: Given a set of n points (coordinate in 2d plane) within a rectangular spac ...