123456

摘要： https://vjudge.net/problem/UVA-10870 裸的矩阵快速幂 注意系数矩阵在前面 因为系数矩阵为d*d 方程矩阵为d * 1 放反了就是d * 1 d * d 不符合矩阵乘法 #include<bits/stdc++.h> using namespace std; typ 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=3884 拓展欧拉定理 http://blog.csdn.net/Pedro_Lee/article/details/51458773这篇写的不错 我不会用latex。。。就不写了 #include<bi 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=4870 矩阵快速幂。。。 人话题意：从nk个物品里选模k余r个物品，问方案数模P 那么我们有方程 f[i][j]=f[i-1][j]+f[i-1][j-1] 跟组合数一个样子 j∈(0,k) 这个物品选 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=4868 三分+贪心 我们可以知道这是一个单峰函数 当A>B那么我们每次调整一个的价钱是最佳的，所以全部用B，如果A<B，那么我们尽量用A，不行再用B #include<bits/stdc++.h> us 阅读全文

该文被密码保护。 阅读全文

摘要： #include<bits/stdc++.h> using namespace std; int main() { int T; scanf("%d", &T); while(T--) { system("data > input.txt"); system("a+b < input.txt > 1 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=4518 水了一发斜率优化 好久没写了 调了一个晚上 dp[i][j]:第i天走到了第j个休息站 dp[i][j] = min(dp[i-1][x] + (s[i]-s[x]-v)^2)v是平均数=s[n 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=4810 问题就在于怎么快速查询 我们先用莫队转移，但是没办法快速地查询，那么我们就用bitset这个东西快速查询。 cnt是一个权值数组，记录每个数出现次数，a，c是两个bitset 减法：a[i]-a 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=4407 以前写过一次线性筛 发现不是很理解 写了个欧拉筛的 t了 其实每次推式子，都会先推出一组的解法，然后推出多组的 化简的方法大概是 先反演 然后把某个里层的sigma提出 然后有两种情况 1.复杂 阅读全文

摘要： http://www.lydsy.com/JudgeOnline/problem.php?id=3669 lct维护最小生成树 裸题 最小的边一定在最小生成树上 如果我们能用其他边调整 那么我们从能调整的边中选一条 因为肯定有一条比替换掉的小 那么就矛盾了 #include<bits/stdc++. 阅读全文