本文通过一系列有趣的数学问题和编程挑战，展示了数学与编程的紧密联系。从高斯日记中计算获得博士学位的日期，到寻找马虎算式中的巧合规律，再到解决阶梯上的不同上法、黄金分割数的精确表示，以及各种算法设计。文章深入浅出，带领读者体验数学之美和编程之妙。
1.题目标题: 高斯日记
大数学家高斯有个好习惯：无论如何都要记日记。
他的日记有个与众不同的地方，他从不注明年月日，而是用一个整数代替，比如：4210
后来人们知道，那个整数就是日期，它表示那一天是高斯出生后的第几天。这或许也是个好习惯，它时时刻刻提醒着主人：日子又过去一天，还有多少时光可以用于浪费呢？
高斯出生于：1777年4月30日。
在高斯发现的一个重要定理的日记上标注着：5343，因此可算出那天是：1791年12月15日。
高斯获得博士学位的那天日记上标着：8113
请你算出高斯获得博士学位的年月日。
提交答案的格式是：yyyy-mm-dd, 例如：1980-03-21
请严格按照格式，通过浏览器提交答案。
注意：只提交这个日期，不要写其它附加内容，比如：说明性的文字。
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long LL;
int mon[13] = {0,31,28,31,30,31,30,31,31,30,31,30,31};//因为数组默认开始为0，但now.m是从1开始，所以mon数组第一个为0
struct Date{
int y,m, d;
}now;
bool is_prime(int year)//特判闰年
{
return (year%4 == 0 && year % 100 != 0) || (year % 400 == 0);
}
int main()
{
now.y = 1777,now.m = 4, now.d = 30;
int n ;
cin >> n;
n--;//不包含那一天
while(n--)
{
now.d++;
if(now.m == 2&& now.d == 28)
{
if(is_prime(now.y))
mon[2] = 29;
else mon[2] = 28;
}
if(now.d > mon[now.m]){
now.m++;
if(now.m > 12){
now.m= 1;
now.y++;
}
now.d = 1;
}
}
cout << now.y << " " << now.m <<" " << now.d;
return 0;
}
答案：1977-07-16
2、马虎的算式
小明是个急性子，上小学的时候经常把老师写在黑板上的题目抄错了。
有一次，老师出的题目是：36 x 495 = ?
他却给抄成了：396 x 45 = ?
但结果却很戏剧性，他的答案竟然是对的！！
因为 36 * 495 = 396 * 45 = 17820
类似这样的巧合情况可能还有很多，比如：27 * 594 = 297 * 54
假设 a b c d e 代表1~9不同的5个数字（注意是各不相同的数字，且不含0）
能满足形如： ab * cde = adb * ce 这样的算式一共有多少种呢？
请你利用计算机的优势寻找所有的可能，并回答不同算式的种类数。
满足乘法交换律的算式计为不同的种类，所以答案肯定是个偶数。
解题方法1：暴力枚举
#include <iostream>
#include<algorithm>
using namespace std;
int main()
{
int a, b,c,d,e;
int i, j,sum= 0;
for(a = 1; a <=9; a++)
for(b = 1; b <=9; b++)
for(c = 1; c <= 9; c++)
for(d = 1; d <= 9; d++)
for(e = 1; e <= 9; e++)
{
i = (a * 10 + b)*(c*100 + d * 10 + e);
j = (a * 100 + d * 10 + b)*(c * 10 + e);
if( i == j && a != b && a!=c&& a!=d && a != e && b != c && b != d && b != e && c != d && c != e && d != e)
sum++;
}
cout << sum;
return 0;
}
方法二：全排列
#include <iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
int main()
{
int s[9]={1,2,3,4,5,6,7,8,9};
int a, b,c,d,sum = 0;
do{
a = s[0]*10 + s[1];
b = s[2]*100 + s[3]*10 + s[4];
c = s[0]*100 + s[3]*10 + s[1];
d = s[2]*10 + s[4];
if(a * b ==c * d)
sum++;
}while(next_permutation(s, s+9));
cout <<sum/24;//除以24是因为题目要求只选五个数字，而计算的时候我是按9个数字算的，所以会由重复的，且重复次数是（9-5）！=4*3*2*1=24；
return 0;
}
3、第39级台阶
小明刚刚看完电影《第39级台阶》，离开电影院的时候，他数了数礼堂前的台阶数，恰好是39级!
站在台阶前，他突然又想着一个问题：
如果我每一步只能迈上1个或2个台阶。先迈左脚，然后左右交替，最后一步是迈右脚，也就是说一共要走偶数步。那么，上完39级台阶，有多少种不同的上法呢？
请你利用计算机的优势，帮助小明寻找答案。
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
int res = 0;
void find1(int n, int step)//递归
{
if(n <0) return;
if(n == 0)
{
if(step%2 == 0) res++;//阶梯走完后，判断是否为偶数步，若是则这种方法成立
else return;
}
find1(n - 1, step + 1);
find1(n - 2,step + 1);
}
int main()
{
find1(39, 0);
cout <<res;
return 0;
}
4.黄金连分数
黄金分割数0.61803… 是个无理数，这个常数十分重要，在许多工程问题中会出现。有时需要把这个数字求得很精确。
对于某些精密工程，常数的精度很重要。也许你听说过哈勃太空望远镜，它首次升空后就发现了一处人工加工错误，对那样一个庞然大物，其实只是镜面加工时有比头发丝还细许多倍的一处错误而已，却使它成了“近视眼”!!
言归正传，我们如何求得黄金分割数的尽可能精确的值呢？有许多方法。
比较简单的一种是用连分数：
1
黄金数 = ---------------------
1
1 + -----------------
1
1 + -------------
1
1 + ---------
1 + ...
这个连分数计算的“层数”越多，它的值越接近黄金分割数。
请你利用这一特性，求出黄金分割数的足够精确值，要求四舍五入到小数点后100位。
小数点后3位的值为：0.618
小数点后4位的值为：0.6180
小数点后5位的值为：0.61803
小数点后7位的值为：0.6180340
（注意尾部的0，不能忽略）
你的任务是：写出精确到小数点后100位精度的黄金分割值。
注意：尾数的四舍五入！ 尾数是0也要保留！
解题思考：
- 转化为斐波那契数列的n和n+1项
- n取多少？再增加n,直到小数点后的101位没有变化
- 不能用C语言定义的整数型直接运算，自己写出相对应的高精对加法和除法
#include<iostream>
#include<string>
#include<algorithm>
using namespace std;
//思路：由黄金连分数可转化为求解斐波那契数，运用大数运算求解小数点后一百位
//n表示斐波那契数列第一百项
int n = 100;
int cmp(string a, string b)
{
if (a == b) return 0;
else if (a.length() < b.length()) return -1;
else if (a.length() > b.length()) return 1;
else
{
if (a > b) return 1;
if (a < b) return -1;
return -1;
}
}
//大数相加，较大斐波那契数无法算出，只能通过大数计算
string add(string a, string b)
{
//去掉开头的0,substr(a,n) 返回第a位开始长度为n的字符串
//find_first_not_of（"abc"）返回最先匹配到abc任意一个字符的最后位置
a = a.substr(a.find_first_not_of('0'));
b = b.substr(a.find_first_not_of('0'));
long long lenA = a.length();
long long lenB = b.length();
long long len = max(lenA, lenB) + 10;
//反转， 便于从最低位逐位求和
reverse(a.begin(), a.end());
reverse(b.begin(), b.end());
string ans(len, '0');
for (int i = 0; i < lenA; i++)
{
ans[i] = a[i];
}
//通过模拟手工计算来进行大数相加，temp用来存储上下两位相加的结果
int temp = 0;
for (int i = 0; i < len; i++)
{
if (i < lenB)
{
temp += ((ans[i] - '0') + (b[i] - '0'));
ans[i] = temp % 10 + '0';
temp /= 10;
}
else
{
temp += (ans[i] - '0');
ans[i] = temp % 10 + '0';
temp /= 10;
}
}
//ans存储的结果最后要反转回来去掉开头的0
reverse(ans.begin(), ans.end());
return ans.substr(ans.find_first_not_of('0'));
}
//大数相减
string substract(string a, string b)
{
//a必须大于b,实际上在divide（）中，a已经>=b了
reverse(a.begin(), a.end());
reverse(b.begin(), b.end());
//将a复制给ans
string ans = a;
//依旧是模拟手工减法计算
for (int i = 0; i < b.length(); i++)
{
if (ans[i] >= b[i])
{
ans[i] = ans[i] - b[i] + '0';
}
//注意：在被减数不够减的情况下，需要向前借位，
//可能被借位为0也不够借，继续往前寻找不为0的位，通过一个while循环实现
else
{
int k = 1;
while (ans[i + k] == '0')
{
//0被借位变为9
ans[i + k] = '9';
k++;
}
//最终i+k位不为0可以被借-1，i位+10
ans[i + k] = ans[i + k] - 1;
ans[i] = ans[i] + 10 - b[i] + '0';
}
}
reverse(ans.begin(), ans.end());
return ans.substr(ans.find_first_not_of('0'));
}
//大数相除
string divide(string a, string b)
{
//前提 a < b 实际题目已满足
//大数除法实际上是模拟除法运算结合大数减法
string ans = "0.";
for (int i = 0; i < 101; i++) //保留101项，保证四舍五入
{
// (a*10)/b = t 不过其中除法用减法substract代替
a.append(1,'0');
int t = 0;
while (cmp(a, b) >= 0)
{
a = substract(a, b);
t++;
}
ans.append(1,t + '0');
}
return ans;
}
int main()
{
//斐波那契数列前两项
string a = "1";
string b = "1";
for (int k = 0; k < 40; k++)
{
//求解斐波那契数列100-140各项与前一项的比值（所谓黄金连分数）
for (int i = 3; i < n+k; i++)
{
string temp = b;
b = add(a, b);
a = temp;
}
string ans = divide(a, b);
cout << 100+k << "项 " <<ans << endl;
}
return 0;
}
5.前缀判断
如下的代码判断 needle_start指向的串是否为haystack_start指向的串的前缀，如不是，则返回NULL。
比如："abcd1234" 就包含了 "abc" 为前缀
char* prefix(char* haystack_start, char* needle_start)
{
char* haystack = haystack_start;
char* needle = needle_start;
while(*haystack && *needle){
if(______________________________) return NULL; //填空位置
}
if(*needle) return NULL;
return haystack_start;
}
#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
char* prefix(char* haystack_start, char* needle_start)
{
char* haystack = haystack_start;
char* needle = needle_start;
while(*haystack && *needle)//两个串都没有越界
{
if(*(haystack++) !=*(needle++)) return NULL; //填空位置,先取内容，再指针移动
}
if(*needle) return NULL;//如果needle的长度大于haystack
return haystack_start;
}
int main()
{
char a[10] = "abcd1234";
char b[4] = "abc";
printf("%s", prefix(a, b));
return 0;
}
6.三部排序
一般的排序有许多经典算法，如快速排序、希尔排序等。
但实际应用时，经常会或多或少有一些特殊的要求。我们没必要套用那些经典算法，可以根据实际情况建立更好的解法。
比如，对一个整型数组中的数字进行分类排序：
使得负数都靠左端，正数都靠右端，0在中部。注意问题的特点是：负数区域和正数区域内并不要求有序。可以利用这个特点通过1次线性扫描就结束战斗!!
以下的程序实现了该目标。
其中x指向待排序的整型数组，len是数组的长度。
void sort3p(int* x, int len)
{
int p = 0;
int left = 0;
int right = len-1;
while(p<=right){
if(x[p]<0){
int t = x[left];
x[left] = x[p];
x[p] = t;
left++;
p++;
}
else if(x[p]>0){
int t = x[right];
x[right] = x[p];
x[p] = t;
right--;
}
else{
__________________________; //填空位置
}
}
}
如果给定数组：
25,18,-2,0,16,-5,33,21,0,19,-16,25,-3,0
则排序后为：
-3,-2,-16,-5,0,0,0,21,19,33,25,16,18,25
#include<iostream>
#include<algorithm>
#include<cstdio>
void sort3p(int* x, int len)
{
int p = 0;
int left = 0;
int right = len-1;
//不超过边界的情况下，比较x[p]是否大于0，如果大于0，交换到右边移动p
while(p<=right){
if(x[p]<0){
int t = x[left];
x[left] = x[p];
x[p] = t;
left++;
p++;
}
else if(x[p]>0){
int t = x[right];
x[right] = x[p];
x[p] = t;
right--;
}
else{
p++; //填空位置
}
}
}
int main()
{
int i;
int a[20]={25,18,-2,0,16,-5,33,21,0,19,-16,25,-3,0};
sort3p(a,14);
for(i=0;i<14;i++)
printf("%d ",a[i]);
return 0;
}
7.错误票据
某涉密单位下发了某种票据，并要在年终全部收回。
每张票据有唯一的ID号。全年所有票据的ID号是连续的，但ID的开始数码是随机选定的。
因为工作人员疏忽，在录入ID号的时候发生了一处错误，造成了某个ID断号，另外一个ID重号。
你的任务是通过编程，找出断号的ID和重号的ID。
假设断号不可能发生在最大和最小号。
要求程序首先输入一个整数N(N<100)表示后面数据行数。
接着读入N行数据。
每行数据长度不等，是用空格分开的若干个（不大于100个）正整数（不大于100000）
每个整数代表一个ID号。
要求程序输出1行，含两个整数m n，用空格分隔。
其中，m表示断号ID，n表示重号ID
例如：
用户输入：
2
5 6 8 11 9
10 12 9
则程序输出：
7 9
再例如：
用户输入：
6
164 178 108 109 180 155 141 159 104 182 179 118 137 184 115 124 125 129 168 196
172 189 127 107 112 192 103 131 133 169 158
128 102 110 148 139 157 140 195 197
185 152 135 106 123 173 122 136 174 191 145 116 151 143 175 120 161 134 162 190
149 138 142 146 199 126 165 156 153 193 144 166 170 121 171 132 101 194 187 188
113 130 176 154 177 120 117 150 114 183 186 181 100 163 160 167 147 198 111 119
则程序输出：
105 120
注意：本题难点在于获取用户输入的票据的ID号，因为每行输入的个数不同
解题思路：
- 首先用字符串来存储每一行的数据，然后通过计算每一行的每个数据（空格为标志）
- 将处理出来的数据存入数组中，同时记录数据出现次数
- 对其进行排序
- 若出现的次数>1 重复的数据+1 ；若出现的次数=0 断号+1
#include<iostream>
#include<algorithm>
#include<cstdio>
using namespace std;
/*将每一行当做字符串读入，然后经过计算放入对应的数组当中，*/
int f[100001];
int main()
{
int N, max = 0, min = 100000,start, num, i , miss, copy;
char a[700], c;
cin >> N;//数据行数
gets(a);
for(i = 0; i < N;i++)
{
gets(a);
start = 0;
while(a[start] != '\0')//判断一个数字是否输入完
{
num = 0;
while(a[start] == ' ') start++;//去除空格
for(c = a[start]; c != ' ' && c != '\0';)
{
num = num * 10 + c - '0';//计算数的大小
c = a[++start];
}
f[num]++;//记录数据出现次数
if(num <min) min = num;//确定最小数
if(num >max) max = num;
while(a[start] == ' ') start++;
}
}
for(i = min ; i < max; i++)
{
if(f[i] == 0) miss = i;
if(f[i] > 1) copy = i;
}
printf("%d %d\n", miss, copy);
return 0;
}
8.翻硬币
小明正在玩一个“翻硬币”的游戏。
桌上放着排成一排的若干硬币。我们用 * 表示正面，用 o 表示反面（是小写字母，不是零）。
比如，可能情形是：**oo***oooo
如果同时翻转左边的两个硬币，则变为：oooo***oooo
现在小明的问题是：如果已知了初始状态和要达到的目标状态，每次只能同时翻转相邻的两个硬币,那么对特定的局面，最少要翻动多少次呢？
我们约定：把翻动相邻的两个硬币叫做一步操作，那么要求：
程序输入：
两行等长的字符串，分别表示初始状态和要达到的目标状态。每行的长度<1000
程序输出：
一个整数，表示最小操作步数
例如：
用户输入：
**********
o****o****
程序应该输出：
5
再例如：
用户输入：
*o**o***o***
*o***o**o***
程序应该输出：
1
思路：
1.模拟翻硬币的过程
2.记录不相同的位置，规律：相邻两个不相同硬币的位置差即为这两个之间需要翻转的次数
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
int main()
{
int len,i,a[10]={0},j=0,sum=0;
char s1[10],s2[10];
scanf("%s",s1);
scanf("%s",s2);
len=strlen(s1);
for(i=0;i<len;i++)
{
if(s1[i]!=s2[i])//当两个串的出现不相等时，记录该点的位置
a[j++]=i;
}
for(i=0;i<j;i+=2)//两个两个地翻转硬币，翻转的次数为两个不相等位置之间的间距
{
sum+=a[i+1]-a[i];
}
printf("%d\n",sum);
return 0;
}
9.带分数
100 可以表示为带分数的形式：100 = 3 + 69258 / 714
还可以表示为：100 = 82 + 3546 / 197
注意特征：带分数中，数字1~9分别出现且只出现一次（不包含0）。
类似这样的带分数，100 有 11 种表示法。
题目要求：
从标准输入读入一个正整数N (N<1000*1000)
程序输出该数字用数码1~9不重复不遗漏地组成带分数表示的全部种数。
注意：不要求输出每个表示，只统计有多少表示法！
例如：
用户输入：
100
程序输出：
11
再例如：
用户输入：
105
程序输出：
6
深搜+全排列
#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstring>
#include<cmath>
using namespace std;
int a[10] = {0};
int b[10] = {0};
int x,y,z, cnt = 0;//cnt记录能够表示的方法总数
void dfs(int s, int n)
{
int i;
if(s >9)//当9位数全部排列完成后
{
//第一种表示方法 100 = 3 + 69258 / 714
x = a[1];
y = a[2]*10000 + a[3]*1000 + a[4]*100 + a[5]*10 + a[6];
z = a[7]* 100 + a[8]*10 + a[9];
if(y % z == 0)
{
if(x + y/z == n)
{
cnt++;
printf(" %d %d %d \n", x, y, z);
}
}
//第一种表示方法 100 = 82 + 3546 / 197
x =a[1]*10+a[2];
y = a[3]*1000+a[4]*100+a[5]*10+a[6];
z = a[7]*100+a[8]*10+a[9];
if(y % z == 0)
{
if(x + y/z == n)
{
cnt++;
printf(" %d %d %d \n", x, y, z);
}
}
}
else//若s<= 9
{
for(i=1;i<10;i++)
if( b[i]==0 )
{
b[i]=1;
a[s]=i;
dfs(s+1,n);
b[i]=0;
}
}
}
int main()
{
int n;
cin >> n;
dfs(1, n);
cout << cnt;
return 0;
}
10.连号区间数
小明这些天一直在思考这样一个奇怪而有趣的问题：
在1~N的某个全排列中有多少个连号区间呢？这里所说的连号区间的定义是：
如果区间[L, R] 里的所有元素（即此排列的第L个到第R个元素）递增排序后能得到一个长度为R-L+1的“连续”数列，则称这个区间连号区间。
当N很小的时候，小明可以很快地算出答案，但是当N变大的时候，问题就不是那么简单了，现在小明需要你的帮助。
输入格式：
第一行是一个正整数N (1 <= N <= 50000), 表示全排列的规模。
第二行是N个不同的数字Pi(1 <= Pi <= N)， 表示这N个数字的某一全排列。
输出格式：
输出一个整数，表示不同连号区间的数目。
示例：
用户输入：
4
3 2 4 1
程序应输出：
7
用户输入：
5
3 4 2 5 1
程序应输出：
9
解释：
第一个用例中，有7个连号区间分别是：[1,1], [1,2], [1,3], [1,4], [2,2], [3,3], [4,4]
第二个用例中，有9个连号区间分别是：[1,1], [1,2], [1,3], [1,4], [1,5], [2,2], [3,3], [4,4], [5,5]
/* 分析：
将所有的区间情况遍历出来，
然后判断在每个区间上，数字是不是连续的。
判断方法：求出区间上数字的最大值和最小值，
然后判断他们的差值是不是等于区间上元素个数-1。*/
#include <iostream>
using namespace std;
int main(){
int n;
cin>>n;
int a[n];
for(int i=0;i<n;i++){
cin>>a[i];
}
long long ans = 0;
int max,min;
for(int i=0;i<n-1;i++){
max = min = a[i];
for(int j=i+1;j<n;j++){
if(a[j]>max)
max = a[j];
if(a[j]<min)
min = a[j];
if((max-min)==(j-i))
ans++;
}
}
cout<<ans+n<<endl; //注意所有的单个数都是满足条件的，所以当第二层循环若是从i+1开始遍历，则最后结果要加上n。 如果第二层循环从i开始遍历，则不需要加n.
return 0;
}
2954
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
