Triangle, square, pentagonal, hexagonal, heptagonal, and octagonal numbers are all figurate (polygonal) numbers and are generated by the following formulae:

Triangle   P3,n=n(n+1)/2   1, 3, 6, 10, 15, ...
Square   P4,n=n2   1, 4, 9, 16, 25, ...
Pentagonal   P5,n=n(3n−1)/2   1, 5, 12, 22, 35, ...
Hexagonal   P6,n=n(2n−1)   1, 6, 15, 28, 45, ...
Heptagonal   P7,n=n(5n−3)/2   1, 7, 18, 34, 55, ...
Octagonal   P8,n=n(3n−2)   1, 8, 21, 40, 65, ...

The ordered set of three 4-digit numbers: 8128, 2882, 8281, has three interesting properties.

  1. The set is cyclic, in that the last two digits of each number is the first two digits of the next number (including the last number with the first).
  2. Each polygonal type: triangle (P3,127=8128), square (P4,91=8281), and pentagonal (P5,44=2882), is represented by a different number in the set.
  3. This is the only set of 4-digit numbers with this property.

Find the sum of the only ordered set of six cyclic 4-digit numbers for which each polygonal type: triangle, square, pentagonal, hexagonal, heptagonal, and octagonal, is represented
by a different number in the set.

又暴力破解了一次ㄟ( ▔, ▔ )ㄏ

一開始没看清题意,我以为这些数依次是满足triangle, square, pentagonal, hexagonal, heptagonal, and octagonal。结果发现无解┑( ̄Д  ̄)┍

#include <iostream>
#include <string>
#include <vector>
#include <unordered_map>
#include <time.h>
using namespace std; int triangle[100];
int pentagonal[10000];
int hextagonal[10000];
int heptagonal[10000];
int octagonal[10000];
int tri_count = 0; void getTriangle()
{
int count = 0;
for (int i = 1; i <= 200; i++)
{
int num = i*(i + 1) / 2;
if (num >1000&&num<10000)
triangle[count++] = num;
}
tri_count = count;
} bool isSqure(int n)
{
int i = sqrt(n);
if (i*i == n&&n>1000&&n<10000)
return true;
return false;
} void getPentagonal()
{
for (int i = 1; i <= 200; i++)
{
int num = i*(3 * i - 1) / 2;
if (num > 1000 && num < 10000)
pentagonal[num] = 1;
}
} bool isPentagonal(int n)
{
if (pentagonal[n] == 1)
return true;
return false;
} void getHexagonal()
{
for (int i = 1; i <= 200; i++)
{
int num = i*(2 * i - 1);
if (num>1000 && num < 10000)
hextagonal[num] = 1;
}
} bool isHexagonal(int n)
{
if (hextagonal[n] == 1)
return true;
return false;
} void getHeptagonal()
{
for (int i = 1; i <= 200; i++)
{
int num = i*(5 * i - 3) / 2;
if (num > 1000 && num < 10000)
heptagonal[num] = 1;
}
} bool isHeptagonal(int n)
{
if (heptagonal[n] == 1)
return true;
return false;
} void getOctagonal()
{
for (int i = 1; i <= 200; i++)
{
int num = i*(3 * i - 2);
if (num > 1000 && num < 10000)
octagonal[num] = 1;
}
} bool isOctagonal(int n)
{
if (octagonal[n] == 1)
return true;
return false;
} bool(*figurate[5])(int) = { isSqure, isPentagonal, isHexagonal, isHeptagonal, isOctagonal }; vector<int> GetRandomSequence()
{
unordered_map<int, int>tab;
vector<int>res;
int num;
for (int i = 0; i < 5; i++)
{
do{
num = rand() % 5;
} while (tab.find(num) != tab.end());
tab.insert(make_pair(num, 1));
res.push_back(num);
}
return res;
} int check()
{
int sum = 0;
srand((int)time(0));
vector<int>rs = GetRandomSequence();
for (int i = 0; i < tri_count; i++)
{
int a = triangle[i] / 100;
int b = triangle[i] % 100;
for (int s = 10; s <= 99; s++)
{
if ((*figurate[rs[0]])(b * 100 + s))
{
for (int p = 10; p <= 99; p++)
{
if ((*figurate[rs[1]])(s * 100 + p))
{
for (int hx = 10; hx <= 99; hx++)
{
if ((*figurate[rs[2]])(p * 100 + hx))
{
for (int hp = 10; hp <= 99; hp++)
{
if ((*figurate[rs[3]])(hx * 100 + hp))
{
if ((*figurate[rs[4]])(hp * 100 + a))
{
sum = triangle[i] + b * 100 + s + s * 100 + p + p * 100 + hx + hx * 100 + hp + hp * 100 + a;
return sum;
}
}
}
}
}
}
}
}
}
}
return -1;
} int main()
{
memset(pentagonal, 0, sizeof(pentagonal));
memset(hextagonal, 0, sizeof(hextagonal));
memset(heptagonal, 0, sizeof(heptagonal));
memset(octagonal, 0, sizeof(octagonal)); getTriangle();
getPentagonal();
getHexagonal();
getHeptagonal();
getOctagonal(); int flag;
while (true)
{
flag = check();
if (flag != -1)
break;
} cout << flag << endl; system("pause");
return 0;
}

把那个随机生成全排列换成next_permutation也是能搞出来的。

Project Euler:Problem 61 Cyclical figurate numbers的更多相关文章

  1. Project Euler:Problem 42 Coded triangle numbers

    The nth term of the sequence of triangle numbers is given by, tn = ½n(n+1); so the first ten triangl ...

  2. Project Euler:Problem 55 Lychrel numbers

    If we take 47, reverse and add, 47 + 74 = 121, which is palindromic. Not all numbers produce palindr ...

  3. Project Euler:Problem 88 Product-sum numbers

    A natural number, N, that can be written as the sum and product of a given set of at least two natur ...

  4. Project Euler:Problem 87 Prime power triples

    The smallest number expressible as the sum of a prime square, prime cube, and prime fourth power is ...

  5. Project Euler:Problem 89 Roman numerals

    For a number written in Roman numerals to be considered valid there are basic rules which must be fo ...

  6. Project Euler:Problem 93 Arithmetic expressions

    By using each of the digits from the set, {1, 2, 3, 4}, exactly once, and making use of the four ari ...

  7. Project Euler:Problem 28 Number spiral diagonals

    Starting with the number 1 and moving to the right in a clockwise direction a 5 by 5 spiral is forme ...

  8. Project Euler:Problem 47 Distinct primes factors

    The first two consecutive numbers to have two distinct prime factors are: 14 = 2 × 7 15 = 3 × 5 The ...

  9. Project Euler:Problem 63 Powerful digit counts

    The 5-digit number, 16807=75, is also a fifth power. Similarly, the 9-digit number, 134217728=89, is ...

随机推荐

  1. css--css选择器,伪类

    前戏 前面我们说过CSS规则由选择器和声明组成,我们要给标签设置属性,那我们就要找到对应的标签,CSS选择器可以帮我们找到我们需要的标签 css选择器有: 标签选择器 类选择器 ID选择器 全局选择器 ...

  2. Ztree 多选,显示勾选的路径

    项目要求,需要向后台传递已经勾选的路径,如 l1-a, l1-l3-c,l1-l3-d;(如果是全选状态则只传递全选状态的路径,不传子节点). 具体可以参考jQ  Ztree 的 v3.5 版本 Me ...

  3. Yii1 获取当前请求的url

    echo Yii::app()->getRequest()->getUrl();

  4. Django reverse函数

    1.总urls.py内容如下: from django.contrib import admin from django.urls import path from django.conf.urls ...

  5. NLog在asp.net core中的应用

    Asp.net core中,自带的Log是在当selfhost运行时,在控制台中输出,不便于查阅,如果用一个log架框,把日志持久化,便于查询. NLog是一个免费的日志记录框架,专门为.net平台下 ...

  6. uC/OSii之任务划分

    满足以下几个指标将会使软件设计比较简洁高效. 满足实时性指标 任务数目合理 简化软件系统 降低资源需求 设备依赖性划分:主要是对于一些输入输出设备进行划分,输入输出设备分为主动型和被动性.主动型有一个 ...

  7. shell for mysql backup in linux

    今天上班只有一台linux系统,就学着在linux上写了个脚本,没啥技术含量 省得每天敲代码备份 没有设置自动备份时间,这里可以参照 http://www.th7.cn/db/mysql/201305 ...

  8. XTUOJ 15503 - C

    15503 - C Accepted: 6    Submissions: 27    Time Limit: 3000 ms    Memory Limit: 1048576 KB 在解决了小女孩的 ...

  9. [Go]条件语句

    package main import ( "io/ioutil" "fmt" ) //条件语句 //if的条件语句不需要括号 //if的条件里可以赋值,if的 ...

  10. vscode安装插件

    十分简单,知道名字叫啥后,直接搜索,安装,就完了,还可以查看自己已经安装了哪些插件. step1 如图.png step2 image.png step 3 去网上查找想要安装的插件的名字 step ...