2019银川F,ccpc威海D - Sternhalma 2022
创始人
2024-01-20 18:31:32
0

1401D - Maximum Distributed Tree

求每个边经过的次数,假设求u,v这条边的次数,边的左端是u这个集合一共有n-siz[v]个点,右端是v这个集合有siz[v]个端点,经过这条边的次数就是siz[v]*(n-siz[v]),然后再按照次数多的乘以大的质因数就可以了,注意m可能大于n-1

D. Maximum Distributed Tree(贪心+树dfs)_小菜鸡加油的博客-CSDN博客

#include 
using namespace std;
#define endl '\n'
#define pause system("pause")
#define int long long
const int mod=1e9+7;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-10;int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)b;}
void dfs(int u,int fa)
{siz[u]=1;for(int i=head[u];i;i=e[i].next){int j=e[i].to;if(j==fa) continue;dfs(j,u);siz[u]+=siz[j];a[++ct]=siz[j]*(n-siz[j]);}
}
signed main()
{//ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);//freopen("in.txt","r",stdin);cin>>t;while(t--){cin>>n;for(int i=1;i<=n;i++) siz[i]=head[i]=0;cnt=ct=0;for(int i=1;i>u>>v;addedge(u,v);addedge(v,u);}dfs(1,0);cin>>m;for(int i=1;i<=m;i++) cin>>p[i];int ans=0;if(m<=ct){sort(a+1,a+ct+1,cmp);sort(p+1,p+m+1,cmp);for(int i=m+1;i<=ct;i++) p[i]=1;for(int i=1;i<=ct;i++){ans=(ans+(p[i]*a[i]%mod))%mod;//cout<=1;i--){ans=(ans+(p[i]*a[i]%mod))%mod;// cout<

F. Function! 2019银川,类似整除分块

因为当b>a的时候,\log_{b}{a}都是小于1的,上取整之后就是1,所以整个式子就变成

\sum_{a}^{n}a\sum_{b=a}^{n}\log_{a}{b},当a>\sqrt{n}时,1<=\log_{a}{b}<2=>\log_{a}{b}=1,所以右边的求和其实就是(n-a+1),这玩意是可以化简得,\sum_{a}^{n}a(n-a+1)=\sum_{a}^{n}(n+1)a-a^{2}=(n+1)\sum_{a}^{n}a-\sum_{a}^{n}a^{2},

一个是等差数列求和,一个是平方和,这就可以o(1)得算出来了;

然后a<=\sqrt{n}时直接暴力算,但发现对于一个b,会有一段连续的a log值时一样的,所以可以利用类似整除分块的思想来优化一下;

2019ICPC(银川) - Function!(数论+数学分块)_Frozen_Guardian的博客-CSDN博客

#include 
using namespace std;
#define endl '\n'
#define pause system("pause")
#define int long long
const int mod=998244353;
const int inf=1e18;
const int N = 4e5+100;
const double eps=1e-10;int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)>n;int ans=0,a;for(a=2;a*a<=n;a++){int tmp=0;int qp=a;int x=1;for(int b=a;b<=n;b++){int d=min(n,qp*a-1LL);tmp=(tmp+x*((d-b+1)%mod)%mod)%mod;b=d;qp*=a;x++;//cout<

D - Sternhalma 2022ccpc威海

一共就19个格子,并且每个格子的权值是不会变的,所以可以记忆化加状压,这题就是一个带状压的记忆化搜索,但是实现雀氏有点难,直接看代码就可以

2022CCPC威海站 铜牌题解 A C D E G I J - 知乎 (zhihu.com)

#include 
using namespace std;
#define endl '\n'
#define lowbit(x) ((x)&(-x))
#define int long long
#define pause system("pause")
const int mod=998244353;
const int inf=1e18;
const int N = 1e6+100;
const double eps=1e-10;int qpow(int a,int b)
{int res=1;while(b){if(b&1) res=res*a%mod;a=a*a%mod;b>>=1;}return res;
}
int sgn(double x)
{if(fabs(x)>coord=
{{1,3},{1,5},{1,7},{2,2},{2,4},{2,6},{2,8},{3,1},{3,3},{3,5},{3,7},{3,9},{4,2},{4,4},{4,6},{4,8},{5,3},{5,5},{5,7}
};
int s[10][10],id[10][10],vis[N],f[N],n;
int tran(string s)
{int res=0;for(int i=0;i>i)&1;if(x==0) continue;int nstate=state&(~(1<>i&1;}for(int i=0;i<19;i++){if((state>>i&1)==0) continue;auto [x,y]=coord[i];for(int j=0;j<6;j++){int ax=x+d[j][0],ay=y+d[j][1];int bx=x+d2[j][0],by=y+d2[j][1];if(ax<0||ay<0||bx<0||by<0) continue;if(id[ax][ay]==-1||id[bx][by]==-1) continue;if(g[ax][ay]==0||g[bx][by]==1) continue;int nstate=state;g[ax][ay]=g[x][y]=0;g[bx][by]=1;nstate=nstate&(~(1<>s[x][y];}vis[0]=1;f[0]=0;cin>>n;for(int i=1;i<=n;i++){string t="",g;for(int j=1;j<=5;j++) cin>>g,t+=g;int ans=dfs(tran(t));cout<

相关内容

热门资讯

喜欢穿一身黑的男生性格(喜欢穿... 今天百科达人给各位分享喜欢穿一身黑的男生性格的知识,其中也会对喜欢穿一身黑衣服的男人人好相处吗进行解...
发春是什么意思(思春和发春是什... 本篇文章极速百科给大家谈谈发春是什么意思,以及思春和发春是什么意思对应的知识点,希望对各位有所帮助,...
网络用语zl是什么意思(zl是... 今天给各位分享网络用语zl是什么意思的知识,其中也会对zl是啥意思是什么网络用语进行解释,如果能碰巧...
为什么酷狗音乐自己唱的歌不能下... 本篇文章极速百科小编给大家谈谈为什么酷狗音乐自己唱的歌不能下载到本地?,以及为什么酷狗下载的歌曲不是...
华为下载未安装的文件去哪找(华... 今天百科达人给各位分享华为下载未安装的文件去哪找的知识,其中也会对华为下载未安装的文件去哪找到进行解...
怎么往应用助手里添加应用(应用... 今天百科达人给各位分享怎么往应用助手里添加应用的知识,其中也会对应用助手怎么添加微信进行解释,如果能...
家里可以做假山养金鱼吗(假山能... 今天百科达人给各位分享家里可以做假山养金鱼吗的知识,其中也会对假山能放鱼缸里吗进行解释,如果能碰巧解...
一帆风顺二龙腾飞三阳开泰祝福语... 本篇文章极速百科给大家谈谈一帆风顺二龙腾飞三阳开泰祝福语,以及一帆风顺二龙腾飞三阳开泰祝福语结婚对应...
美团联名卡审核成功待激活(美团... 今天百科达人给各位分享美团联名卡审核成功待激活的知识,其中也会对美团联名卡审核未通过进行解释,如果能...
四分五裂是什么生肖什么动物(四... 本篇文章极速百科小编给大家谈谈四分五裂是什么生肖什么动物,以及四分五裂打一生肖是什么对应的知识点,希...