Almost equilateral triangles

It is easily proved that no equilateral triangle exists with integral length sides and integral area. However, the almost equilateral triangle 5-5-6 has an area of 12 square units.

We shall define an almost equilateral triangle to be a triangle for which two sides are equal and the third differs by no more than one unit.

Find the sum of the perimeters of all almost equilateral triangles with integral side lengths and area and whose perimeters do not exceed one billion (1,000,000,000).


几乎等边的三角形

可以证明,不存在边长为整数的等边三角形其面积也是整数。但是,存在几乎等边的三角形 5-5-6,其面积恰好为12。

我们定义几乎等边的三角形是有两条边一样长,且第三边与这两边最多相差1的三角形。

对于周长不超过10亿的三角形中,找出边长和面积都为整数的近似等边三角形的周长和。

解题

这个直接暴力是可以的

先说:几乎等边三角形有两种,边长可以是:a、a、a+1 和a、a、a-1都满足三角形任意两边之和大于第三边,两边之差小于第三边。

下面就是这么根据三边怎么计算面积,并判断面积是整数

我第一想到的是海伦公式

当边长是:a、a、a+1的时候

判断根号下面的平方数,再判断结果S能够被4整除

当边长是:a、a、a-1的时候

这个和上面也一样了。

这个一个变量时间比较长,表示很伤不起

JAVA

package Level3;

public class PE094{
public static void run(){
long MAX = 1000000000;
long L1 = 0;
long L2 = 0;
for(long a =2;;a++){
if(3*a+1>MAX || 3*a-1>MAX){
break;
}
if(a>=2){
if(areaIsInteger(a,a+1) == true)
L1 +=3*a + 1;
}
if(a>=3){
if(areaIsInteger(a,a-1) == true)
L2 +=3*a - 1;
}
} System.out.println("L1: "+L1);
System.out.println("L2: "+L2);
System.out.println(L1 + L2);
}
// L1: 109552588
// L2: 408855758
// 518408346
// running time=124s662ms
public static boolean areaIsInteger(long a ,long c){
// c = a + 1
if( c == a+1){
long s1 = (3*a + 1)*(a-1);
if((long)Math.sqrt(s1) * (long)Math.sqrt(s1) !=s1)
return false;
long s2 = (long) ((a+1)*Math.sqrt(s1));
if(s2%4 != 0)
return false;
else{
return true;
}
}
if( c== a-1){
long s1 = (3*a -1)*(a+1);
if((long)Math.sqrt(s1)*(long)Math.sqrt(s1) !=s1)
return false;
long s2 = (long)((a-1)*Math.sqrt(s1));
if(s2%4 != 0)
return false;
else{
return true;
}
}
return true;
}
public static void main(String[] args) {
long t0 = System.currentTimeMillis();
run();
long t1 = System.currentTimeMillis();
long t = t1 - t0;
System.out.println("running time="+t/1000+"s"+t%1000+"ms"); }
}

上面根据海伦公式判断面积是否是整数,比较复杂。

直接根据三角形底乘高除以2比较简单点。

公式:

判断面积是否是整数的程序如下

    public static boolean areaIsInteger2(long a,long c){
if(c*c%4!=0)
return false;
long s0 = a*a-c*c/4;
long sqrt = (long)(Math.sqrt(s0));
if(sqrt * sqrt != s0)
return false;
long s = c*sqrt;
if(s%2!=0)
return false;
return true;
}

这个比较简单了,运行时间还是比较快的

//    L1: 109552588
// L2: 408855758
// 518408346
// running time=48s879ms

Project Euler 94:Almost equilateral triangles 几乎等边的三角形的更多相关文章

  1. Python练习题 040:Project Euler 012:有超过500个因子的三角形数

    本题来自 Project Euler 第12题:https://projecteuler.net/problem=12 # Project Euler: Problem 12: Highly divi ...

  2. Project Euler 91:Right triangles with integer coordinates 格点直角三角形

    Right triangles with integer coordinates The points P (x1, y1) and Q (x2, y2) are plotted at integer ...

  3. hackerrank Project Euler #210: Obtuse Angled Triangles

    传送门 做出一个好几个星期屯下来的题目的感觉就是一个字: 爽! 上图的黄点部分就是我们需要求的点 两边的部分很好算 求圆的地方有一个优化,由于圆心是整数点,我们可以把圆分为下面几个部分,阴影部分最难算 ...

  4. Project Euler 39 Integer right triangles( 素勾股数 )

    题意:若三边长 { a , b , c } 均为整数的直角三角形周长为 p ,当 p = 120 时,恰好存在三个不同的解:{ 20 , 48 , 52 } , { 24 , 45 , 51 } , ...

  5. Python练习题 039:Project Euler 011:网格中4个数字的最大乘积

    本题来自 Project Euler 第11题:https://projecteuler.net/problem=11 # Project Euler: Problem 10: Largest pro ...

  6. [project euler] program 4

    上一次接触 project euler 还是2011年的事情,做了前三道题,后来被第四题卡住了,前面几题的代码也没有保留下来. 今天试着暴力破解了一下,代码如下: (我大概是第 172,719 个解出 ...

  7. Python练习题 029:Project Euler 001:3和5的倍数

    开始做 Project Euler 的练习题.网站上总共有565题,真是个大题库啊! # Project Euler, Problem 1: Multiples of 3 and 5 # If we ...

  8. Project Euler 9

    题意:三个正整数a + b + c = 1000,a*a + b*b = c*c.求a*b*c. 解法:可以暴力枚举,但是也有数学方法. 首先,a,b,c中肯定有至少一个为偶数,否则和不可能为以上两个 ...

  9. Project Euler 44: Find the smallest pair of pentagonal numbers whose sum and difference is pentagonal.

    In Problem 42 we dealt with triangular problems, in Problem 44 of Project Euler we deal with pentago ...

随机推荐

  1. 当年的笔记_apache配置虚拟主机

    下午需要,在网上找了一堆,没找到合适的,翻出来自己当年的笔记,还是自己记的容易理解. 解决方案1:通过端口来区分 1>添加一个虚拟主机1.在d盘下新建www目录,如:d:/www. 2.修改ht ...

  2. 团队自动化环境搭建与管理--php博弈

    我是方少,很开心与大家日后与大家交流技术上面的一些想法和一些业务上的分享.以前从来没写过博客,因为觉得不重要吧,如今觉得有必要沉淀一些想法和回忆.好了费话不多说. 先上图: 业务问题:在每次新伙伴加入 ...

  3. phpcms v9 源码解析(4)content模块下的index.php文件的init()方法解析

    在了解index.php中的init函数的时候,让我们先看看最开始的几行代码 1-5  第二行, defined('IN_PHPCMS') or exit('Nopermission resource ...

  4. ECSHOP订单自动确认

    1.运行sql代码,生成数据库 CREATE TABLE `ecs_order_auto_confirm` ( `id` INT() UNSIGNED NOT NULL AUTO_INCREMENT, ...

  5. python with语句上下文管理的两种实现方法

    在编程中会经常碰到这种情况:有一个特殊的语句块,在执行这个语句块之前需要先执行一些准备动作:当语句块执行完成后,需要继续执行一些收尾动作.例如,文件读写后需要关闭,数据库读写完毕需要关闭连接,资源的加 ...

  6. 【netstream】探索数据传输对象1

    什么是“从当前流中读取一个字符串.字符串有长度前缀,一次 7 位地被编码为整数.” 来探索一下: 写一段简单的程序: FileStream fs= new FileStream("d:\\q ...

  7. 【python】 入门 - 函数式编程

    函数式编程的一个特点就是,允许把函数本身作为参数传入另一个函数,还允许返回一个函数 http://www.liaoxuefeng.com/wiki/001374738125095c955c1e6d8b ...

  8. java collections

    http://www.docjar.com/docs/api/java/util/technotes/guides/collections/changes4.html http://www.docja ...

  9. Jquery ajax请求导出Excel表格

    直接贴代码吧 $("#btn-export").click(function(){ var exportExcel = "export_excel"; data ...

  10. 3243: [Noi2013]向量内积 - BZOJ

    Description 两个d 维向量A=[a1,a2,...,ad]与B=[b1,b2,...,bd]的内积为其相对应维度的权值的乘积和,即: 现有 n 个d 维向量x1,...,xn ,小喵喵想知 ...