1//f[]：可以取走的石子个数 //sg[]:0~n的SG函数值 //mex[]:mex{} int f[N],sg[N],mex[N]; void getSG(int n){ int i,j;memset(sg,0,sizeof(sg)); for(i=1;i<=n;i++){ memset(mex,0,sizeof(mex)); for(j=1;f[j]<=i;j++)mex[sg[i-f[j]]]=1; for(j=0;j<=n;j++)//求mes{}中未出现的最小的非负整数 if(mex[j]==0){sg[i]=j; break;} } } 2//注意 S数组要...

RMQ#include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define scnaf scanf #define itn int using namespace std; const int o_o=5e4+5; const int mod=1e9+7; const int oo=0x7fffffff; const int sup=0x80000000; typedef long long ll; typedef unsigned long long ull; int dpmax[o_o][21]; int d...

序列自动机for(int i=n;i;i--){ for(int j=0;j<26;j++) next[i-1][j]=next[i][j]; next[i-1][s[i]-'a']=i; } 最大最小表示法int min_max_express(char *s,int len,bool flag) { int i=0,j=1,k=0; while(i<len&&j<len&&k<len) { int t=s[(j+k)%len]-s[(i+k)%len]; if(t==0) k++; else { if(flag) { if(t>...

query       #include<bits/stdc++.h> using namespace std; const int N=2e6+5; int n,m=0,q,ans[N],pos[N]; struct node{ int x,y,id; bool friend operator<(node l,node r){ return l.y<r.y; } }qx[N],qy[N]; int sum[N]={0}; void up(int x){for(int i=x;i<=n;i+=i&-i)sum[i]++;} int Q(int x){ in...

HDU-4059 The Boss on Mars       #include<bits/stdc++.h> using namespace std; const int mod=1e9+7; const int inv=233333335; typedef long long ll; ll f(ll n){ ll ans=n*(n+1)%mod*(2LL*n+1)%mod; ll x=(3LL*n*n%mod+3LL*n-1)%mod*inv%mod; return ans*x%mod; } int mu(ll n){ if(n==1)return 1; int cnt=0; ...

The beautiful values of the palace  二维子矩阵和(坐标太大).      #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define sc scanf using namespace std; const int N=1e6+5; typedef long long ll; struct matrix{ ll TD(ll x,ll y,ll n){ ll r=0; if(x<=y&&x+y<=n+1){ r=x;return 4LL*(r-...

Easy Math        #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define sc scanf #define itn int #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << double...

Sum        #include<bits/stdc++.h> using namespace std; const int N=2e7+7; typedef long long ll; int prime[N],tot=0; ll f[N]={0}; bool vis[N]={0}; void pre(){ f[1]=1; for(int i=2;i<N;i++){ if(!vis[i])prime[++tot]=i,f[i]=2; for(int j=1;j<=tot&&i*prime[j]<N;j++){ vis[i*prime[j]]...

Convex Hull        #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define sc scanf #define itn int #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << doub...

K Sum        #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define sc scanf #define itn int #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << double(end...

G. tsy's number        #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define sc scanf #define itn int #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << ...

HDU-6390 GuGuFishtion    求逆元时线性递推即可。   #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define sc scanf #define itn int #define STR clock_t startTime = clock(); #define END clock_t endTime = clock(...

BZOJ-3529 数表       #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define sc scanf #define itn int #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << doub...

BZOJ-2440 完全平方数       #include<bits/stdc++.h> #define me(a,x) memset(a,x,sizeof(a)) #define sc scanf #define itn int #define IN freopen("in.txt","r",stdin); #define OUT freopen("out.txt","w",stdout); #define STR clock_t startTime = clock(); #define END clock_t endTime = clock();cout << d...