A.数列
题目描述
下面数列的第 n 项：
(f(0)=a0,f(1)=a1,f(2)=a2)
(f(n)=b×f(n−1)+c×f(n−2)+d×f(n−3)+e(n≥3))
输入格式
包含 1 行，共 8 个整数：(a0、a1a2、b、c、d、e、n)
输出格式
输出 (f(n)) 的后 (18) 位（后 (18) 位的前缀 (0) 需要输出，不足 (18) 位用 (0) 补齐）。
样例输入
1 2 3 4 5 6 7 3
样例输出
000000000000000035
数据范围
对于 (30 \%) 的数据，(0 leq a_0,a_1,a_2,b,c,d,e,n≤10^6)
对于 $100 % $ 的数据，(0≤a_0,a_1,a_2,b,c,d,e,n≤10^{18})
这个题，考场上直接敲了暴力递推，yy能水30分，结果妥妥爆零
正解是用矩阵快速幂+高精（说好不考高精度，某虎说一时失误选了这个题）
代码除了难调没别的了。。。。
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define int long long
using namespace std;
inline int read(){
int x = 0, w = 1;
char ch = getchar();
for(; ch > '9' || ch < '0'; ch = getchar()) if(ch == '-') w = -1;
for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
return x * w;
}
const int maxn=1e6+50;
const int mod=1e18;
int a0,a1,a2,b,c,d,e,n;
int ans;
inline int Calc(int x,int y){
int a[30] = {0}, b[30] = {0}, c[30] = {0};
int lena = 0, lenb = 0;
while(x != 0 || lena == 0){
a[++lena] = x % 10;
x /= 10;
}
while(y != 0 || lenb == 0){
b[++lenb] = y % 10;
y /= 10;
}
for(int i = 1; i <= lena; i++){
for(int j = 1; j <= lenb; j++){
if(i + j - 1 <= 18){
c[i + j - 1] += a[i] * b[j];
c[i + j] += c[i + j - 1] / 10;
c[i + j - 1] %= 10;
}
}
}
int ans = 0;
for(int i = 18; i; i--){
ans = ans * 10 + c[i];
}
return ans;
}
struct Matrix{
int a[5][5];
void One(){
memset(a,0,sizeof(a));
for(int i=1;i<=4;i++)
a[i][i]=1;
}
inline void Chu(){
memset(a,0,sizeof(a));
a[1][1] = b;
a[1][2] = c;
a[1][3] = d;
a[1][4] = e;
a[2][1] = a[3][2] = a[4][4] = 1;
}
inline void Cheng(Matrix &b){
Matrix C;
memset(C.a,0,sizeof(C.a));
for(int i = 1; i <= 4; i++){
for(int j = 1; j <= 4; j++){
for(int k = 1; k <= 4; k++){
C.a[i][j] = C.a[i][j] + Calc(a[i][k], b.a[k][j]);
if(C.a[i][j] >= mod)
C.a[i][j] -= mod;
}
}
}
memcpy(a,C.a,sizeof(a));
}
};
Matrix qpow(int x){
Matrix last;
last.One();
Matrix now;
now.Chu();
while(x){
if(x & 1)
last.Cheng(now);
now.Cheng(now);
x >>= 1;
}
return last;
}
inline void Out(int x){
int b[30];
for(int i=0;i<18;i++){
b[i]=x%10;
x/=10;
}
for(int i=17;i>=0;i--){
printf("%lld",b[i]);
}
}
signed main(){
a0 = read(), a1 = read(), a2 = read(), b = read(), c = read(), d = read(), e = read(), n = read();
Matrix w = qpow(n - 2);
ans = Calc(w.a[1][1], a2) + Calc(w.a[1][2], a1) + Calc(w.a[1][3], a0) + w.a[1][4];
Out(ans);
return 0;
}
B.旗木双翼
题目描述
菲菲和牛牛在一块n行m列的棋盘上下棋，菲菲执黑棋先手，牛牛执白棋后手。
棋局开始时，棋盘上没有任何棋子，两人轮流在格子上落子，直到填满棋盘时结束。落子的规则是：一个格子可以落子当且仅当这个格子内没有棋子且这个格子的左侧及上方的所有格子内都有棋子。
_Itachi听说有不少学弟在省选现场AC了D1T1，解决了菲菲和牛牛的问题，但是_Itachi听说有的人认为复杂度玄学，_Itachi并不想难为学弟学妹，他想为大家节约时间做剩下的题，所以将简化版的D1T1带给大家。
_Itachi也在一块n行m列的棋盘上下棋，不幸的是_Itachi只有黑棋，不过幸好只有他一个人玩。现在，_Itachi想知道，一共有多少种可能的棋局（不考虑落子顺序，只考虑棋子位置）。
_Itachi也不会为难学弟学妹们去写高精度，所以只需要告诉_Itachi答案mod 998244353（一个质数）的结果。
输入格式
第一行包括两个整数n,m表示棋盘为n行m列。
输出格式
一个整数表示可能的棋局种数。
样例输入
10 10
样例输出
184756
考场上模拟出来杨辉三角，码了暴力，RE30，刚刚讲完发现很像一个卡特兰数，所以一个组合式就是答案了
对于除以n和m的阶乘，采用费马小定理，用快速幂处理(p - 2)次方就行了
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define int long long
using namespace std;
inline int read(){
int x = 0, w = 1;
char ch = getchar();
for(; ch > '9' || ch < '0'; ch = getchar()) if(ch == '-') w = -1;
for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
return x * w;
}
const int mod = 998244353;
inline int power(int a, int b){
int ans = 1;
while(b){
if(b & 1) ans = (ans % mod) * (a % mod) % mod;
a = (a % mod) * (a % mod) % mod;
b >>= 1;
}
ans %= mod;
return ans;
}
signed main(){
int n = read(), m = read();
int tmp = 1;
for(int i = m + n; i >= 1; i--){
tmp *= i;
tmp %= mod;
}
int n1 = 1;
int m1 = 1;
for(int i = 1; i <= n; i++){
n1 *= i;
n1 %= mod;
}
for(int i = 1; i <= m; i++){
m1 *= i;
m1 %= mod;
}
int ans = (tmp % mod) * power(n1 ,mod - 2) % mod * power(m1, mod - 2) % mod;
ans %= mod;
cout << ans << endl;
return 0;
}
C.乌龟棋
题目描述
小明过生日的时候，爸爸送给他一副乌龟棋当作礼物。
乌龟棋的棋盘是一行N个格子，每个格子上一个分数（非负整数）。棋盘第1格是唯一的起点，第N格是终点，游戏要求玩家控制一个乌龟棋子从起点出发走到终点。
(1 2 3 4 5……N)
乌龟棋中M张爬行卡片，分成4种不同的类型（M张卡片中不一定包含所有4种类型的卡片，见样例），每种类型的卡片上分别标有1、2、3、4四个数字之一，表示使用这种卡片后，乌龟棋子将向前爬行相应的格子数。游戏中，玩家每次需要从所有的爬行卡片中选择一张之前没有使用过的爬行卡片，控制乌龟棋子前进相应的格子数，每张卡片只能使用一次。
游戏中，乌龟棋子自动获得起点格子的分数，并且在后续的爬行中每到达一个格子，就得到该格子相应的分数。玩家最终游戏得分就是乌龟棋子从起点到终点过程中到过的所有格子的分数总和。
很明显，用不同的爬行卡片使用顺序会使得最终游戏的得分不同，小明想要找到一种卡片使用顺序使得最终游戏得分最多。
现在，告诉你棋盘上每个格子的分数和所有的爬行卡片，你能告诉小明，他最多能得到多少分吗？
输入格式
第1行2个正整数N和M，分别表示棋盘格子数和爬行卡片数。
第2行N个非负整数，a1, a2,……, an，其中ai表示棋盘第i个格子上的分数。
第3行M个整数，b1,b2,……, bm，表示M张爬行卡片上的数字。
输入数据保证到达终点时刚好用光M张爬行卡片，即N−1=所有卡片上的数字之和。
输出格式
输出只有1行，1个整数，表示小明最多能得到的分数。
样例输入
9 5
6 10 14 2 8 8 18 5 17
1 3 1 2 1
样例输出
73
样例1解释
小明使用爬行卡片顺序为1，1，3，1，2，得到的分数为6+10+14+8+18+17=73。注意，由于起点是1，所以自动获得第1格的分数6。
作为考场上唯一想出来并且 能够 A掉的题，爷忘了删调试了，忘了，，忘了，，，结果就是，从前十掉到了后十，自闭极了，就就就当是长教训了，以后注意吧，希望看我博客的筒子们引以为戒吧，哎
#include <iostream>
#include <cstring>
#include <cstdio>
#include <algorithm>
#define int long long
using namespace std;
inline int read(){
int x = 0, w = 1;
char ch = getchar();
for(; ch > '9' || ch < '0'; ch = getchar()) if(ch == '-') w = -1;
for(; ch >= '0' && ch <= '9'; ch = getchar()) x = x * 10 + ch - '0';
return x * w;
}
const int maxn = 52;
int dp[maxn][maxn][maxn][maxn];
int step[520];
int c[maxn];
signed main(){
int n = read(), m = read();
for(int i = 1; i <= n ;i++) step[i] = read();
for(int i = 1; i <= m; i++){
int u = read();
c[u]++;
}
dp[c[1]][c[2]][c[3]][c[4]] = step[1];
for(int i = c[1]; i >= 0; i--)
for(int j = c[2]; j >= 0; j--)
for(int l = c[3]; l >= 0; l--)
for(int r = c[4]; r >= 0; r--){
int u = 1 * (c[1] - i) + 2 * (c[2] - j) + 3 * (c[3] - l) + 4 * (c[4] - r) + 1;
dp[i][j][l][r] = max(dp[i][j][l][r], dp[i + 1][j][l][r] + step[u]);
dp[i][j][l][r] = max(dp[i][j][l][r], dp[i][j + 1][l][r] + step[u]);
dp[i][j][l][r] = max(dp[i][j][l][r], dp[i][j][l + 1][r] + step[u]);
dp[i][j][l][r] = max(dp[i][j][l][r], dp[i][j][l][r + 1] + step[u]);
}
/*
for(int i = 0; i <= c[1]; i++, cout << endl)
for(int j = 0; j <= c[2]; j++)
for(int l = 0; l <= c[3]; l++, cout << endl)
for(int r = 0; r <= c[4]; r++){
cout << "dp[";
cout << i;
cout << "][";
cout << j;
cout << "][";
cout << l;
cout << "][";
cout << r;
cout << "]=";
cout << dp[i][j][l][r] << " ";
}
*/
cout << dp[0][0][0][0] << endl;
return 0;
}
D. 假面舞会
题目描述
一年一度的假面舞会又开始了，栋栋也兴致勃勃的参加了今年的舞会。
今年的面具都是主办方特别定制的。每个参加舞会的人都可以在入场时选择一 个自己喜欢的面具。每个面具都有一个编号，主办方会把此编号告诉拿该面具的人。为了使舞会更有神秘感，主办方把面具分为k (k≥3)类，并使用特殊的技术将每个面具的编号标在了面具上，只有戴第i 类面具的人才能看到戴第i+1 类面具的人的编号，戴第k 类面具的人能看到戴第1 类面具的人的编号。 参加舞会的人并不知道有多少类面具，但是栋栋对此却特别好奇，他想自己算出有多少类面具，于是他开始在人群中收集信息。 栋栋收集的信息都是戴第几号面具的人看到了第几号面具的编号。如戴第2号面具的人看到了第5 号面具的编号。栋栋自己也会看到一些编号，他也会根据自己的面具编号把信息补充进去。由于并不是每个人都能记住自己所看到的全部编号，因此，栋栋收集的信 息不能保证其完整性。现在请你计算，按照栋栋目前得到的信息，至多和至少有多少类面具。由于主办方已经声明了k≥3，所以你必须将这条信息也考虑进去。
输入格式
第一行包含两个整数n, m，用一个空格分隔，n 表示主办方总共准备了多少个面具，m 表示栋栋收集了多少条信息。
接下来m 行，每行为两个用空格分开的整数a, b，表示戴第a 号面具的人看到了第b 号面具的编号。相同的数对a, b 在输入文件中可能出现多次。
输出格式
包含两个数，第一个数为最大可能的面具类数，第二个数为最小可能的面具类数。
如果无法将所有的面具分为至少3 类，使得这些信息都满足，则认为栋栋收集的信息有错误，输出两个-1。
样例输入
6 5
1 2
2 3
3 4
4 1
3 5
样例输出
4 4
考场想到了判环，写了一手spfa判了环然后想不出来咋用了，就弃掉了，然后输出-1水了10分
正解：如果有一个环，那么最长就是环的长度，最小就是最小的约数
对于多个环的情况，最长就是最大公约数，最少就是答案为不小于3的，同时是所有环中点的个数的最小的约数
对于全是链的情况，长度之和即为答案
#include <cstdio>
#include <iostream>
#include <cstring>
#include <algorithm>
using namespace std;
const int maxn = 1000000 + 10, maxm = 1000000 + 10;
#define ll long long
int head[maxn], len = -1;
struct node{
int nex, to, w;
}e[maxm << 1];
void Add(int u, int v, int w){
e[++len].to = v;
e[len].nex = head[u];
e[len].w = w;
head[u] = len;
}
int n, m, dep[maxn], maxx, minn, ans;
bool vis[maxn], flag[maxm << 1];
int gcd(int x, int y){
return y == 0 ? x : gcd(y, x % y);
}
void Dfs1(int u){
vis[u] = 1;
for(int i = head[u]; ~i; i = e[i].nex){
int v = e[i].to;
if (!vis[v]){
dep[v] = dep[u] + e[i].w;
Dfs1(v);
}
else {
ans = gcd(ans, abs(dep[u] + e[i].w - dep[v]));
}
}
}
void Dfs2(int u){
maxx = max(maxx, dep[u]);
minn = min(minn, dep[u]);
vis[u] = 1;
for(int i = head[u]; ~i; i = e[i].nex) {
if (!flag[i]){
flag[i] = flag[i ^ 1] = 1;
int v = e[i].to;
dep[v] = dep[u] + e[i].w;
Dfs2(v);
}
}
}
int main(){
cin >> n >> m;
memset(head, -1, sizeof head);
int u, v;
for(int i = 1; i <= m; i++){
scanf("%d %d", &u, &v);
Add(u, v, 1);
Add(v, u, -1);
}
for(int i = 1; i <= n; i++){
if(!vis[i]) Dfs1(i);
}
if(ans){
if(ans < 3) {printf("-1 -1
"); return 0;}
else{
int i = 3;
for(; i <= ans; i++) if(ans % i == 0) break;
printf("%d %d
", ans, i);
return 0;
}
}
memset(vis, 0, sizeof vis);
ans = 0;
for(int i = 1; i <= n; i++){
if(!vis[i]){
maxx = minn = dep[i] = 0;
Dfs2(i);
ans += maxx - minn + 1;
}
}
if(ans >= 3) printf("%d 3
", ans);
else printf("-1 -1
");
return 0;
}