第1题：统计不含4的数字
题目大意
统计10000至99999中，不包含4的数值个数。
题解
(@『落』常笑鹰 提供)直接推导出数学公式
最高位除了0、4不能使用，其余8个数字(1,2,3,5,6,7,8,9)均能使用，剩下的四位(千位、百位、十位、个位)可以使用除了4以外的所有数字，所以共有 8*9*9*9*9种解，计算得答案为：52488。
第2题：计算1千天后的日期
题目大意
2014-11-09再过1000天是哪一日？
题解
直接打开Excel计算：
(记得结果要按2017-08-05的格式提交)
扩展些Excel知识
Excel其实是用一个实数来存储日期，这个实数的整数部分，就是1900年1月1日到当前日期的间隔天数，所以日期与整数间，可以直接做加减运算。
有时候日期小数部分的值非0，小数部分就是对该天24小时的精确指定。
Excel是很强大的，闲着无聊的同学，也可以用下图的方式把日期改成题目要求的格式：
一些元老级的参赛选手，估计还做过“高斯日记”这题，那题也是可以用Excel来解：
扩展的扩展：Excel的日期计算其实是有个很小的bug的，这在Joel的《软件随想录》中有提到他在微软的经历，挺有趣的一段历史遗留问题，有兴趣的同学不妨看看他的书。
第3题：竖式加法
题目大意
祥 瑞 生 辉
三 羊 献 瑞
−−−−−−−−
三 羊 生 瑞 气
题目用了8个不同的汉字，表示0~9里八种不同的数字。组成两个数值相加，等于第三个数值。
题解
定义一个数组int a[10]，初始化用a[i]存储i。将a[2]~a[9]与各个汉字对应，然后用next_permutation全排列暴力解，注意三个数值开头不能为0，能解出第2个数值，即“三羊献瑞”对应的数字是1085。
#include #include using namespace std;
int main() {
int a[10];
for (int i = 0; i < 10; i++) a[i] = i;
do {
if (!a[2] || !a[6]) continue;
int x = a[2]*1000 + a[3]*100 + a[4]*10 + a[5];
int y = a[6]*1000 + a[7]*100 + a[8]*10 + a[3];
int z = a[6]*10000 + a[7]*1000 + a[4]*100 + a[3]*10 + a[9];
if (x + y == z) printf("%d + %d = %d\n", x, y, z);
} while (next_permutation(a, a+10));
return 0;
}
当然，这题也可以用数学知识简单的推导出一部分数字后，再暴力解剩下的部分。
第4题：古怪的星号修饰符
题目大意
是道代码填空题，主要是完成一个字符串s，按宽度width截断后，在固定宽度为width，两边为符号’|’的居中输出。
难点是题目给出了printf(“%*s%s%*s”,___)，要求填printf的实参列表。
题解
有些童鞋可能知道scanf里用*修饰符，是起到过滤读入的作用。比如一个有三列数值的数据，我只想得到第2列数值，可以在循环里用scanf(“%*d%d%*d”, a[i])来读入第i行的第2个数值到a[i]。
* 修饰符在printf中的含义完全不同。如果写成printf(“%6d”, 123)，很多童鞋应该就不会陌生了，这是设置域宽的意思。同理，%6s也是域宽。* 修饰符正是用来更灵活的控制域宽。使用%*s，表示这里的具体域宽值由后面的实参决定，如printf(“%*s”, 6, “abc”)就是把”abc”放到在域宽为6的空间中右对齐。
明白了 * 是用变量来控制域宽，那么这题就简单了，这里应该填写5个实参。然后字符长度的计算应该用buf而不是s，因为buf才是截断后的长度，用s的话，如果s长度超过了width-2，效果就不对了。
简单还原现场及参考答案：
#include #include int main() {
char s[100] = "abcd1234";
char buf[1000];
int width = 20;
strcpy(buf, s);
buf[width-2] = 0;
printf("|%*s%s%*s|\n", (width-strlen(buf)-2)/2, "", buf, (width-strlen(buf)-2)/2, "");
return 0;
}
第5题：补充全排列的回溯算法
题目大意
1,2,3…9 这九个数字组成一个分数，其值恰好为1/3，如何组法？
题解
与第3题一样可以用全排列暴力，但这题是代码填空题，且已经手写一部分全排列算法，我们负责补充其中一行代码。
写全排就是用回溯的思想，直接猜到for循环里的第三行，应该就是把第一行的交换操作再交换回来~~复制for里的第一行代码，运行下程序，还不放心就调试下，看看数组是不是有按字典序在变化就行了。
答案：{t=x[k]; x[k]=x[i]; x[i]=t;}
第6题：加号改乘号
题目大意
把1+2+3+…+48+49中的两个加号改成乘号(修改位置不能相邻)，使得式子的结果由1225变为2015。
题解
用双循环暴力两个乘号的位置，计算在数字i、j后的加号改为乘号，式子数值的变化即可，注意j的起始位置为i+2。
#include int main() {
for (int i = 1; i <= 48; i++) {
for (int j = i + 2; j <= 48; j++) {
if (1225 - i - (i+1) - j - (j+1) ==
2015 - i*(i+1) - j*(j+1))
printf("%d\n", i);
}
}
return 0;
}
除了输出的一个10是题目已提示的解，另一个16就是答案。
第7题：牌型种数
题目大意
去掉大小王的52张扑克牌，不考虑花色及顺序，取13张，一共有多少种取法？
题解
大部分人应该都能想到用13重循环暴力解(每层循环都是0~4)的方法，当累加和为13时，ans++。
我在赛中想到的是动态规划，因为我想到了“小明爬楼梯”，且跟赛前，__M子__给我们讲的去年赛题“地宫取宝”的dp思路非常像。
假设牌是从1到13按顺序取的，dp[i][j]表示取到第i号的牌，目前总共有j张牌的取法总数，那么有状态转移方程(注意公式没考虑边界处理)：
dp[i][j]=∑k=j−4jdp[i−1][k]
#include using namespace std;
typedef long long LL;
LL dp[14][14]; // dp[i][j]: 当前到第i张牌，总共有j张牌时的解的个数
int main() {
dp[1][0] = dp[1][1] = dp[1][2] = dp[1][3] = dp[1][4] = 1;
for (int i = 2; i <= 13; i++) {
for (int k = 0; k <= 13; k++) {
if (k - 4 >= 0) dp[i][k] += dp[i-1][k-4];
if (k - 3 >= 0) dp[i][k] += dp[i-1][k-3];
if (k - 2 >= 0) dp[i][k] += dp[i-1][k-2];
if (k - 1 >= 0) dp[i][k] += dp[i-1][k-1];
dp[i][k] += dp[i-1][k];
}
}
cout << dp[13][13] << endl;
return 0;
} 答案：3598180
第8题：计算房子间曼哈顿距离
题目大意
房子按S形摆放，如
1 2 3
6 5 4
7 8 ……
现输入每行的宽度w，计算出任意两个房子m、n的曼哈顿距离(横向距离+竖向距离)。
题解
直接看代码
#include #include using namespace std;
// 输入宽度w和房子编号n，返回房子所在行x，列y
void GetPos(int w, int n, int& x, int& y) {
x = (n-1) / w + 1;
y = n % w;
if (y == 0) y = w;
if (x%2 == 0) { // 偶数行要倒着数
y = w - y + 1;
}
}
int main() {
int w, m, n;
int x1, y1, x2, y2;
cin >> w >> m >> n;
GetPos(w, m, x1, y1);
GetPos(w, n, x2, y2);
cout << abs(x1-x2) + abs(y1-y2) << endl;
return 0;
}
第9题：垒骰子(矩阵快速幂模)
题目大意
输入n,m，表示用n个骰子，在m个约束条件下，从下往上叠成一列。接下来m行，每行有两个数a,b，表示骰子之间数字a、b两个所在面不能拼接。
问这n个骰子一共有多少种垒法？答案对1e9+7取余。
注意骰子摆好后，四个方向均能转动，即样例输入：
2 1
1 2
的样例输出是：
544
矩阵快速幂基础
请读者先看这篇文章：hdu5171 GTY’s birthday gift(BestCoder Round #29 1002)，或阅读其他矩阵快速幂模文章，自行补基础。我给出的这篇文章，写了个矩阵快速幂模的类模板，将在下面的代码中使用。实际解题中，可以根据特殊问题写特殊代码，不必写出通用的矩阵类，加快解题速度。
题解
比赛中，不使用模板，建议数组从下标1开始使用。这里为了和模板配合，所以假设骰子从0开始编号。可以发现规律，如果骰子的a面为0,1,2，那么a的对面数值就是a+3；相反，如果骰子的a面为3,4,5，那么a的对面数值就是a-3。这个由当前面变换到对面的操作可以用公式完成：(a+3)%6。
现在用一个bool数组isLimit[6][6]来存储，isLimit[i][j]为真，代表当前骰子i面朝上，不能接j面朝上的骰子。
矩阵快速幂模，主要是找到A^(n-1) * X中的A和X。那么这里的X向量，X[i]初值就是第一个骰子i面朝上的方法数，一共是4种。即X[0]~X[5]均为4。
A[i][j]就是当前骰子i面朝上，接j面朝上后，原方法数要乘上的因子。故A[i][j]可以由isLimit直接得到。如果i->j被限制，乘因子就是0，否则就是4。
#include #include #include using namespace std;
/*矩阵类(用于快速幂模)
1、T是整型类型，N是方阵大小，MOD是取余的值
2、Eye是单位阵
*/
template class Matrix {
T val[N][N];
public:
Matrix() { memset(val, 0, sizeof(val)); }
Matrix(T a[N][N]) { memcpy(val, a, sizeof(val)); }
Matrix operator*(const Matrix& c) const {
Matrix res;
for (int i = 0; i < N; ++i)
for (int j = 0; j < N; ++j)
for (int k = 0; k < N; ++k) {
res.val[i][j] += val[i][k] * c.val[k][j];
//防止矩阵元素变为负数，若不需要，去掉"+MOD"
res.val[i][j] = (res.val[i][j] + MOD) % MOD;
}
return res;
}
Matrix& operator*=(const Matrix& c) {
*this = *this * c;
return *this;
}
Matrix operator^(int k) const { //返回*this^k
Matrix res = Eye();
Matrix step(*this);
while (k) {
if (k & 1) res *= step;
k >>= 1;
step *= step;
}
return res;
}
Matrix Eye() const {
Matrix a;
for (int i = 0; i < N; i++) a.val[i][i] = 1;
return a;
}
void out() const {
for (int i = 0; i < N; i++) {
for (int j = 0; j < N; j++)
cout << val[i][j] << " ";
cout << "\n";
}
}
T* operator[](int i) { return val[i]; }
};
typedef long long LL;
const int MOD = 1e9 + 7;
int main() {
int n, m;
MatrixA, X, Y;
bool isLimit[6][6];
// 输入处理
cin >> n >> m;
memset(isLimit, false, sizeof(isLimit));
for (int i = 0; i < m; i++) {
int a, b;
cin >> a >> b;
isLimit[a-1][ (b+2)%6 ] = true;
isLimit[b-1][ (a+2)%6 ] = true;
}
// 初始化A
for (int i = 0; i < 6; i++) {
for (int j = 0; j < 6; j++) {
A[i][j] = (!isLimit[i][j]) * 4;
}
}
// 初始化X
for (int i = 0; i < 6; i++) X[i][0] = 4;
// 计算Y及答案
Y = (A^(n-1)) * X;
LL sum = 0;
for (int i = 0; i < 6; i++) sum += Y[i][0];
cout << sum % MOD << endl;
return 0;
}
第10题：树形DP
题目大意
在一颗树中，每个点都有一个权值，有正有负。
要求找一个连通的点集S，使得累加和最大。
题解
树形DP不是我熟悉的模块，这个坑就交给别人填了
。
树形DP参考博客：“树形DP(HDOJ1011 2196 4003 5148 POJ2342)”，及下方 __M子__提供的代码：
#include#include#include#includeusing namespace std;
typedef long long LL;
const int MAXN=100000+10;
const int INF=-1000000-10;
LL dp[MAXN][2];
vectortree[MAXN];
LL w[MAXN];
void dfs(int p,int fa){
dp[p][0]=max(dp[p][0],w[p]);
dp[p][1]=max(dp[p][1],w[p]);
for(int i=0;i0)dp[p][1]+=dp[son][1];
}
}
dp[p][0]=max(dp[p][0],dp[p][1]);
}
int main(){
int n,a,b;
cin>>n;
for(int i=1;i<=n;i++)cin>>w[i];
for(int i=1;i<=n;i++){dp[i][0]=INF;dp[i][1]=INF;}
for(int i=1;i>a>>b;
tree[a].push_back(b);
tree[b].push_back(a);
}
dfs(1,-1);
cout<
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
