题目链接 题意： 思路： #include <bits/stdc++.h> using namespace std; #define LL long long const int N = 1e7; const int mod = 20101009; int n, m, mu[N + 5], p[N Read More
\n与endl都有换行的含义，但有所区别。 对于cout输出流来说，并不是简单的输出在屏幕那么简单，而是先将后面接着的字符串存在一个缓冲槽。而对于\n来说，属于其中的字符串。对于endl来说，他会控制缓冲槽刷新清空缓冲槽，并且换行，相当于 "\n"<<std::flush; 那cout什么时候会自动 Read More
题目链接 题意： 思路：首先由容斥定理可得 红色箭头为反演步骤。 最后由数论分块来写。 #include <bits/stdc++.h> using namespace std; #define ll long long const int N = 50000; int mu[N + 5], p[N Read More
题目链接 题意：求∑k mod i 思路：可以化解为 所以，如何求k/i向下取整才是关键，对于一个块k/i=j，我们可以直接求i到j区间的值，i到j区间都为k/i，时间复杂度为sqrt(n)。大体如下证明方法可行性和复杂度 #include <bits/stdc++.h> using namespa Read More
题目链接 题意：求斐波那契第n项，只不过是最后m位。 思路：矩阵快速幂板子 #include <bits/stdc++.h> using namespace std; #define ll long long const int N=52; int mod=1000; struct Matrix { Read More
题目链接 题意：给定一个序列的变换，求变换r次后各项的值。 思路：矩阵快速幂+矩阵变换知识。 #include <bits/stdc++.h> using namespace std; #define ll long long const int N=52; const int mod=1000; Read More
题目链接 题意：给一个A矩阵大小为n*K，B矩阵大小为K*n。求（A*B）^(n*n)的矩阵的所有元素和。n<=1000,K<=6。 思路：首先结构体中最大开800*800，如果直接算就无法运行。所以要想到式子转化为A*(B*A)^(n*n-1)*B。这样只要最多开6*6了 剩下的就都是手写矩阵乘法 Read More
题目链接 题意：给定进制base，和分数score，求在base进制下，有多少个数的值为score，要求不能有连续相同的数字以及前导0.计算一个数的值即为相邻两位数的差平方的和。 思路：因为score很大，所以直接dp肯定超时，但是即使对于base=6的情况，每次新添一个数score最大增加25（0 Read More
题目链接 题意：求出代码中的值。 思路：两种：第一种打表发现fi=fi-1+2*fi-2+1 第二种： //a[n] = a[n-1] + n%2//a[n] = a[n-1] + (n+1)/2 - n/2 #include <bits/stdc++.h> #include<unordered_m Read More
题目链接 题意：第一行矩阵为23,233,2333...给出第一列的矩阵值，求an,m。 思路：很容易想到矩阵求解 大概就是这样中间一堆都是1. #include <bits/stdc++.h> using namespace std; #define ll long long const int Read More