代码大全

代码大全

常数优化

编译器卡常优化（CF上比较有用）

#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#pragma GCC target("sse,sse2,sse3,ssse3,sse4,popcnt,abm,mmx,avx,tune=native")

读入优化：

inline int getint() {
    register char ch;
    while(!isdigit(ch=getchar()));
    register int x=ch^'0';
    while(isdigit(ch=getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}

mmap读入优化：

#include<sys/mman.h>
#include<sys/stat.h>
class MMapInput {
    private:
        char *buf,*p;
        int size;
    public:
        MMapInput() {
            register int fd=fileno(stdin);
            struct stat sb;
            fstat(fd,&sb);
            size=sb.st_size;
            buf=reinterpret_cast<char*>(mmap(0,size,PROT_READ,MAP_PRIVATE,fileno(stdin),0));
            p=buf;
        }
        char getchar() {
            return (p==buf+size||*p==EOF)?EOF:*p++;
        }
};
MMapInput mmi;
inline int getint() {
    register char ch;
    while(!isdigit(ch=mmi.getchar()));
    register int x=ch^'0';
    while(isdigit(ch=mmi.getchar())) x=(((x<<2)+x)<<1)+(ch^'0');
    return x;
}

数论

扩欧求逆元：

void exgcd(const int &a,const int &b,int &x,int &y) {
    if(!b) {
        x=1,y=0;
        return;
    }
    exgcd(b,a%b,y,x);
    y-=a/b*x;
}
inline int inv(const int &x) {
    int ret,tmp;
    exgcd(x,mod,ret,tmp);
    return (ret%mod+mod)%mod;
}

多项式

FFT

inline void init_ru(const int &n) {
    for(register int i=0;i<n;i++) {
        ru[i]=(comp){cos(2*pi*i/n),sin(2*pi*i/n)};
        iru[i]=conj(ru[i]);
    }
}
inline void dft(comp f[],comp w[],const int &n) {
    for(register int i=0,j=0;i<n;i++) {
        if(i>j) std::swap(f[i],f[j]);
        for(register int l=n>>1;(j^=l)<l;l>>=1);
    }
    for(register int i=2;i<=n;i<<=1) {
        const int m=i>>1;
        for(register int j=0;j<n;j+=i) {
            for(register int k=0;k<m;k++) {
                const comp z=f[j+m+k]*w[n/i*k];
                f[j+m+k]=f[j+k]-z;
                f[j+k]+=z;
            }
        }
    }
}

FWT异或：

inline void fwt(int a[]) {
    for(register int j=0;j<len;j++) {
        for(register int i=0;i<n;i++) {
            if(i>>j&1) continue;
            const int x=a[i],y=a[(1<<j)|i];
            a[i]=(x+y)%mod;
            a[(1<<j)|i]=(x-y+mod)%mod;
        }
    }
}

字符串

后缀数组求sa[i]和rank[i]：

inline bool cmp(const int &i,const int &j) {
    if(rank[i]!=rank[j]) return rank[i]<rank[j];
    const int ri=i+k<=n?rank[i+k]:-1;
    const int rj=j+k<=n?rank[j+k]:-1;
    return ri<rj;
}
inline void suffix_sort() {
    for(register int i=0;i<=n;i++) {
        sa[i]=i;
        rank[i]=s[i];
    }
    for(k=1;k<=n;k<<=1) {
        std::sort(&sa[0],&sa[n]+1,cmp);
        tmp[sa[0]]=0;
        for(register int i=1;i<=n;i++) {
            tmp[sa[i]]=tmp[sa[i-1]]+!!cmp(sa[i-1],sa[i]);
        }
        std::copy(&tmp[0],&tmp[n]+1,rank);
    }
}

后缀数组求lcp[i]：

inline void init_lcp() {
    for(register int i=0,h=0;i<n;i++) {
        if(h>0) h--;
        const int j=sa[rank[i]-1];
        while(j+h<n&&i+h<n&&s[j+h]==s[i+h]) h++;
        lcp[rank[i]-1]=h;
    }
}

计算几何：

自适应辛普森法：

inline double simpson(const double &a,const double &b) {
    const double c=(a+b)/2;
    return (F(a)+4*F(c)+F(b))*(b-a)/6;
}
inline double asr(const double &a,const double &b,const double &eps,const double &s) {
    const double c=(a+b)/2,ls=simpson(a,c),rs=simpson(c,b);
    if(fabs(ls+rs-s)<15*eps) return ls+rs+(ls+rs-s)/15;
    return asr(a,c,eps/2,ls)+asr(c,b,eps/2,rs);
}

转载于:https://www.cnblogs.com/skylee03/p/9154114.html