【联合省选 2020 A】题解

最新推荐文章于 2022-04-22 21:04:00 发布
最新推荐文章于 2022-04-22 21:04:00 发布
阅读量743 收藏
点赞数 1
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42555009/article/details/106947663
版权

冰火战士

线段树上二分即可
具体实现可以参见代码

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}
char obuf[(int)(5e7+7)],*oh=obuf,ch[23];
template<typename T>void write(T a,char c){
	if(a<0)*oh++='-',a=-a;int tl=0;
	do ch[++tl]=a%10;while(a/=10);
	while(tl)*oh++=ch[tl--]^48;*oh++=c; 
}
void writes(char T[],char c){
	for(int i=0;T[i]!='\0';i++)*oh++=T[i];
	*oh++=c;
}
struct obuf_flusher{~obuf_flusher(){fwrite(obuf,1,oh-obuf,stdout);}}Flusher;

}
using IO::read;
using IO::write;
using IO::writes;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}

cs int N=2e6+5;
int n,m,a[N],b[N],fr[N],pf[N];
namespace seg{

cs int N=::N<<2;
int sa[N],sb[N],fir[N];
#define lc (u<<1)
#define rc ((u<<1)|1)
int p,kd,k;
void upd(int u,int l,int r){
	if(l==r){
		if(!kd)sa[u]+=k;
		else sb[u]+=k;
		fir[u]=sb[u]?l:-1;
		return;
	}int mid=(l+r)>>1;
	if(p<=mid)upd(lc,l,mid);
	else upd(rc,mid+1,r);
	sa[u]=sa[lc]+sa[rc],sb[u]=sb[lc]+sb[rc];
	fir[u]=sb[lc]?fir[lc]:fir[rc];
}
void update(int _p,int _kd,int _k){
	p=_p,kd=_kd,k=_k,upd(1,1,n);
}
pii query(int u,int l,int r,int _sa,int _sb){
	if(l==r)return pii(l,2*min(_sa+sa[u],_sb+sb[u]));
	int mid=(l+r)>>1,ls=_sa+sa[lc],rs=_sb+sb[rc];
	if(_sa+sa[lc]<=_sb+sb[rc])return query(rc,mid+1,r,ls,_sb);
	else {
		pii x=query(lc,l,mid,_sa,rs);
		if(x.se==2*min(ls,rs))x.fi=fir[rc];
		return x;
	}
}
#undef lc
#undef mid
}
struct opt{int kd,x,y;}p[N];
int lower(int x){
	int l=1,r=n,res=n;
	while(l<=r){
		int mid=(l+r)>>1;
		if(b[mid]>=x)res=mid,r=mid-1;
		else l=mid+1;
	}return res;
}
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	freopen("my.out","w",stdout);
	#endif
	m=read();
	for(int i=1;i<=m;i++){
		int op=read();
		if(op==1){
			p[i].kd=read(),p[i].x=read(),p[i].y=read();
			b[++n]=p[i].x;
		}
		else fr[i]=read(),p[i]=p[fr[i]],p[i].y=-p[i].y;
	}sort(b+1,b+n+1);
	n=unique(b+1,b+n+1)-b-1;
	for(int i=1;i<=m;i++){
		int ps;
		if(p[i].y>=0)ps=lower(p[i].x);
		else ps=pf[fr[i]];pf[i]=ps;
		seg::update(ps,p[i].kd,p[i].y);
		pii x=seg::query(1,1,n,0,0);
		if(!x.se)writes("Peace",'\n');
		else write(b[x.fi],' '),write(x.se,'\n');
	}return 0;
}

组合数问题

直接斯特林展开然后二项式定理合并组合数即可
如果转成下降幂多项式做可以两个log，如果利用组合意义/推导可以一个log

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}
}

using IO::read;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}

int mod;
inline int add(int a,int b){return a+b>=mod?a+b-mod:a+b;}
inline int dec(int a,int b){return a<b?a-b+mod:a-b;}
inline int mul(int a,int b){return (ll)a*b%mod;}
inline void Add(int &a,int b){a=a+b>=mod?a+b-mod:a+b;}
inline void Dec(int &a,int b){a=a<b?a-b+mod:a-b;}
inline void Mul(int &a,int b){a=(ll)a*b%mod;}
inline int ksm(int a,int b,int res=1){for(;b;b>>=1,Mul(a,a))if(b&1)Mul(res,a);return res;}
inline int Inv(int x){return ksm(x,mod-2);}
cs int N=1005;
int n,x,m,a[N],s[N][N];
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	#endif
	n=read(),x=read(),mod=read(),m=read();
	s[0][0]=1;
	for(int i=0;i<=m;i++)a[i]=read();
	for(int i=1;i<=m;i++)
	for(int j=1;j<=i;j++)s[i][j]=add(s[i-1][j-1],mul(s[i-1][j],j));
	int res=0;
	for(int t=0,mtx=1,dn=1;t<=m;Mul(dn,n-t),Mul(mtx,x),t++){
		int nw=0;
		for(int j=t;j<=m;j++)Add(nw,mul(a[j],s[j][t]));
		Add(res,mul(nw,mul(dn,mul(mtx,ksm(add(1,x),n-t)))));
	}cout<<res;return 0;
}

魔法商店

首先可以利用拟阵推出偏序关系
其实就是不同元素在满秩线性基中的相互替换
然后就变成 L 2 L_2 L2​保序回归问题了
由于是整数，直接整数整体二分然后最小割 c h e c k check check即可

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
typedef unsigned long long ul;
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}
inline ul readul(){
	char ch=gc();
	ul res=0;
	while(!isdigit(ch))ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return res;
}

}

using IO::read;
using IO::readul;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}

cs int N=1055,M=66,inf=1e9;
int S,T,tot;
bool vis[N];
namespace flow{

struct edge{
	int v,cap,r;
	edge(int _v=0,int _c=0,int _r=0):v(_v),cap(_c),r(_r){}
};
vector<edge>e[N];
typedef vector<edge>::iterator It;
It tp[N];int lev[N];
inline void addedge(int u,int v,int cap){
	e[u].pb(edge(v,cap,e[v].size()));
	e[v].pb(edge(u,0,e[u].size()-1));
}
void clearflow(){
	for(int i=1;i<=tot;i++)e[i].clear();
	memset(vis,0,sizeof(bool)*(tot+1));
	tot=0;
}
bool bfs(){
	memset(lev,-1,sizeof(int)*(tot+1));
	queue<int> q;lev[S]=0,q.push(S);
	while(q.size()){
		int u=q.front();q.pop();
		for(cs edge &x:e[u])if(lev[x.v]==-1&&x.cap>0){
			lev[x.v]=lev[u]+1,q.push(x.v);
		}
	}return lev[T]!=-1;
}
int dfs(int u,int flow){
	if(u==T||!flow)return flow;
	int res=0;
	for(It &it=tp[u];it!=e[u].end();++it){
		if(lev[it->v]==lev[u]+1&&it->cap>0){
			int now=dfs(it->v,min(it->cap,flow-res));
			res+=now,it->cap-=now,e[it->v][it->r].cap+=now;
			if(flow==res)break;
		}
	}return res;
}
inline void dinic(){
	while(bfs()){
		for(int i=1;i<=tot;i++)tp[i]=e[i].bg();
		dfs(S,inf);
	}
}
void color(int u){
	vis[u]=1;
	for(cs edge &x:e[u])if(x.cap>0&&!vis[x.v])color(x.v);
}

}
using flow::addedge;
using flow::dinic;
using flow::color;
using flow::clearflow;

namespace Bas{

ul b[M];
void clearbas(){memset(b,0,sizeof(b));}
void insert(ul x){
	for(int i=63;~i;i--)if((x>>i)&1){
		if(b[i])x^=b[i];
		else{b[i]=x;break;}
	}
}
bool check(ul x){
	for(int i=63;~i;i--)if((x>>i)&1)x^=b[i];
	return (!!x);
}

}
using Bas::clearbas;
using Bas::insert;
using Bas::check;
int n,m,p[N],v[N],a[M],b[M],inst[N];
ul c[N];
vector<int> lk[N];
ll ans;
ll P(int x){return (ll)x*x;}
void solve(int l,int r,int ql,int qr){
	if(l>r)return;
	if(ql==qr){
		for(int i=l;i<=r;i++)ans+=P(v[p[i]]-ql);
		return;
	}
	int mid=(ql+qr)>>1;
	clearflow(),tot=r-l+1,S=++tot,T=++tot;
	for(int i=l;i<=r;i++)inst[p[i]]=i-l+1;
	for(int t=l,i;t<=r;t++){
		i=p[t];int vl=2*(v[i]-mid)-1;
		if(vl>0)addedge(inst[i],T,vl);
		else addedge(S,inst[i],-vl);
//		if(v[i]>mid)addedge(inst[i],T,P(v[i]-mid)-P(v[i]-mid-1));
//		else addedge(S,inst[i],P(mid+1-v[i])-P(mid-v[i]));
		for(int v:lk[i])if(inst[v])addedge(inst[i],inst[v],inf);
	}
	for(int i=l;i<=r;i++)inst[p[i]]=0;
	dinic(),color(S);int ps=l;
	for(int i=l;i<=r;i++)
	if(vis[i-l+1])swap(p[i],p[ps++]);
	solve(l,ps-1,ql,mid),solve(ps,r,mid+1,qr);
}
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	#endif
	n=read(),m=read();
	for(int i=1;i<=n;i++)c[i]=readul();
	for(int i=1;i<=n;i++)v[i]=read();
	for(int i=1;i<=m;i++)a[i]=read();
	for(int i=1;i<=m;i++)b[i]=read();
	for(int i=1;i<=m;i++)inst[a[i]]=1;
	for(int i=1;i<=m;i++){
		clearbas();
		for(int j=1;j<=m;j++)if(i^j)insert(c[a[j]]);
		for(int j=1;j<=n;j++)if(!inst[j]&&check(c[j]))lk[j].pb(a[i]);
		//,cout<<j<<" "<<a[i]<<'\n';
	}
//	cout<<'\n';
	for(int i=1;i<=m;i++)inst[a[i]]=0;
	for(int i=1;i<=m;i++)inst[b[i]]=1;
	for(int i=1;i<=m;i++){
		clearbas();
		for(int j=1;j<=m;j++)if(i^j)insert(c[b[j]]);
		for(int j=1;j<=n;j++)if(!inst[j]&&check(c[j]))lk[b[i]].pb(j);
		//,cout<<b[i]<<" "<<j<<'\n';
	}
	for(int i=1;i<=m;i++)inst[b[i]]=0;
	for(int i=1;i<=n;i++)p[i]=i;
	solve(1,n,0,1e6);
	cout<<ans<<'\n';
	return 0;
}

信号传递

状压 d p dp dp，预处理贡献卡空间，分块一下即可

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}

}

using IO::read;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}

cs int N=24,M=(1<<23)|5,B=(1<<12)|5;
int lk[N][N],coef11[N][B],coef12[N][B],coef21[N][B],coef22[N][B];
int ls,n,vk,str[100005],lim,lm1,lm2,ln1,ln2;
int f[M],ct[M],cnt[M];
inline int coef1(int i,int s){
	return coef11[i][s&(lm1-1)]+coef12[i][(s>>12)&(lm2-1)];
}
inline int coef2(int i,int s){
	return coef21[i][s&(lm1-1)]+coef22[i][(s>>12)&(lm2-1)];
}
#define lb(x) (x&(-x))
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	#endif
	ls=read(),n=read(),vk=read(),lim=1<<n;
	ln1=min(n,12),ln2=n-ln1;
	lm1=1<<ln1,lm2=1<<ln2;
	for(int i=1;i<=n;i++)ct[1<<(i-1)]=i;
	for(int i=0;i<lim;i++)cnt[i]=cnt[i>>1]+(i&1);
	for(int i=1;i<=ls;i++)
		str[i]=read();
	for(int i=1;i<ls;i++)
		lk[str[i]][str[i+1]]++;
	for(int i=1;i<=n;i++){
		for(int s=0;s<lm1;s++)
		for(int j=1;j<=ln1;j++)if(s&(1<<(j-1))){
			coef11[i][s]+=lk[i][j]*vk+lk[j][i];
			coef21[i][s]+=lk[j][i]*vk-lk[i][j];
		}
		for(int s=0;s<lm2;s++)
		for(int j=1;j<=ln2;j++)if(s&(1<<(j-1))){
			coef12[i][s]+=lk[i][j+ln1]*vk+lk[j+ln1][i];
			coef22[i][s]+=lk[j+ln1][i]*vk-lk[i][j+ln1];
		}
	}memset(f,127/2,sizeof(f));
	f[0]=0;
	for(int s=1;s<lim;s++){
		int &res=f[s];
		for(int t=s,i=ct[lb(t)];t;t-=lb(t),i=ct[lb(t)]){
			chemn(res,f[s-(1<<(i-1))]+(coef1(i,s-(1<<(i-1)))+coef2(i,(lim-1)^s))*cnt[s]);
		}
	}cout<<f[lim-1]<<'\n';
	return 0;
}

树

直接 0 / 1 t r i e 0/1trie 0/1trie合并维护子树答案，+1用之前agc的套路低位到高位维护，直接 s h i f t shift shift即可

另一种做法是
考虑每个点每一位对祖先有贡献的是一段一段的，在倍增数组上差分贡献最后统计答案也可以做到 O ( n l o g n ) O(nlogn) O(nlogn)

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}

}

using IO::read;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}
char xx;
cs int N=525015;
namespace trie{

cs int N=::N*20;
int lc[N],rc[N],xo[N],sm[N],tot;
bool sz[N];
void pushnow(int u,int dep,int k){
	xo[u]^=k;if(k&(1<<dep))swap(lc[u],rc[u]);
}
void pushdown(int u,int dep){
	if(!xo[u])return;
	pushnow(lc[u],dep+1,xo[u]);
	pushnow(rc[u],dep+1,xo[u]);
	xo[u]=0;
}
void merge(int &u,int r1,int r2,int dep){
	if(!r1||!r2){u=r1+r2;return;}
	pushdown(r1,dep),pushdown(r2,dep);
	u=r1,sz[u]=sz[r1]^sz[r2];
	merge(lc[u],lc[r1],lc[r2],dep+1);
	merge(rc[u],rc[r1],rc[r2],dep+1);
	sm[u]=(sm[lc[u]]^sm[rc[u]])+(1<<dep)*sz[rc[u]];
}
void insert(int &u,int k,int dep){
	if(!u)u=++tot;sz[u]^=1;if(dep==21)return;
	if(!((k>>dep)&1))insert(lc[u],k,dep+1);
	else insert(rc[u],k,dep+1);
	sm[u]=(sm[lc[u]]^sm[rc[u]])+(1<<dep)*sz[rc[u]];
}
void shift(int u,int dep){
	if(!u)return;
	pushdown(u,dep);
	swap(lc[u],rc[u]);
	shift(lc[u],dep+1);
	sm[u]=(sm[lc[u]]^sm[rc[u]])+(1<<dep)*sz[rc[u]];
}

}
using trie::insert;
using trie::merge;
using trie::pushnow;
using trie::shift;
using trie::sm;
char yy;
int rt[N],n,val[N];
ll ans;
vector<int> e[N];
void dfs(int u){
	for(int v:e[u]){
		dfs(v);
		merge(rt[u],rt[u],rt[v],0);
	}
	shift(rt[u],0);
	insert(rt[u],val[u],0);
	ans+=sm[rt[u]];
//	cout<<u<<" "<<ans<<'\n';
//	if(u==2)exit(0);
}
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	#endif
	n=read();
	for(int i=1;i<=n;i++)val[i]=read();
	for(int i=2;i<=n;i++)e[read()].pb(i);
	dfs(1);
	cout<<ans<<'\n';
	return 0;
}

作业题

被玩烂的套路
枚举每个权值把倍数的边提出来矩阵树，最后容斥回答案
矩阵树如果对占位多项式插值是 O ( n 4 ) O(n^4) O(n4)的
如果用二项式定理做，即每个是 1 + w x 1+wx 1+wx在 % x 2 \%x^2 %x2下直接求多项式即可单次 O ( n 3 ) O(n^3) O(n3)

#include<bits/stdc++.h>
using namespace std;
#define cs const
#define pb push_back
#define pii pair<int,int>
#define fi first
#define se second
#define ll long long
#define bg begin
namespace IO{

cs int rlen=1<<22|1;
char ibuf[rlen],*ib,*ob;
inline char gc(){
	(ib==ob)&&(ob=(ib=ibuf)+fread(ibuf,1,rlen,stdin));
	return (ib==ob)?EOF:*ib++;
}
inline int read(){
	char ch=gc();
	int res=0;bool f=0;
	while(!isdigit(ch))f=ch=='-',ch=gc();
	while(isdigit(ch))res=(res*10)+(ch^48),ch=gc();
	return f?-res:res;
}

}

using IO::read;
template<typename tp>inline void chemx(tp &a,tp b){(a<b)?(a=b):0;}
template<typename tp>inline void chemn(tp &a,tp b){(a>b)?(a=b):0;}

cs int N=34,M=152511;
cs int mod=998244353;
inline int add(int a,int b){return a+b-(mod&-(a+b>=mod));}
inline int dec(int a,int b){return a-b+(mod&-(a<b));}
inline int mul(int a,int b){static ll r;r=(ll)a*b;return r>=mod?r%mod:r;}
inline void Add(int &a,int b){a=a+b>=mod?a+b-mod:a+b;}
inline void Dec(int &a,int b){a=a<b?a-b+mod:a-b;}
inline void Mul(int &a,int b){a=(ll)a*b%mod;}
inline int ksm(int a,int b,int res=1){for(;b;b>>=1,Mul(a,a))if(b&1)Mul(res,a);return res;}
inline void ex_gcd(int a,int b,int &x,int &y){
	if(!b){x=1,y=0;return;}ex_gcd(b,a%b,y,x);y-=a/b*x;
}
inline int Inv(int a){
	int x,y;ex_gcd(mod,a,y,x);
	return x+(x>>31&mod);
}
int gcd(int a,int b){
	return b?gcd(b,a%b):a;
}
bool vis[M];
int pri[M],pcnt,lim;
inline void init_sieve(cs int n){
	for(int i=2;i<=n;i++){
		if(!vis[i])pri[++pcnt]=i;
		for(int j=1;j<=pcnt&&i*pri[j]<=n;j++){
			vis[i*pri[j]]=1;
			if(i%pri[j]==0)break;
		}
	}
}
inline void fgt(int *f){
	for(int i=1;i<=pcnt;i++)
	for(int j=1,l=lim/pri[i];j<=l;j++)Dec(f[j],f[j*pri[i]]);
}
struct node{
	int a,b;
	node(int _a=0,int _b=0):a(_a),b(_b){}
	bool is_z()cs{return !a&&!b;}
	friend inline node operator +(cs node &a,cs node &b){
		return node(add(a.a,b.a),add(a.b,b.b));
	}
	friend inline node operator -(cs node &a,cs node &b){
		return node(dec(a.a,b.a),dec(a.b,b.b));
	}
	friend inline node operator *(cs node &a,cs int &b){
		return node(mul(a.a,b),mul(a.b,b));
	}
	friend inline node operator *(cs node &a,cs node &b){
		return node(mul(a.a,b.a),add(mul(a.b,b.a),mul(a.a,b.b)));
	}
	void operator +=(cs node &b){*this=*this+b;}
	void operator -=(cs node &b){*this=*this-b;}
};
node Inv(node x){
	if(!x.b)return node(Inv(x.a),0);
	int c=Inv(x.a),d=dec(0,mul(x.b,mul(c,c)));
	return node(c,d);
}
node a[N][N];
int u[N*N],v[N*N],w[N*N];
int n,m,f[M];
int pr[N],cnt[N],tot;
vector<int> eid[M];
void dfs(int ps,int mt,int id){
	if(ps>tot)return eid[mt].pb(id);
	for(int i=0;i<=cnt[ps];i++){
		dfs(ps+1,mt,id),mt*=pr[ps];
	}
}
void divid(int k,int id){tot=0;int x=k;
	for(int i=2;i*i<=x;i++)if(x%i==0){
		pr[++tot]=i;cnt[tot]=0;
		while(x%i==0)cnt[tot]++,x/=i;
	}if(x>1)pr[++tot]=x,cnt[tot]=1;
	dfs(1,1,id);
}
int fa[N];
int find(int x){return fa[x]==x?x:fa[x]=find(fa[x]);}
void merge(int u,int v){
	u=find(u),v=find(v);
	if(u==v)return;
	fa[u]=v;
}
inline int Det(int n){
	node res(1,0);
	for(int i=1;i<=n;i++){
		int ps=i;
		for(;ps<=n;ps++)if(!a[ps][i].is_z()){ps=i;break;}
		if(ps!=i)swap(a[ps],a[i]),res=node()-res;
		node iv=Inv(a[i][i]);res=res*a[i][i];
		for(int j=i+1;j<=n;j++){
			node mt=a[j][i]*iv;
			for(int k=i;k<=n;k++)
			a[j][k]-=a[i][k]*mt;
		}
	}return res.b;
}
inline void solve(int k){
	int vl=0;
	for(int i=1;i<=n;i++)fa[i]=i;
	for(int i:eid[k])merge(u[i],v[i]),vl=gcd(vl,w[i]);
	for(int i=1;i<=n;i++)if(find(i)!=find(1))return;
	if(vl!=k)return;
	for(int i=1;i<=n;i++)
	for(int j=1;j<=n;j++)a[i][j]=node(0,0);
	for(int i:eid[k]){
		node vl(1,w[i]);
		a[u[i]][u[i]]+=vl,a[v[i]][v[i]]+=vl;
		a[u[i]][v[i]]-=vl,a[v[i]][u[i]]-=vl;
	}
	f[k]=Det(n-1);
	for(int i=k+k;i<=lim;i+=k)Dec(f[k],f[i]);
}
int main(){
	#ifdef Stargazer
	freopen("lx.in","r",stdin);
	#endif
	n=read(),m=read();
	for(int i=1;i<=m;i++){
		u[i]=read(),v[i]=read(),w[i]=read();
		divid(w[i],i),chemx(lim,w[i]);
	}init_sieve(lim);
	for(int i=lim;i;i--)if(eid[i].size()>=n-1)solve(i);
	int res=0;
	for(int i=1;i<=lim;i++)Add(res,mul(f[i],i));
	cout<<res<<'\n';return 0;
}
stargazer.
关注
  • 1
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
2020美赛A题解题思路（Moving North）
01-06
文前推广 《计算机保研经历分享》 —— 大学期间参加了有含金量的比赛，该如何利用自身的优势保研、考研或求职？或者仅仅就是为了做一份好看的简历“以备不时之需”？ —— 关注我的文章，我会后期开始更新文章，以...
联合省选A卷2020总结
EasternCountry的博客
07-04 363
比赛前 我们学校比较特殊，初三的老师狠抓文化课，导致比赛前一天还在做语文试卷…… Day1 很紧张，上来第一题就是长长长的题面……看了大概30分钟，猛地发现这原来是一道线段树上分治的简单题？ 第一题打完后大概10:00 看第二题，送了40分，后面展开后发现ksk^sks怎么也展不开，想起以前学斯特林数的时候好像碰到过这个展开，但是已经忘却，开始找规律，最后什么也没有找到。 看第三题，什么想法都没有 考试出来后大家的估分大部分都是140.感觉好像还没有凉的太厉害。 下午赶快把斯特林数补了补，发现第二题就是一道
2020联合省选总结
Deep_Thinking的博客
06-22 673
结束了，这么久以来的努力 Day1 T1：这不是裸的数据结构，线段树上二分吗，2e6???，出题人怕不是脑子有坑，专门卡常？？？打了3K多，极限数据要跑5s，输入输出几千万个字符，我当场裂开。 T2：感觉m这么小肯定是推导它的组合意义，想了很久背包什么的东西想直接表示出来，并没有结果。也曾想过斯特林数拆掉k的次幂，但是不记得长什么样了！！！！！！！40分灰溜溜地走人。 T3：怎么拆线性基啊啊啊啊啊！！！我一直在想怎么记录选择的线性基，最后也是暴力走人。 感觉丝毫没有竞争力，T1也没有高明的树状数组二分.
题解 [联合省选 2020 A] 作业题（LOJ #3304 / 洛谷 P6624）【矩阵树定理】
破壁人五号的博客
06-24 395
题目链接：洛谷 P6624 / LOJ #3304 题意 给定无向带权简单图，其生成树的价值定义为树内边权之和乘以其 gcd⁡\gcdgcd，求其所有生成树的价值和。n≤30n\leq 30n≤30，w（边权）≤152501w{\tiny\text{（边权）}} \leq 152501w（边权）≤152501。 题解 考虑当 gcd⁡\gcdgcd 确定时所有生成树的边权和之和，然后把它转化为 gcd⁡\gcdgcd 为 ddd 的倍数时所有生成树的边权和之和再容斥一下。 先枚举 ddd，把所有 ddd 的
2020统一省选A卷题解
Master.Yi的博客
06-22 341
T2 ∑k=0nf(k)∗xk∗(nk)=∑k=0n∑i=0maikixk(nk)=∑i=0mai∑k=0nkixk(nk)=∑i=0mai∑k=0n∑j=0i{ij}(kj)j!∗xk(nk)=∑i=0mai∑j=0i{ij}j!∑k=0n(kj)xk(nk)=∑i=0mai∑j=0i{ij}j!∑k=0n(nj)(n−jk−j)xk=∑i=0mai∑j=0i{ij}j!(nj)∑k=0n(n−jk−j)xk=∑i=0mai∑j=0i{ij}j!(nj)xj∑k=0n−j(n−jk)xk=∑i=0mai∑
【2020省选模拟】题解
Stargazer的博客
04-02 205
NoipNoipNoip模拟？ T1 直接树剖即可 维护各凭本事 我是维护的轻儿子之和，所以细节要多很多 直接维护所有儿子之和似乎好做一些 标记永久化了，但似乎还是跑最慢。。。 #include<bits/stdc++.h> using namespace std; #define cs const #define pb push_back #define pii pair<i...
2020CCPC网络赛题解.pdf
07-12
2020CCPC网络赛题解
2016山东省NOI省队集训DAY6 题目题解数据标程 NOIP省选集训
04-18
2016山东省NOI省队集训DAY6 题目题解数据标程 NOIP省选集训
电赛-2020电赛A题代码整理-优质题解.zip
最新发布
02-10
电赛_2020电赛A题代码整理_优质题解.zip
2020美赛A题解题思路方法：向北运动
01-07
2020 MCM周末1 问题A：向北移动 全球海洋温度影响某些海洋生物的栖息地质量。 当温度变化太大而不利于它们继续繁衍时，这些物种就开始寻找 其他更适合他们现在和将来生活和繁殖成功的栖息地。一个 美国缅因州的龙虾...
【游记】联合省选2022
七代目鸽影
04-22 977
既退役，便不再去想。
联合省选2022游记
StaroForgin的博客
04-18 1642
这就是我们的现实。
2020 联合省选 杂题选讲（基本都不会）
woshidalaocxy的博客
06-23 928
B卷 D1T1 随便口胡了一个贪心就A了。 题目意思不难理解，这题目直接就能想到贪心求解，然后随便口胡了一个算法： 比如说我们现在手里拿到的牌的数是xxx，位置为ppp，则我们要求最小的rrr，满足x+∑k=p+1rakx+\sum_{k=p+1}^ra_kx+∑k=p+1r​ak​。大概的口胡证明如下： 假设我们现在有个牌数为AAA，旁边的牌为1,21,21,2，即A,1,2A,1,2A,1,2（A∈NA \in \mathbb{N}A∈N） 如果我们取A,1A,1A,1，那么就变成A+1,2A+1,2A
NOI2022联合省选 题解
liuzhangfeiabc的博客
04-17 2573
NOI2022联合省选 题解 目前只有day2
2022 联合省选游记 & 信息奥赛五年退役记
04-18 1605
再见我的青春，再见我的热爱。
联合省选2021总结——彻底的失败
wangyuda123456789的博客
04-15 355
两天的比赛策略出现了严重的问题。 该切的Day1T1没有切掉 没写对拍为了Day1T3暴力 一堆脑残错误送我上天，最后没时间了 Day2就是菜，没什么好说的 再见
省选联考 2020 A 卷
lahlah的博客
11-30 288
省赛虽然炸了，但题解还是要写滴（笑） 第200篇blog耶 组合数问题 题意简单明了 首先要看出要把多项式的每一项分别算 第ccc项的答案就是 ac∗∑k=0nkc×xk×Cnk\large a_c*\sum\limits_{k=0}^n k^c \times x^k \times C_{n}^{k}ac​∗k=0∑n​kc×xk×Cnk​ 先把aca_cac​扔了，看后面的那个式子，发现ccc比较小 考虑第二类斯特林数nm=∑i=0mS(m,i)×n!(n−i)!n^m=\sum\limits_{i=0}
联合省选2022复习计划
jwg2732的专栏
04-07 469
尽量通过做题的方式复习 已复习： 莫反 杜教筛 李超线段树 网络流 线段树扫描线 线段树合并 决策单调性分治 树形DP 线性基 根号分治 状压dp 组合
2020全国统一省选day1 魔法商店
A1847225889的博客
06-27 462
题目 正解 据说是一道论文题…… 论文：2018集训队论文高睿泉《浅谈保序回归问题》 保序回归问题： 有一个正整数ppp，给出一个有向无环图，点iii有权值(wi,yi)(w_i,y_i)(wi​,yi​)，需要调整yiy_iyi​的值使得yiy_iyi​满足有向无环图的偏序关系。调整的代价为前后yiy_iyi​的差的ppp次方乘wiw_iwi​，求最小的代价。 形式化地说：给每个点赋一个新的权值fif_ifi​，使得每条边(u,v)∈E(u,v)\in E(u,v)∈E满足fu≤fvf_u\leq f
洛谷p1036 选数python题解
08-31
对于洛谷上的p1036题目，我们可以使用Python来解决。下面是一个可能的解法： ```python def dfs(nums, target, selected_nums, index, k, sum): if k == 0 and sum == target: return 1 if index >= len(nums) ...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值