对拍
对拍是什么?
打个比方。比如你有一个不一定对的正解代码,和一个一定不会 WA 的错解代码,那么你可以造样例来确认正解代码有没有正确性。
对拍就是这个原理。在正常 OJ 刷题的过程中,如果很久调不出代码,可以复制一篇思路相同的题解过来,输出中间变量进行调试。
对拍的前提是你能够写出暴力代码。比如 dfs、bfs(开无限栈)、暴力 for 循环等。
过程
造数据可以完全随机生成。C++ 在 <stdlib>
头文件里面有专属的 rand()
函数可以返回一个随机值。
但是,rand()
函数需要初始化随机数种子,通常的说法是 srand(time(nullptr))
。这时,我们还需要引入 <ctime>
头文件。你也可以根据刷题的经验来发现一个极品种子。(据 CJ 同学所述,
114514
114514
114514 很容易被卡)。
rand()
函数的值在
[
1
,
2
16
−
1
]
[1,2^{16}-1]
[1,216−1] 之间,范围较小。本人通常多来几个 rand()
乘一起或者再来一个 rand()
作指数取模。
之后,顺着题目范围去跑,用 freopen()
放到一个 test.in
里面。(注意,可以写作 freopen("test.in","w",stdout);
)。之后,用暴力跑一份 test.ans
文件,正解跑一遍 test.out
文件。再弄一个程序 check.cpp
,分别从两个文件里面开 freopen("test.ans/out","r",stdin);
,最后 fclose(stdin);
来切换文件。把输出放在容器里面(比如 vector
),一一对比。
如果赶时间,可以只造一组(但是要很极品),如果时间充裕并且求稳,可以多造几组。如果没有差异,那么说明正解很大的概率是对的。反之,就把没过的数据扒出来看。
例子1
用 A+B Problem 做个示范。
比如,有暴力代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int a,b,ans=0;
scanf("%d %d",&a,&b);
while(a--){
++ans;
}
while(b--){
++ans;
}
printf("%d",ans);
return 0;
}
有正解代码:
#include<bits/stdc++.h>
using namespace std;
int main(){
int a,b;
scanf("%d %d",&a,&b);
printf("%d",a+b);
return 0;
}
造样例:
#include<bits/stdc++.h>
#define MOD 1000000001//不超过 1e9 范围
using namespace std;
typedef long long ll;
string a[5]={"1","2","3","4","5"};//假如造 5 组数据
inline ll power(ll x,ll y){//快速幂指数
ll res=1;
while(y){
if(y&1){
(res*=x)%=MOD;
}
(x*=x)%=MOD;
y>>=1;
}
return res;
}
inline ll randf(){//可能对多个 rand() 相乘比较有帮助
return rand();
}
inline ll randp(){
return power(randf(),randf());//指数
}
int main(){
for(int i=0;i<5;++i){
string inf="test"+a[i]+".in";//拼接
freopen(inf.c_str(),"w",stdout);
int a=randp(),b=randp();
printf("%d %d\n",a,b);
fclose(stdout);//关闭程序
}
return 0;
}
接着,用各种方法跑文件:
#include<bits/stdc++.h>
using namespace std;
string a[5]={"1","2","3","4","5"};
int pmain(){
int a,b,ans=0;
scanf("%d %d",&a,&b);
while(a--){
++ans;
}
while(b--){
++ans;
}
printf("%d",ans);
return 0;
}
int fmain(){
int a,b;
scanf("%d %d",&a,&b);
printf("%d",a+b);
return 0;
}
int main(){
for(int i=0;i<5;++i){
string inf="test"+a[i]+".in";
string ans="test"+a[i]+".ans";
string ouf="test"+a[i]+".out";
freopen(inf.c_str(),"r",stdin);
freopen(ans.c_str(),"w",stdout);
pmain();//跑暴力
fclose(stdout);
freopen(ouf.c_str(),"w",stdout);
fmain();//跑正解
fclose(stdin);
fclose(stdout);
}
}
用 check
比较文本:
#include<bits/stdc++.h>
using namespace std;
string a[5]={"1","2","3","4","5"};
int main(){
for(int i=0;i<5;++i){
string ans="test"+a[i]+".ans";
string ouf="test"+a[i]+".out";
vector<int> ans,ouf;
freopen(ans.c_str(),"r",stdin);
int n;
while(~scanf("%d",&n)){
ans.push_back(n);
}
fclose(stdin);
freopen(ouf.c_str(),"r",stdin);
int n;
while(~scanf("%d",&n)){
ouf.push_back(n);
}
fclose(stdin);
for(int i=0;i<ans.size();++i){
printf("%d %d\n",ans[i],ouf[i]);
}
}
}
数据点分治
这个用于防写正解全挂。
比如,你有一个全排列暴力 dfs
代码,能过
1
≤
n
≤
9
1\le n\le 9
1≤n≤9 的数据,然后还有一个正解代码,正确性未知。如果求稳,那就当
1
≤
n
≤
9
1\le n\le 9
1≤n≤9 的情况用 dfs
跑,别的地方用正解。就算正解挂了,也能够得到 dfs
的部分分。
本算法源于洛谷 @limingyuan333 巨佬在这场比赛的 T2 中,发现有一份代码没过 Subtask2,过了 Subtask3,有一份代码正好过了其他点。然后就用数据点分治 AC 了。
经过 JW 大佬们的不断演化,最终形成了数据点分治的算法。
规模分析法
规模分析法指的是根据数据规模进行估算时间复杂度多项式,来确定算法。
比如,输入一个 n ( 1 ≤ n ≤ 5 × 1 0 5 ) n(1\le n\le 5\times 10^5) n(1≤n≤5×105),并且输入若干个数,那么这个很可能是 n log n n\log n nlogn 算法或者 n n n\sqrt{n} nn 算法。比如二分、线段树、平衡树、分块等。
再比如一个 n ( 1 ≤ n ≤ 1 0 12 ) n(1\le n\le 10^{12}) n(1≤n≤1012),并且要求求出 n n n 的多项式,那么可能是数论分块的 n \sqrt{n} n 算法。如果还需要莫比乌斯反演或者欧拉反演并且求的 φ ( n ) \varphi(n) φ(n) 或者 μ ( n ) \mu(n) μ(n) 中的 n ( 1 ≤ n ≤ 2 31 − 1 ) n(1\le n\le 2^{31}-1) n(1≤n≤231−1),那么应该还需要杜教筛。
如果 1 ≤ n ≤ 1 0 6 1\le n\le 10^6 1≤n≤106,那么可能需要线性算法。
比如输入一个 n ( 1 ≤ n ≤ 1 0 18 ) n(1\le n\le 10^{18}) n(1≤n≤1018) 并且要求输出一个数,那么大概率是 O ( 1 ) \operatorname{O}(1) O(1) 算法。
卡常
CCF 不卡常!但是还是需要规避常数大的情况。
正常情况下,要规避不断地乘或者取模。
一般情况下,<cstdio>
比 <iostream>
快。但是有一些玄学情况,解绑 <iostream>
跑得飞快,快读却 T 了。(@limingyuan333 考试亲测)。
在 for
循环内,可以开 register
加速。常量可以加一个 const
或者 const&
,可以加速。(本人亲测,T 的代码用 const
立大功)。
通常,位运算比正常运算快。所以,n*2=n<<1
,n/2=n>>1
,n%2=n&1
。