博客
关于我
bzoj 2956: 模积和
阅读量:259 次
发布时间:2019-03-01

本文共 1174 字,大约阅读时间需要 3 分钟。

题意:求
nimj(n mod i)(m mod j)(i!=j)

题解:

我各种sb错误荒废了一上午啊啊啊。

先无视限制:即

ni(n mod i)mj(m mod j)
同bzoj 1257()
然后要减去相等的情况。
min(n,m)i(n mod i)(m mod i)
=min(n,m)i(nnii)(mmii)
=min(n,m)i(nm)(nim+min)+nimii2
然后就能分块了。

注意:乘一步mod一步啊啊啊!!!

code:

#include
#include
#include
#include
#define LL long longusing namespace std;const LL mod=19940417,inv=3323403;LL sum(LL l,LL r){ return ((LL)(r-l+1)*(LL)(l+r)/2)%mod;}LL sum_2(LL n){ return (LL)n*(LL)(n+1)%mod*(LL)(2*n+1)%mod*inv%mod;}LL solve(LL n){ LL ans=0;LL pos; for(LL i=1;i<=n;i=pos+1) { pos=(n/(n/i)); ans=(ans+((LL)n*(pos-i+1)-(LL)(n/i)*sum(i,pos))%mod)%mod; } return ans;}LL solve2(LL n,LL m){ if(n>m) swap(n,m); LL ans=0;LL pos; for(LL i=1;i<=n;i=pos+1) { pos=min(n/(n/i),m/(m/i)); LL t1=n*m%mod*(pos-i+1)%mod,t2=sum(i,pos)%mod*(n*(m/i)%mod+m*(n/i)%mod)%mod,t3=(m/i)*(n/i)%mod*(sum_2(pos)-sum_2(i-1)+mod)%mod; t3=(t3+mod)%mod; ans=((ans+(t1-t2+t3)%mod)+mod)%mod; } return ans;}LL n,m;int main(){ scanf("%lld %lld",&n,&m); printf("%lld",((solve(n)*solve(m)%mod-solve2(n,m))%mod+mod)%mod);}

转载地址:http://rdza.baihongyu.com/

你可能感兴趣的文章
Mysql学习总结(65)——项目实战中常用SQL实践总结
查看>>
Mysql学习总结(66)——设置MYSQL数据库编码为UTF-8
查看>>
Mysql学习总结(67)——MYSQL慢查询日志
查看>>
Mysql学习总结(68)——MYSQL统计每天、每周、每月、每年数据 SQL 总结
查看>>
Mysql学习总结(69)——Mysql EXPLAIN 命令使用总结
查看>>
Mysql学习总结(6)——MySql之ALTER命令用法详细解读
查看>>
Mysql学习总结(70)——MySQL 优化实施方案
查看>>
Mysql学习总结(71)——MySQL 重复记录查询与删除总结
查看>>
Mysql学习总结(71)——数据库介绍(MySQL安装 体系结构、基本管理)再回顾
查看>>
Mysql学习总结(72)——MySQL 开发者开发,设计规范再总结
查看>>
Mysql学习总结(73)——MySQL 查询A表存在B表不存在的数据SQL总结
查看>>
Mysql学习总结(74)——慢SQL!压垮团队的最后一根稻草!
查看>>
Mysql学习总结(75)——并发量大、数据量大的互联网业务数据库设计军规
查看>>
Mysql学习总结(76)——MySQL执行计划(explain)结果含义总结
查看>>
Mysql学习总结(77)——温故Mysql数据库开发核心原则与规范
查看>>
Mysql学习总结(78)——MySQL各版本差异整理
查看>>
Mysql学习总结(79)——MySQL常用函数总结
查看>>
Mysql学习总结(7)——MySql索引原理与使用大全
查看>>
Mysql学习总结(80)——统计数据库的总记录数和库中各个表的数据量
查看>>
Mysql学习总结(81)——为什么MySQL不推荐使用uuid或者雪花id作为主键?
查看>>