1. 程式人生 > >NOIP考試技巧及注意事項&&資訊學競賽常用函式/模板

NOIP考試技巧及注意事項&&資訊學競賽常用函式/模板


OI 中轉站
OI演算法大全


[NOIP提高組(2018)考試技巧及注意事項] (https://blog.csdn.net/hi_ker/article/details/81181615)

1.保持好的心態

2.審題

3.考慮嚴謹

4.程式碼正確

5.freopen

6.檢查

還有注意檢查你的程式是否在考試指定資料夾裡面,而且不能存在指定資料夾的子資料夾裡面

7.防爆(防止爆陣列)

8.標頭檔案

9.資料型別

10.增加程式的魯棒性(Robust)

11.注意隱性錯誤(非常重要)


以下討論主要是基於windows考試環境

1.時間複雜度的選擇

2.利用Dev C++

3.對拍

4.巧用cmd(選學)

5.查錯

6.騙分




[資訊學競賽常用函式/模板] (https://blog.csdn.net/Hi_KER/article/details/82973124)

  • 擴充套件歐幾里得演算法
void exgcd(int a,int b,int& d,int& x,int& y){
    if(!b){d=a,x=1,y=0;return;}
    exgcd(b,a%b,d,y,x);y-=x*(a/b);
}
  • 乘法逆元
void
exgcd(int a,int b,int& d,int& x,int& y){ if(!b){d=a,x=1,y=0;return;} exgcd(b,a%b,d,y,x);y-=x*(a/b); } int inverse(int a,int m){ int x,y,d; exgcd(a,m,d,x,y); return d==1?(x%m+m)%m:-1; }
  • 高精度模板

https://blog.csdn.net/Hi_KER/article/details/80982710

  • O
    N ) O(\sqrt {N})
    計算尤拉函式
int phi(int x){
	int ans=x,m=sqrt(x+0.01);
	for(int i=2;i<=m;i++) if(x%i==0){
		ans=ans/i*(i-1);
		while(x%i==0) x/=i;
	}
	if(x>1) ans=ans/x*(x-1);//x為質數
	return ans;
}

  • 線性時間複雜度生成尤拉函式表+質數表
/* 
p|n &&  p*p|n phi(n)=phi(n/p)*p
p|n && !p*p|n phi(n)=phi(n/p)*(p-1) 
*/ 
int v[maxn],pri[maxn],phi[maxn];
void eular(int n){
	int num=0;//number of prime 
	for(int i=2;i<=n;i++) {
		if(v[i]==0){//v[i]:i 的最小質因子
			v[i]=i,pri[++num]=i;
			phi[i]=i-1;
		}
		for(int j=1;j<=num;j++){
			if(pri[j]>v[i]||pri[j]>n/i) break;
			v[i*pri[j]]=pri[j];
			phi[i*pri[j]]=phi[i]*(i%pri[j]?pri[j]-1:pri[j]);
		}
	}
	for(int i=1;i<=n;i++) printf("phi(%d)=%d\n",i,phi[i]);
}

  • BSGS演算法(說明:用於計算高次同餘方程 a x b   ( m o d   p ) a ^ x \equiv b\ (mod\ p) 的最小非負數解,其中 p p 為素數,時間複雜度 O ( P ) O(\sqrt{P} ) ,返回 1 -1 表示無解)
int pow_mod(int a,int n,int p){
	int ret=1;
	while(n){
		if(n&1) ret=(long long)ret*a%p;
		a=(long long)a*a%p,n>>=1;
	}
	return ret;
}
map<int,int>mp;
int bsgs(int a,int b,int p){
	if(a%p==0) return b?-1:0;
	mp.clear();
	int m=ceil(sqrt(p)),T=b%p;
	mp[T]=0;
	for(int i=1;i<=m;i++)
		mp[T=(long long)T*a%p]=i;
	int t=pow_mod(a,m,p);T=1;
	for(int i=1;i<=m;i++){
		T=(long long)T*t%p;
		if(mp.count(T)) return i*m-mp[T];
	}
	return -1;
}

  • 盧卡斯(Lucas)定理
int C(int N,int M,int p){
    if(N<M) return 0;
    if(!n) return 1;
    return (ll)fac[N]*inv[M]*inv[N-M]%p;
}
int lucas(int N,int M,int p){
    if(M==0) return 1;
    return (ll)lucas(N/p,M/p,p)*C(N%p,M%p,p)%p;
}

  • 擴充套件中國剩餘定理
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int maxn=100005;
#define ll __int128
void exgcd(ll a,ll b,ll& d,ll& x,ll& y){
    if(!b){x=1,y=0,d=a;return;}
    exgcd(b,a%b,d,y,x),y-=x*(a/b);
}
ll gcd(ll a,ll b){
    return b?gcd(b,a%b):a;
}
ll lcm(ll a,ll b){
    return a*b/gcd(a,b);
}
ll a[maxn],b[maxn];
long long a__,b__;
int n;
int main(){
    scanf("%d",&n);
    for(int i=1;i<=n;i++) 
        scanf("%lld%lld",&a__,&b__),a[i]=a__,b[i]=b__;
    for(int i=1;i<n;i++){
        ll a_=a[i],b_=-a[i+1],c_=b[i+1]-b[i],d,x,y,md;
        exgcd(a_,b_,d,x,y);
        if(c_%d){break;/*No answer*/}
        x=c_/d*x,md=b_/d;
        if(md<0) md=-md;
        x=(x%md+md)%md;
        b[i+1]=a[i]*x+b[i];
        a[i+1]=lcm(a[i],a[i+1]);
    }
    printf("%lld\n",(long long)b[n]);
    return 0;
}

  • 矩陣運算
struct matrix{
	#define T int
	#define MOD 0x7fffffff
    int N,M;
	T mtx[maxn][maxn];
    matrix(){M=N=0;memset(mtx,0,sizeof(mtx));}
    void resize(int n,int m){N=n,M=m;}
    void unit_matrix(int n){
    	memset(mtx,0,sizeof(mtx));
    	N=M=n;
    	for(int i=0;i<n;i++) mtx[i][i]=1;
	}
    void input(){//just for int
        for(int i=0;i<N;i++)
          for(int j=0;j<M;j++)
            scanf("%d",&mtx[i][j]);
    }
    void output(){//just for int
        for(int i=0;i<M;i++){
        	for(int j=0;j<N;j++)
        		printf("%d ",mtx[i][j]);
        	putchar('\n');
        }
    }
    friend matrix operator * (matrix x,matrix y){
        matrix ans;
    	if(x.N!=y.M) return ans; //Error
        ans.N=x.N,ans.M=y.M;
        for(int i=0;i<x.N;i++)
          for(int j=0;j<y.M;j++){
            T S=0;
            for(int k=0;k<x.N;k++)
                S=(S+x.mtx[i][k]*y.mtx[k][j])%MOD;
            ans.mtx[i][j]=S%MOD;
        }
        return ans;
    }
    friend matrix operator ^ (matrix x,int n){
    	matrix ans;
    	if(x.N!=x.M) return ans;//Error
    	ans.unit_matrix(x.N);
    	while(n){
    		if(n&1) ans=ans*x;
    		x=x*x,n>>=1;
		}
		return ans;
	}
    #undef T
    #undef MOD
};


資料結構模組

  • 並查集
struct union_set{
	int f[maxn];
	void init(int n){for(int i=1;i<=n;i++) f[i]=i;}
	int find(int x){return x==f[x]?x:f[x]=find(f[x]);}
	bool merge(int x,int y){
		int fx=find(x),fy=find(y);
		if(fx!=fy) f[fx]=fy;
		return fx==fy;
	}
};

  • 連結串列 (陣列版)
struct node{
    int v,pre,nxt;
};
struct link_list{
    node E[maxn];
    int now,head,tail;
    link_list(){
        now=2,head=1,tail=2;
        E[head].nxt=tail,E[tail].pre=head;
    }
    int insert(int pos,int v){
        now++,E[now].v=v;
        E[E[pos].nxt].pre=now,E[now].nxt=E[pos].nxt;
        E[pos].nxt=now,E[now].pre=pos;
        return now;
    }
    int insert_to_tail(int v){
        return insert(E[tail].pre,v);
    }
    void del(int pos){
        E[E[pos].pre].nxt=E[pos].nxt;
        E[E[pos].nxt].pre=E[pos]