#include <queue> #include <cstdlib> #include <cmath> #include <cstdio> #include <string> #include <cstring> #include <iostream> #include <algorithm> using namespace std; typedef long long ll; const int maxn = 1e6+9; int trie[maxn][26]; int cntword[maxn...

#include<bits/stdc++.h> using namespace std; typedef long long int ll; const int maxn=3e4+10; ll sun3(ll x,ll y) //异或 { if(abs(x-y)<4) { if(x>y)swap(x,y);ll sum=0; for(ll i=x;i<=y;i++)sum=sum^i; return sum; } return sun3(x-1-(x-1)%4,x-1)^sun3(y-y%4,y); } ll sun2(ll l,ll r){ ll sum=0; ...

#include<iostream> #include<cstdio> #include<cstring> #include<algorithm> #include<cmath> #define N 100010 using namespace std; const double pi=acos(-1); struct node{double x,y;}b[N]; node O; double R; int n; double sqr(double x){return x*x;} double dis(node x,node y){ ...

#include<bits/stdc++.h> #define PI acos(-1) #define ios ios::sync_with_stdio(false) #define endl '\n' #define srand srand((unsigned)time(NULL)) using namespace std; typedef long long int ll; typedef pair<int,int> PII; const int maxn=1e5+10; struct node{ int x,y,z; node(int fir,int sec,in...

#include <bits/stdc++.h> typedef long long ll; using namespace std; const int N=1e5+5; int n, cnt; int ls[N<<5], rs[N<<5], sum[N<<5], rt[N<<5]; int a[N], b[N]; int build(int l, int r){ int root = ++ cnt; if(l == r) return root; int mid = (l + r) >> 1; ls[root] = b...

#include<bits/stdc++.h> using namespace std; typedef long long int ll; const int maxn=400+10; const int char_long=400; typedef pair<int,int> pii; int num[maxn]; string add(string x,string y); string minuse(string x,string y); string int_string(int s){ //int->string if(s==0)return &quo...

#include <iostream> #include <algorithm> using namespace std; int main() { int num[3]={1,2,3}; do { cout<<num[0]<<" "<<num[1]<<" "<<num[2]<<endl; }while(next_permutation(num,num+3)); return 0; }

inline __int128 read() { __int128 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; } inline void write(__int128 x) { if (x < 0) { put...

#include<bits/stdc++.h> using namespace std; typedef long long int ll; const int maxn=1e5+10; int ch[maxn*31][2];//树（ch[i][0]表示） int val[maxn];//存数值 int a[maxn];//初始数据 int bx[maxn],tail=1; int tol=1; void insert(int x,int k){ //插入数据 int u=0; for(int i=31;i>=0;--i){ int v=(x>>i)&1;...

球体的一般方程：        #include<bits/stdc++.h> using namespace std; const double Pi=acos(-1); typedef double da; double solve(da x1,da y1,da z1,da x2,da y2,da z2,da r1,da r2){ double d=sqrt((x1-x2)*(x1-x2)+(y1-y2)*(y1-y2)+(z1-z2)*(z1-z2)); if(d>=r1+r2)return 0; if(r1+d<=r2||r2+d<=r1){ return...

int lar(int heights[],int n) { heights[++n]=-1; stack<int> st; int ret = 0, top; for (int i = 1; i <= n; i++) { if (st.empty() || heights[st.top()] <= heights[i]) { st.push(i); } else { while (!st.empty() && heights[st.top()] > heights[i]) { top = st.top(); st.pop(); int tmp =...

#include<bits/stdc++.h> #define fx for(i=0;i<30;i++)cout<<a[i]<<" ";cout<<endl; using namespace std; typedef complex<double> CP; const int lim = 1<<21; const int P = 500001; const double Pi = acos(-1); bool vis[lim]; //结果 void FFT(CP *x,int lim,int in...

#include<bits/stdc++.h> using namespace std; const int maxn=1e5+10; typedef long long int ll; int tree[maxn]; int n; void update(int x,int v){ //在x位置加上v while(x<=n){ tree[x]+=v; x+=(x&(-x)); } } int getsum(int x){ //求和前x位置 int ans = 0; while(x){ ans+=tree[x]; x-=(x&(-x)); } return a...

#include<bits/stdc++.h> using namespace std; const int maxn=1e4; int shu[maxn*4],bis[maxn*4];//bis标记 int qu[maxn]; int n; //1~~n //st=1,en=n,node=1; void built(int node,int st,int en){ int mod=(st+en)/2; if(st==en){ shu[node]=qu[st]; return ; } built(node*2,st,mod);built(node*2+1,mod+1,en); sh...