NewOJ Week 3---BCDE

B.Radical of N

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
using namespace std;
#define int long long
const int N = 1000010 , mod = 1e18;

int p[N],rad[N],s[N];
void init()
{
	for(int i=1;i<N;i++)rad[i]=1;
	for(int i=2;i<N;i++)
	{
		if(!p[i])
		{
			rad[i]=i;
			for(int j=i+i;j<N;j+=i)
				rad[j]*=i,p[j]=1;
		}
	}
	for(int i=1;i<N;i++)s[i]=s[i-1]+rad[i];
}
signed main()
{
	init();
	int T;cin>>T;
	while(T--)
	{
		int l,r;cin>>l>>r;
		cout<<s[r]-s[l-1]<<'\n';
	}
}

C.二进制中的1—数位DP+扩展欧拉定理

#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define int long long
const int N = 55 , mod = 10000007;
int qmi(int a,int b){int res=1;a=a%mod;while(b){if(b&1)res=res*a%mod;a=a*a%mod;b>>=1;}return res%mod;}
int get_oula(int x){int res=x;for(int i=2;i<=x/i;i++){if(x%i==0){res=res/i*(i-1);while(x%i==0)x/=i;}}if(x>1)res=res/x*(x-1);return res;}

//组合数
int c[N][N];
void init(){
    for(int i=0;i<N;i++){
        for(int j=0;j<=i;j++){
            if(!j)c[i][j]=1;
            else c[i][j]=(c[i-1][j-1]+c[i-1][j])%get_oula(mod);
        }
    }
}
//数位dp 
int dp(int n)
{
	vector<int>nums;
	int res=1,pre_sum=0;
	while(n)nums.push_back(n&1),n>>=1;
	for(int i=nums.size()-1;i>=0;i--)
	{
		int x=nums[i];
		if(x)
		{
			for(int j=0;j<=i;j++)
				if(j+pre_sum)res=res*qmi(j+pre_sum,c[i][j])%mod;
			pre_sum++;
		}
	}
	res=res*pre_sum%mod;
	return res;
}

signed main()
{
	init();
	int n;cin>>n;
	cout<<dp(n)<<'\n';
}

D. 最小生成树—Kruskal+贪心

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
using namespace std;
#define int long long
const int N = 100010, mod = 10000007;

int p[N];
int find(int x)
{
	if(p[x]!=x)p[x]=find(p[x]);
	return p[x];
}
struct Ed{
	int a,b,w;
	bool operator < (const Ed &W)const
	{
		return w<W.w;
	}
};
pair<int,int> vec[N];
vector<Ed>e;
int n;
int last[15];

void add_edge(int u,int v)
{
	int dis=(vec[u].first-vec[v].first)*(vec[u].first-vec[v].first)+(vec[u].second-vec[v].second)*(vec[u].second-vec[v].second);
	e.push_back({u,v,dis});
}
void kruskal()
{
	for(int i=1;i<=n;i++)p[i]=i;
	sort(e.begin(),e.end());
	int res=0;
	for(int i=0;i<e.size();i++)
	{
		int a=find(e[i].a),b=find(e[i].b);
		if(a!=b)
		{
			p[a]=b;
			res+=e[i].w;
		}
	}
	cout<<res<<'\n';
}
signed main()
{
	cin>>n;
	for(int i=1,a,b;i<=n;i++)
	{
		cin>>a>>b;
		vec[i]={a,b};
	}
	sort(vec+1,vec+n+1);
	memset(last,-1,sizeof last);
	for(int i=1;i<=n;i++)
	{
		for(int y=0;y<=10;y++)
			if(last[y]!=-1)
				add_edge(i,last[y]);
		last[vec[i].second]=i;
	}
	memset(last,-1,sizeof last);
	for(int i=n;i>=1;i--)
	{
		for(int y=0;y<=10;y++)
			if(last[y]!=-1)
				add_edge(i,last[y]);
		last[vec[i].second]=i;
	}
	kruskal();
}

E.黑白边—dfs树上差分+树状数组

#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
#define int long long
const int N = 250010 , mod = 10000007;
int n,m,cnt;
int l[N],r[N];
vector<int>g[N];
int tr[2*N];
void dfs(int u)
{
    //cout<<u<<" ";
	l[u]=++cnt;
	for(auto v : g[u])
		dfs(v);
	r[u]=++cnt;
	//cout<<u<<" ";
}
int lowbit(int x)
{
	return x&-x;
}
void add(int x,int v)
{
	for(int i=x;i<=cnt;i+=lowbit(i))tr[i]+=v;
}
int query(int x)
{
	int res=0;
	for(int i=x;i>0;i-=lowbit(i))res+=tr[i];
	return res;
}
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0);cout.tie(0);
	cin>>n;
	for(int i=1,u,v;i<n;i++)
	{
		cin>>u>>v;
		g[u].push_back(v);
	}
	
	dfs(1);
	for(int i=1;i<=n;i++)
		add(l[i],1),add(r[i],-1);
	
	cin>>m;
	m=m+n-1;
	while(m--)
	{
		char op[2];
		cin>>op;
		if(op[0]=='A')
		{
			int x,y;cin>>x>>y;
			if(x>y)swap(x,y);
			add(l[y],-1);
			add(r[y],1);
		}
		else
		{
			int x;cin>>x;
			cout<<query(l[x])-1<<'\n';
		}
	}
}