2017北京国庆刷题Day4 afternoon
期望得分:100+100+0=200
实际得分:5+0+0=5
每加入一个数,x的因数位置++
注意:根号x枚举时,如果x是完全平方数,根号x会重复累计2次,要减去
考场上没减,5分 /(ㄒoㄒ)/~~
#include<cmath>
#include<cstdio>
#include<iostream>
using namespace std;
#define N 40001
int sum[N];
void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
}
int main()
{
freopen("a.in","r",stdin);
freopen("a.out","w",stdout);
int n,ans=;
read(n);
int x,y,z;
while(n--)
{
read(x); read(y);
if(x==)
{
z=sqrt(y);
for(int i=;i<=z;i++)
if(y%i==) sum[i]++,sum[y/i]++;
if(z*z==y) sum[z]--;
}
else ans^=sum[y];
}
printf("%d",ans);
}
树形DP
令col[i] 表示 i与i的父节点之间的连边
f[i] 表示 在节点i的子树中,链的一个端点为i,且与i直接相邻的边的边的颜色 不与col[i] 相同 的方案数
F[i] 表示对应的权值和
g[i] 表示在节点i的子树中,链的一个端点为i 的方案数
G[i] 表示对应的权值和
将i的子树中经过i的链分为两种
① 以i为端点,用 g G 统计
②i为中间的一个点,用 f F 统计
状态转移方程:
设j为i的子节点,val[i]表示节点i的点权
① g[i]=Σ(f[j]+1) G[i]=Σ(F[j]+val[j])
② f[i]= Σ(f[j]+1) F[i]=Σ(F[j]+val[j]) 其中,j满足col[j]!=col[i]
③ G[i]+=g[i]*val[i]
④ F[i]+=f[i]*val[i]
解释:
1、g[i]、f[i] 的转移 就是把子树中的 g,f 累加起来,再加上子树的个数,因为i和i的子树的根节点(即i的子节点)构成一条新的合法的链
2、①、② 中 G、F 的转移 就是把子树中的 G、F累加起来,再加上子节点的权值和,因为i和的i的子节点构成一条新的合法的链
3、③、④ 所有的以i为一个端点的链都累加一个 i的权值
统计答案:
1、以i为端点的链,就是ans+=G[i]
2、以i为中间一个点的链,显然是要拿以i为端点的两条链拼起来
①对于i的每个子树j,假设j会被使用sum次,子树j的合法链的权值总和为V,那么 这个子树j 对 答案的贡献就是 sum*(V+val[j])。
加val[j]是因为 子树j的根节点的权值不属于V,但以这个点为链的一个端点,以i的其他子树的一个点为链的另一个端点, 这就是一条合法的链
如何统计sum?
设s[k]表示当前i的子树中,以i为链的一个端点 且 与i直接相连的边的颜色为k 的链的条数
颜色编号大至1e9? ——离散化
那么sum=g[i]-s[col[j]] 即只要与i直接相连的边的颜色 不等于 col[j] ,就可以与j的子树中的链 以及 j 拼接
② 考虑了链的拼接,还差两条链的交点的权值没有加,即i的权值。
设两条链拼接的总方案数位 tot,那么最后再加上 tot*val[i] 就行了
tot=Σ ( (f[[j]+1) * (g[i]-s[col[j]]) ) / 2 原理同上
除2是因为 一条链被枚举了两次
#include<cstdio>
#include<iostream>
#include<algorithm> using namespace std; #define N 300001 typedef long long LL; int val[N];
int front[N],to[N<<],nxt[N<<],tot,col[N<<];
int has[N];
LL f[N],s[N],F[N];
LL ans; void read(int &x)
{
x=; char c=getchar();
while(!isdigit(c)) c=getchar();
while(isdigit(c)) { x=x*+c-''; c=getchar(); }
} void add(int u,int v,int w)
{
to[++tot]=v; nxt[tot]=front[u]; front[u]=tot; col[tot]=w;
to[++tot]=u; nxt[tot]=front[v]; front[v]=tot; col[tot]=w;
} void init()
{
int n; read(n);
for(int i=;i<=n;i++) read(val[i]);
int u,v,c;
for(int i=;i<n;i++)
{
read(u); read(v); read(c);
has[i]=c;
add(u,v,c);
}
sort(has+,has+n+);
int tot=unique(has+,has+n+)-(has+);
for(int i=;i<=n-<<;i++) col[i]=lower_bound(has+,has+tot+,col[i])-has;
} void dfs(int x,int y,int z)
{
int num=;
for(int i=front[x];i;i=nxt[i])
if(to[i]!=y) num++,dfs(to[i],x,col[i]);
LL g=,G=;
for(int i=front[x];i;i=nxt[i])
if(to[i]!=y)
{
if(col[i]!=z) f[x]+=f[to[i]]+,F[x]+=F[to[i]]+val[to[i]];
g+=f[to[i]]+,G+=F[to[i]]+val[to[i]];
s[col[i]]+=f[to[i]]+;
}
if(!num) return;
F[x]+=1ll*f[x]*val[x];
G+=1ll*g*val[x];
ans+=G;
LL res=,tmp=;
for(int i=front[x];i;i=nxt[i])
if(to[i]!=y) res+=(F[to[i]]+val[to[i]])*(g-s[col[i]]),tmp+=(+f[to[i]])*(g-s[col[i]]);
ans+=res+tmp/*val[x];
for(int i=front[x];i;i=nxt[i]) s[col[i]]=;
} int main()
{
freopen("b.in","r",stdin);
freopen("b.out","w",stdout);
init();
dfs(,,);
printf("%I64d\n",ans);
}
解决本题关键:
同一行/列只有被选中奇数次才有效
假设有i行j列被翻了过来
那么可以得到等式
i*m+j*n-2*i*j=s
解得j=(s-i*m)/(n-2*i)
由此可知,我们只需要枚举i,就可以直接算出j
这个 i,j 合法的条件是:
① 不越界
②(n-i)%2=0,(m-j)%2=0
因为只能再翻偶数次,才能保证当前i,j 合法
如何计算一对合法的i,j的答案?
n行里i行被翻了过来 C(n,i)
m列里j列被翻了过来 C(m,j)
被翻了偶数次的行,就是把(r-i)/2 次机会 分给 n 行 C((r-i)/2+n-1,n-1)
注意不是n行里面选 (r-i)/2 行翻过来,因为同一行可以不翻,也可以翻多次
同理,偶数次列为 C((c-i)/2+m-1,m-1)
把这4个C 乘起来就是这一对i,j 的答案
最后累加所有的i,j 的贡献即可
#include<cstdio>
#include<algorithm> using namespace std; #define N 200001 typedef long long LL; const int mod=1e9+; LL inv[N],fac[N]; LL pow(LL a,int b)
{
LL res=;
for(;b;a=a*a%mod,b>>=)
if(b&) res*=a,res%=mod;
return res;
} int main()
{
freopen("c.in","r",stdin);
freopen("c.out","w",stdout);
int n,m,r,c; LL s;
scanf("%d%d%d%d%I64d",&n,&m,&r,&c,&s);
fac[]=;inv[]=;
for(int i=;i<=max(n,m)+max(r,c);i++) fac[i]=i*fac[i-],fac[i]%=mod,inv[i]=pow(fac[i],mod-);
int j;
int ans=; LL res;
for(int i=(r&);i<=min(r,n);i+=)
{
if(n!=*i)
{
if((s-1ll*i*m)%(n-*i)) continue;
j=(s-1ll*i*m)/(n-*i);
if(j< || j>min(c,m) || (c-j)&) continue;
res=fac[n]*inv[i]%mod*inv[n-i]%mod;
res=res*fac[m]%mod*inv[j]%mod*inv[m-j]%mod;
res=res*fac[n+(r-i>>)-]%mod*inv[n-]%mod*inv[r-i>>]%mod;
res=res*fac[m+(c-j>>)-]%mod*inv[m-]%mod*inv[c-j>>]%mod;
ans+=res; ans%=mod;
//printf("%d %d %I64d\n",i,j,res);
}
}
printf("%d",ans);
}
2017北京国庆刷题Day4 afternoon的更多相关文章
- 2017北京国庆刷题Day1 afternoon
期望得分:100+100+100=300 实际得分:100+100+100=300 T1 一道图论好题(graph) Time Limit:1000ms Memory Limit:128MB 题目 ...
- 2017北京国庆刷题Day5 afternoon
期望得分:100+60+100=260 实际得分:0+60+40=100 设图中有m个环,每个环有si条边,有k条边不在环中 ans= (2^s1 -2)*( 2^s2 -2)* (2^s3 -2)… ...
- 2017北京国庆刷题Day2 afternoon
期望得分:100+100+50=250 实际得分:100+70+50=220 T1 最大值(max) Time Limit:1000ms Memory Limit:128MB 题目描述 LYK有一 ...
- 2017北京国庆刷题Day4 morning
期望得分:0+40+30=70 实际得分:0+10+10=20 题目修改:只能由0变1,只能用一次操作 大模拟 #include<cstdio> #include<cstring&g ...
- 2017北京国庆刷题Day6 afternoon
期望得分:100+100+40=240 实际得分:100+0+40=140 二进制拆分.二进制前缀和 #include<cstdio> #include<iostream> u ...
- 2017北京国庆刷题Day3 afternoon
期望得分:100+0+30=130 实际得分:100+36.5+0=136.5 T3 一个变量写混了,丢了30.. 模拟栈 #include<cstdio> #include<cst ...
- 2017北京国庆刷题Day7 afternoon
期望得分:100+30+100=230 实际得分:60+30+100=190 排序去重 固定右端点,左端点单调不减 考场上用了二分,没去重,60 #include<cstdio> #inc ...
- 2017北京国庆刷题Day7 morning
期望得分:100+0+100=200 实际得分:100+20+0=120 离散化搞搞 #include<cstdio> #include<iostream> #include& ...
- 2017北京国庆刷题Day3 morning
期望得分:100+60+0=160 实际得分:100+30+0=130 考场上用的哈希 #include<cstdio> #include<cstring> #include& ...
随机推荐
- "感应锁屏"Alpha版使用说明
“感应锁屏”Alpha版使用说明 一.产品介绍 感应锁屏是锁屏软件的一大创新.相比传统的锁屏软件,“感应锁屏”从可操作性.方便性.功能全面性都有了很大的提升,可以让用户方便快捷的进行锁屏操作. “感应 ...
- RovingUI组件库-包含堆栈式通知提醒框(Toast)的小程序组件库
RovingUI是个人在开发小程序过程中将用到的组件集合而成的一个UI库,包含一些基本通用组件(按钮.栅格.通用样式.徽标.通知和面包屑). 源起得归于我在开发中没有找到现成的堆栈式提醒框(比如ant ...
- Spark Shuffle之Sort Shuffle
源文件放在github,随着理解的深入,不断更新,如有谬误之处,欢迎指正.原文链接https://github.com/jacksu/utils4s/blob/master/spark-knowled ...
- 使用docker国内镜像解决方案
1:蜂巢镜像 https://c.163yun.com/hub#/m/library/ 例如: docker pull hub.c.163.com/library/nginx:1.8 再次执行dock ...
- Qt多线程-QThread
版权声明:若无来源注明,Techie亮博客文章均为原创. 转载请以链接形式标明本文标题和地址: 本文标题:Qt多线程-QThread 本文地址:http://techieliang.com/2 ...
- 关于初装kali linux 2.0时DEB文件安装失败的问题
kali linux 是一个基于debian 的linux发行版本,支持deb文件格式的图形化安装. 刚装上kali linux时安装程序总是失败,提示处理时错误. 经过一番爬贴,是软件源的原因,解决 ...
- [cnbeta]微软最强数据中心级操作系统
微软近日发表了一篇介绍Windows系统内核的博文,期间为了展示Windows的强大扩展性,放出了一张非常震撼的Windows任务管理器截图:乍一看似乎没啥特别的,几十甚至上百个逻辑核心的系统并不罕见 ...
- Redis 备份数据的两种方式
既然是数据库,那就一定有数据备份方式了,而且 Redis 是内存形式的数据库,更需要数据备份了,要不然断电数据就全都丢失了. Redis 数据备份有两种方式: RDB(数据快照) AOF(记录操作日志 ...
- Delphi ADOQuery多个参数重复 改编技巧
今天看了多年前的一个帖子,发现回答不合理,有些还将其归为delphi的bug.其实主要是没有灵活应用参数. ADOQUERY查询时,这样不行,结果不正确. WITH ADOQUERY1 DO BEGI ...
- SQL优化套路
直奔主题: (1) SELECT(2) DISTINCT <SELECT_LIST>(3) FROM <LEFT_TABLE> <JOIN_TYPE> (4) JO ...