考前整理了一下模板,供自己复习。(虽说考试已经过去了,不过明年你还可以用啊,已退役的OIer当我没说)
考前模板汇总
快速幂
#define ll long long
#include<bits/stdc++.h>
using namespace std;
ll power(ll a,ll b,ll c) //a为底数,b为指数,c 为模数
{
ll ans=1%c;
while(b)
{
if(b&1)
ans=ans*a%c;
b>>=1;
a=a*a%c;
}
return ans%c;
}
矩阵快速幂
#include<iostream>
#include<cstring>
#define mod 1000000007
#define ll long long
using namespace std;
struct Mat{
ll m[101][101];
};//结构体存矩阵
Mat a,e;//a是输入的矩阵,e是单位矩阵
ll n,p;
Mat Mul(Mat x,Mat y) //矩阵乘
{
Mat c;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
c.m[i][j]=0;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
for(int k=1;k<=n;k++)
{
c.m[i][j]=c.m[i][j]%mod+x.m[i][k]*y.m[k][j]%mod;
}
return c;
}
Mat pow(Mat x,ll y) //矩阵快速幂
{
Mat ans=e;
while(y)
{
if(y&1)
ans=Mul(ans,x);
x=Mul(x,x);
y>>=1;
}
return ans;
}
int main()
{
cin>>n>>p;
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
cin>>a.m[i][j];
for(int i=1;i<=n;i++)
e.m[i][i]=1;
Mat ans=pow(a,p);
for(int i=1;i<=n;i++)
{
for(int j=1;j<=n;j++)
cout<<ans.m[i][j]%mod<<" ";
cout<<endl;
}
return 0;
}
扩展欧几里得
#define ll long long
#include<bits/stdc++.h>
using namespace std;
ll exgcd(ll a,ll b,ll &x,ll &y)//x为逆元,ans为最大公因数
{
if(!b)
{
x=1;
y=0;
return a;
}
ll ans=exgcd(b,a%b,y,x);
y-=a/b*x;
return ans;
}
逆元递推
const int maxn=3e6+5;
ll n,p;//n为要求的数,p为模数
ll ivp[maxn];
int main()
{
scanf("%lld%lld",&n,&p);
ivp[0]=0;ivp[1]=1;
for( int i=2;i<=n;i++)
ivp[i]=p-p/i*ivp[p%i]%p;
printf("%lld\n",ivp[n]);
}
//逆元三种求法,快速幂,扩欧,递推
欧拉筛
#include<bits/stdc++.h>
#define re register
#define maxx 10000001
using namespace std;
bool b[maxx];
int n;
int prime[1000001],cnt;
void getprime(int n)
{
b[0]=1;
b[1]=1;
for(re int i=2;i<=n;i++)
{
if(!b[i])
prime[++cnt]=i;
for(re int j=1;j<=cnt&&prime[j]*i<=n;j++)
{
b[prime[j]*i]=1;
if(i%prime[j]==0) break;
}
}
}
void putprime()
{
for(re int i=1;i<=cnt;i++)
printf("%d ",prime[i]);
}
并查集
#define maxx 10000001
int f[maxx];
//查找
int find(int root)
{
if(f[root]==root) return root;
return f[root]=find(f[root]);
}
//合并
inline void uni(int x,int y)
{
f[find(x)]=find(y);
}
前缀和及差分
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e5+9;
int a[maxn],b[maxn];
int main(){
int i,j,k,n,m,p;
cin>>n>>m;
for(i=1;i<=n;i++){
cin>>a[i];
}
for(i=1;i<=m;i++){
int L,R,t;
cin>>t>>L>>R>>p;
if(t==1){
b[L]+=p;b[R+1]-=p;
}
else{
b[L]-=p;b[R+1]+=p;
}
}
int add=0;
for(i=1;i<=n;i++){
add+=b[i];
a[i]+=a[i-1]+add;
}
int x,y;
cin>>x>>y;
cout<<a[y]-a[x-1]<<endl;
}
二维前缀和及差分
#include<bits/stdc++.h>
using namespace std;
const int maxn=1e3+9;
int a[maxn][maxn];
int main(){
int i,j,k,n,m,q;
cin>>n>>m>>q;
for(i=1;i<=n;i++){
for(j=1;j<=m;j++)
cin>>a[i][j];
}
for(i=1;i<=n;i++){
for(j=1;j<=m;j++)
a[i][j]+=a[i][j-1]+a[i-1][j]-a[i-1][j-1];
}
for(i=1;i<=q;i++){
int x1,y1,x2,y2;
cin>>x1>>y1>>x2>>y2;
int ans=a[x2][y2]-a[x1-1][y2]-a[x2][y1-1]+a[x1-1][y1-1];
cout<<ans<<endl;
}
}
//二维差分
for(int i=0;i<m;i++){//m是修改操作次数
int x1,y1,x2,y2,p;
cin>>x1>>y1>>x2>>y2>>p;
b[x1][y1]+=p;b[x2+1][y2+1]+=p;
b[x2+1][y1]-=p;b[x1][y2+1]-=p;
}
组合数递推
C[1][0] = C[1][1] = 1;
for (int i = 2; i < N; i++){
C[i][0] = 1;
for (int j = 1; j < N; j++)
C[i][j] = (C[i - 1][j] + C[i - 1][j - 1])%mod;
}
卡特兰数(Catalan数)
h[0]=1;
for(int i=1;i<=n;i++)
h[i]=h[i-1]*(4*i-2)/(i+1);
//括号化:矩阵连乘:P=a1*a2*a3*…*an,依据乘法结合律,试问有几种括号化的方案?
//出栈次序
//买票找零:有2n个人排成一行进入剧场,入场费5元,其中只有n个人有一张5元钞票,另外n人只有10元钞票,剧院无其它钞票,问有多少种方法使得只要有10元的人买票,售票处就有5元的钞票找零?
//凸多边形三角划分:在一个凸多边形中,通过若干条互不相交的对角线,把这个多边形划分成了若干个三角形。
//一位大城市的律师在她住所以北n个街区和以东n个街区处工作。每天她走2n个街区去上班。如果她从不穿越(但可以碰到)从家到办公室的对角线,那么有多少条可能的道路?
//在圆上选择2n个点,将这些点成对连接起来使得所得到的n条线段不相交的方法数?
//给定节点组成二叉搜索树:给定N个节点,能构成多少种不同的二叉搜索树?
//n对括号正确匹配数目:给定n对括号,求括号正确配对的字符串数
最小生成树
#include<bits/stdc++.h>
#define maxx 10000001
using namespace std;
int f[maxx];
int find(int root) {
if(f[root]==root) return root;
return f[root]=find(f[root]);
}
struct node {
int l,r,c;
} bian[200001];
bool cmp(const node x,const node y) {
return x.c<y.c;
}
int n,m,cnt,ans;
int main() {
cin>>n>>m;
for(int i=1; i<=m; i++)
cin>>bian[i].l>>bian[i].r>>bian[i].c;
sort(bian+1,bian+m+1,cmp);
for(int i=1; i<=n; i++)
f[i]=i;
for(int i=1; i<=m; i++) {
if(find(bian[i].l)!=find(bian[i].r)) {
f[find(bian[i].l)]=find(bian[i].r);
ans+=bian[i].c;
cnt++;
if(cnt==n-1) {
cout<<ans;
return 0;
}
}
}
cout<<"orz";
return 0;
}
图论最短路
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
using namespace std;
struct node {
int to,dis,nxt;
} edg[500002];
int mod;
const int inf=2147483647;
bool vis[10002];
int n,m,s,head[10002],sum[10002],num,a,b,c;
int dis[10002];
inline void add(int x,int y,int z) {
edg[++num].to=y;
edg[num].dis=z;
edg[num].nxt=head[x];
head[x]=num;
}
//dijkstra算法
/*void dij() {
fill(dis+1,dis+n+1,inf);
dis[s]=0;
for(int i=1; i<=n; i++) {
int cur;
int mm=inf;
for(int j=1; j<=n; j++)
if(!vis[j]&&mm>dis[j])
mm=dis[cur=j];
if(mm==inf) break;
vis[cur]=1;
for(int p=head[cur]; p; p=edg[p].nxt) {
if(dis[edg[p].to]>mm+edg[p].dis&&!vis[edg[p].to])
dis[edg[p].to]=mm+edg[p].dis;
}
}
}*/
//spfa算法
/*
void spfa()
{
fill(dis+1,dis+n+1,inf);
dis[s]=0;
vis[s]=1;
q.push(s);
while(!q.empty())
{
int cur=q.front();
vis[cur]=0;
for(int i=head[cur];i;i=edg[i].nxt)
{
if(dis[edg[i].to]>dis[cur]+edg[i].dis)
{
dis[edg[i].to]=dis[cur]+edg[i].dis;
if(!vis[edg[i].to])
{
q.push(edg[i].to);
vis[edg[i].to]=1;
//sum[edg[i].to]=sum[cur];
}
}
//else if(dis[edg[i].to]==dis[cur]+edg[i].dis){
// sum[edg[i].to]+=sum[cur];
// sum[edg[i].to]%=mod;
//}
}
q.pop();
}
}*/
//dijkstra+堆优化
struct nod {
int now,x;
bool operator < (const nod &tmp) const {
return x>=tmp.x;
}
} cur;
priority_queue<nod>q;
void _dij() {
fill(dis+1,dis+n+1,inf);
dis[s]=0;
q.push((nod) {
s,dis[s]
});
while(!q.empty()) {
cur=q.top();
q.pop();
if(vis[cur.now]) continue;
vis[cur.now]=1;
if(dis[cur.now]!=cur.x) continue;
for(int p=head[cur.now]; p; p=edg[p].nxt) {
if(dis[edg[p].to]>dis[cur.now]+edg[p].dis) {
dis[edg[p].to]=dis[cur.now]+edg[p].dis;
q.push((nod) {
edg[p].to,
dis[edg[p].to]
});
//sum[edg[i].to]=sum[cur];
}
}
//else if(dis[edg[i].to]==dis[cur]+edg[i].dis){
// sum[edg[i].to]+=sum[cur];
// sum[edg[i].to]%=mod;
//}
}
}
}
int main() {
scanf("%d%d%d",&n,&m,&s);
for(int i=1; i<=m; i++) {
scanf("%d%d%d",&a,&b,&c);
add(a,b,c);
}
sum[s]=1;
_dij();//spfa();//dij();
for(int i=1; i<=n; i++)
cout<<dis[i]<<' ';
return 0;
}
图论大集合(数量,路径,距离)
#include<iostream>
#include<cstdio>
#include<queue>
#include<algorithm>
#include<stack>
using namespace std;
struct node {
int to,dis,nxt;
} edg[500002];
const int inf=2147483647;
bool vis[10002];
int n,m,s,head[10002],num,sum[10002],a,b,c;
int fa[100002]= {-1};
int dis[10002];
inline void add(int x,int y,int z) {
edg[++num].to=y;
edg[num].dis=z;
edg[num].nxt=head[x];
head[x]=num;
}
queue<int>q;
void spfa() {
fill(dis+1,dis+n+1,inf);
dis[s]=0;
vis[s]=1;
q.push(s);
while(!q.empty()) {
int cur=q.front();
vis[cur]=0;
for(int i=head[cur]; i; i=edg[i].nxt) {
if(dis[edg[i].to]>dis[cur]+edg[i].dis) {
dis[edg[i].to]=dis[cur]+edg[i].dis;
fa[edg[i].to]=cur;
if(!vis[edg[i].to]) {
q.push(edg[i].to);
vis[edg[i].to]=1;
sum[edg[i].to]=sum[cur];}
}
else if(dis[edg[i].to]==dis[cur]+edg[i].dis)
{
sum[edg[i].to]+=sum[cur];
}
}
q.pop();
}
}
int main() {
scanf("%d%d%d",&n,&m,&s);
for(int i=1; i<=m; i++) {
scanf("%d%d%d",&a,&b,&c);
add(a,b,c);
}
sum[s]=1;
spfa();
stack<int>k ;
for(int i=1;i<=n;i++)
{
cout<<"num:"<<sum[i]<<" "<<"path:";
for(int j=i;j!=s;j=fa[j])
k.push(fa[j]);
while(!k.empty())
{
cout<<k.top()<<">-";
k.pop();
}
cout<<i<<' ';
cout<<"dis:"<<dis[i]<<endl;
}
return 0;
}
高精加减乘
#include<bits/stdc++.h>
using namespace std;
char a1[600],b1[600],tmp[600];
int a[600],b[600],c[600],d[1201],e[600];
void jia(char x[],char y[]) {
int lena=strlen(x),lenb=strlen(y);
int z=0;
for(int i=0; i<lena; i++) a[lena-i]=x[i]-'0'; //倒位存进数组
for(int i=0; i<lenb; i++) b[lenb-i]=y[i]-'0'; //倒位存进数组
int lenc=1;
while(lenc<=lena||lenc<=lenb) {
c[lenc]=a[lenc]+b[lenc]+z;//加上进位
z=c[lenc]/10;//更新进位
c[lenc]%=10;//更新位数使其小于10
lenc++;
}
c[lenc]=z ;
if(c[lenc]==0) lenc--;//删除前导零
for(int i=lenc; i>=1; i--)
printf("%d",c[i]);//倒位输出
}
void chen(char x[],char y[]) {
int lena=strlen(x),lenb=strlen(y);
int z=0;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
for(int i=0; i<lena; i++) a[lena-i]=x[i]-'0';
for(int i=0; i<lenb; i++) b[lenb-i]=y[i]-'0';
for(int i=1; i<=lena; i++) {
z=0;
for(int j=1; j<=lenb; j++) {
d[i+j-1]=a[i]*b[j]+z+d[i+j-1];//当前乘积+上次乘积进位+原数
z=d[i+j-1]/10;//进位
d[i+j-1]%=10;//更新
}
d[i+lenb]=z;
}
int lend=lena+lenb;
while(d[lend]==0&&lend>1)//删除前导零
lend--;
for(int i=lend; i>=1; i--)
cout<<d[i];
}
void jian(char x[],char y[]) {
if(strlen(x)<strlen(y)||(strlen(x)==strlen(y)&&strcmp(x,y)<0))
{
strcpy(tmp,x);
strcpy(x,y);
strcpy(y,tmp);//交换x,y
cout<<"-";
}
int lena=strlen(x);
int lenb=strlen(y);
for(int i=0;i<lena;i++)
a[lena-i]=x[i]-'0';
for(int i=0;i<lenb;i++)
b[lenb-i]=y[i]-'0';
int lene=1;
int z=0;//是否借位
while(lene<=lena||lene<=lenb)
{
e[lene]=a[lene]-b[lene]-z;
if(e[lene]<0){
z=1;e[lene]+=10;}
else
z=0;
lene++;
}
e[lene]=z;
while(e[lene]==0&&lene>1)
lene--;
for(int i=lene;i>=1;i--)
printf("%d",e[i]);
}
int main() {
cin>>a1;
cin>>b1;
jian(a1,b1);
cout<<endl;
jia(a1,b1);
cout<<endl;
chen(a1,b1);
cout<<endl;
}
最近公共祖先(LCA)Tarjan
struct edge
{
int to,next;
}e[1000001];//边的储存
struct questions
{
int to,next,same,num;
bool flag;
questions(){flag=false;}
}q[1000001];//询问的储存,flag=false表示还没回答,num表示是第几个询问,same储存与这个询问相同的询问序号。
bool b[500001];
int head[500001],que[500001],father[500001];
int n,m,s,nume=0,numq=0,ans[500001];
void add_edge(int x,int y)
{
e[++nume].to=y;
e[nume].next=head[x];
head[x]=nume;
e[++nume].to=x;
e[nume].next=head[y];
head[y]=nume;
}
void add_que(int x,int y,int k)
{
q[++numq].to=y;
q[numq].same=numq+1;
q[numq].next=que[x];
q[numq].num=k;
que[x]=numq;
q[++numq].to=x;//询问要储存到两个点的链表序列里,删的时候也要一起删
q[numq].same=numq-1;
q[numq].next=que[y];
q[numq].num=k;
que[y]=numq;
}
int find(int x)//并查集
{
if(father[x]!=x)father[x]=find(father[x]);
return father[x];
}
void unionn(int x,int y)//并查集
{
father[find(y)]=find(x);
}
void LCA(int point,int f)//point是当前搜索节点,f是它的父亲
{
for(int i=head[point];i!=0;i=e[i].next)//遍历与point相连的所有边
if(e[i].to!=f&&!b[e[i].to])
{
LCA(e[i].to,point);
unionn(point,e[i].to);//合并
b[e[i].to]=1;
}
for(int i=que[point];i!=0;i=q[i].next)//遍历与point相关的询问
if(!q[i].flag&&b[q[i].to])//如果另一个点遍历过了并且该询问没有回答过
{
ans[q[i].num]=find(q[i].to);//记录下答案
q[i].flag=1;
q[q[i].same].flag=1;//把两个点上的询问都去掉
}
}
int main()
{
read(n);read(m);read(s);
for(int i=1,x,y;i<=n-1;i++)
{
father[i]=i;
read(x);read(y);
add_edge(x,y);
}
father[n]=n;
for(int i=1,x,y;i<=m;i++)
{
read(x);read(y);
add_que(x,y,i);
}
LCA(s,0);
for(int i=1;i<=m;i++)
printf("%d\n",ans[i]);
return 0;
}
字符串转化(插入,删除,替换)
for (int i=1; i<=row; i++)
dp[i][0] = i;//从长度为i的字符串到空串需要变换i次
for (int j=1; j<=col; j++)
dp[0][j] = j;//从长度为kong的字符串到长度为j的字符串需要变换j次
for (int i=1; i<=row; i++){
for (int j=1; j<=col; j++){
if (word1[i-1] == word2[j-1])
dp[i][j] = dp[i-1][j-1];
else
dp[i][j] = min(min(dp[i-1][j-1]+1,dp[i-1][j] +1), dp[i][j-1]+1 );
}
}
return dp[row][col];
}
最长公共子序列
for (int i = 0; i <= len1; i++) {
for( int j = 0; j <= len2; j++) {
if(i == 0 || j == 0) {
c[i][j] = 0;
} else if (str1.charAt(i-1) == str2.charAt(j-1)) {
c[i][j] = c[i-1][j-1] + 1;
} else {
c[i][j] = max(c[i - 1][j], c[i][j - 1]);
}
}
}
cout<<c[len1][len2];
子串
int result = 0; //记录最长公共子串长度
int c[][] = new int[len1+1][len2+1];
for (int i = 0; i <= len1; i++) {
for( int j = 0; j <= len2; j++) {
if(i == 0 || j == 0) {
c[i][j] = 0;
} else if (str1.charAt(i-1) == str2.charAt(j-1)) {
c[i][j] = c[i-1][j-1] + 1;
result = max(c[i][j], result);
} else {
c[i][j] = 0;
}
}
}
cout<< result;
}
最长上升子序列
int main()
{
int n;
while(cin>>n)
{
for(int i=0; i<n; i++)
{
cin>>a[i];
dp[i]=1; // 初始化为1,长度最短为自身
}
int ans=1;
for(int i=1; i<n; i++)//枚举子序列的终点
{
for(int j=0; j<i; j++)//从头向终点检查每一个元素
{
if(a[i]>a[j])
{
dp[i]=max(dp[i],dp[j]+1); // 状态转移
}
}
ans=max(ans,dp[i]); // 比较每一个dp[i],最大值为答案
}
cout<<ans<<endl;
}
最长上升公共子序列
for(int i=1;i<n;i++)
{
int mm=0;
for(int j=1;j<m;j++)
{
if(a[i]>b[j])mm=max(mm,f1[j]);//优化
else if(a[i]==b[j])
f1[j]=mm+1,ans2=max(ans2,f1[j]);
}
}
01背包
for(int i=1;i<=m;i++)
for(int j=1;j<=T;j++)
{
f[i][j]=f[i-1][j];
if(j>=t[i]) f[i][j]=max(f[i-1][j-t[i]]+v[i],f[i][j]);
}
printf("%d",f[m][T]);
常数优化及注意事项
//位运算
x*10 <=> (x<<3)+(x<<1)
x!=y <=> x^y
x!=-1 <=> ~x
x*2 <=> x<<1
x*2+1 <=> x<<1|1
x/2 <=> x>>1
(x+1)%2 <=> x^1
x%2 <=> x&1
x%2==0 <=> ~(x&1)
//语法
inline 在非递归函数前加修饰。
循环变量 int i =>register int i
x+=add;x%=mod;
=> x+=add;x>=mod?x%=mod:1;
三目运算符比if,else快 A?B:C;
//注意事项
有返回值的函数要加return A;
回溯加return ;
using namespace std;
return 0;要注意
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
仔细检查
if while 后不要加分号;
输入输出优化
int read()
{
int x=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-') f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=x*10+ch-'0';
ch=getchar();
}
return x*f;
}
void write(int x)
{
if(x<0)
{
putchar('-');
x=-x;
}
if(x>9) write(x/10);
putchar(x%10+'0');
}
树状数组
#include <iostream>
#include <cstdio>
#define lowbit(x) x & -x
using namespace std;
long long tree[500005];
int n, m;
void add(int x, long long num) {
while (x <= n) {
tree[x] += num;
x += lowbit(x);
}
}
long long query(int x) {
long long ans = 0;
while (x) {
ans += tree[x];
x -= lowbit(x);
}
return ans;
}
int main() {
freopen("in.txt", "r", stdin);
scanf("%d%d", &n, &m);
long long last = 0, now;
for (int i = 1; i <= n; i++) {
scanf("%lld", &now);
add(i, now - last);
last = now;
}
int flg;
while (m--) {
scanf("%d", &flg);
if (flg == 1) {
int x, y;
long long k;
scanf("%d%d%lld", &x, &y, &k);
add(x, k);
add(y + 1, -k);
} else if (flg == 2) {
int x;
scanf("%d", &x);
printf("%lld\n", query(x));
}
}
return 0;
}