https://vjudge.net/contest/321565#problem/C 超时代码
2024-10-07 21:35:19
#include <iostream>
#include <cstdio>
#include <queue>
#include <algorithm>
#include <cmath>
#include <cstring>
#define inf 2147483647
#define N 10100
#define p(a) putchar(a)
#define For(i,a,b) for(register long long i=a;i<=b;++i)
//by war
//2019.8.22
using namespace std;
long long T,n,x,y,cnt,tot;
long long prime[N],mu[N],ans[];
bool vis[N]; struct dian{
long long l,r,t;
}a[N]; struct node{
long long n;
node *next;
}*e[N]; inline void in(long long &x){
long long y=;char c=getchar();x=;
while(c<''||c>''){if(c=='-')y=-;c=getchar();}
while(c<=''&&c>=''){ x=(x<<)+(x<<)+c-'';c=getchar();}
x*=y;
}
inline void o(long long x){
if(x<){p('-');x=-x;}
if(x>)o(x/);
p(x%+'');
} inline void push(long long x,long long y){
node *p;
p=new node();
p->n=y;
if(e[x]==)
e[x]=p;
else{
p->next=e[x]->next;
e[x]->next=p;
}
} void Euler(){
mu[]=;
For(i,,){
if(!vis[i]) prime[++cnt]=i,mu[i]=-;
for(register long long j=;j<=cnt&&i*prime[j]<=;j++){
vis[i*prime[j]]=;
if(i%prime[j]==){
mu[i*prime[j]]=;
break;
}
mu[i*prime[j]]=-mu[i];
}
}
} long long dfs(long long x,long long fa,long long w){
long long res=;
for(node *i=e[x];i;i=i->next)
if(i->n!=fa)
res+=dfs(i->n,x,w*a[i->n].t);
return res+w;
} inline long long F(register long long d){
long long res=,kk=;
For(i,,n){
a[i].t=a[i].r/d-a[i].l/d;
if(a[i].l%d==)
a[i].t++;
}
For(i,,n){
res+=dfs(i,i,a[i].t);
res-=a[i].t;
kk+=a[i].t;
}
return res/+kk;
} inline void clear(){
For(i,,)
e[i]=;
memset(ans,,sizeof(ans));
} signed main(){
in(T);
Euler();
while(T--){
clear();
in(n);
For(i,,n-){
in(x);in(y);
push(x,y);
push(y,x);
}
For(i,,n) in(a[i].l);
For(i,,n) in(a[i].r);
For(i,,)
for(register long long d=i;d<=;d+=i)
ans[i]+=mu[d/i]*F(d);
printf("Case %lld:\n",++tot);
For(i,,){
o(i);p(':');p(' ');
o(ans[i]);p('\n');
}
}
return ;
}
最新文章
- LogBack,升级版的log4J
- 使用rest方式修改服务端xml文件
- ajax返回数据类型为XML数据的处理
- [solr] - Facet - autocomplete
- Linux 学习碎片
- web.xml配置
- Dire Wolf ---hdu5115(区间dp)
- 设计模式 适配器-Adapter
- ili9341 横屏驱动代码
- 安装 php 转
- MongoDB源码分析——mongod数据查询操作
- ubuntu 下的 ftp (gftp)
- Android源码中的FLAG为何使用16进制
- appendChild的用法
- php的命名空间层级与目录层级是一致的吗?
- Activity的Task详解
- efcore 配置链接sqlserver 记录
- BZOJ:4816: [Sdoi2017]数字表格
- php 定时任务
- pythonのsimple_tag