Luogu 2671 求和 NOIP2015T3
题解
20pts
$O(n^3)$枚举$x,y,z$,根据题目要求判断
40pts
$O(n^2)$枚举$x,z$,需要满足$x,z$奇偶相同
20~40pts的代码我都没有写过...就不贴了
70~90pts
尝试对40pts的暴力进行优化
题目对三元组的两个限制我们在40pts的时候已经用来优化过了,还有一个限制是颜色相同
我们可以以颜色为关键字对数组进行排序,对每个数算答案的时候,只需要枚举相同颜色的这一段即可(因为三元组要求有序,所以要从$x+2$开始枚举)
这样的复杂度是$O(n*cnt_{col})$(这里的$cnt_{col}$为出现次数最多的颜色的出现次数)
使用$scanf$将会得到70pts
使用快读或$scanf+O2$将会得到80pts
使用快读+$O2$将会得到90pts
($fread$在这里的效果和快读是差不多的,因为数只有$1e5$个)
#include <bits/stdc++.h> #define ll long long
#define inf 0x3f3f3f3f
#define il inline namespace io { #define in(a) a=read()
#define out(a) write(a)
#define outn(a) out(a),putchar('\n') #define I_int int
inline I_int read() {
I_int x = , f = ; char c = getchar() ;
while( c < '' || c > '' ) { if( c == '-' ) f = - ; c = getchar() ; }
while( c >= '' && c <= '' ) { x = x * + c - '' ; c = getchar() ; }
return x * f ;
}
char F[ ] ;
inline void write( I_int x ) {
I_int tmp = x > ? x : -x ;
if( x < ) putchar( '-' ) ;
int cnt = ;
while( tmp > ) {
F[ cnt ++ ] = tmp % + '' ;
tmp /= ;
}
while( cnt > ) putchar( F[ -- cnt ] ) ;
}
#undef I_int }
using namespace io ; using namespace std ; #define N 100010
const int mod = ; int n = read() , m = read() ;
struct node {
int col , val , id ;
} a[ N ] ; bool cmp( node a , node b ) {
return a.col < b.col ;
} int main() {
for( int i = ; i <= n ; i ++ ) a[ i ].val = read() , a[ i ].id = i ;
for( int i = ; i <= n ; i ++ ) a[ i ].col = read() ;
sort( a + , a + n + , cmp ) ;
int cur = ;
ll ans = ;
for( int i = ; i <= n ; i ++ ) {
ll sum = ;
for( cur = i + ; a[ i ].col == a[ cur ].col && cur <= n ; cur ++ ) {
if( ( a[ i ].id + a[ cur ].id ) % == )
sum = ( sum + ( 1ll * ( a[ i ].id + a[ cur ].id ) % mod * 1ll * ( a[ i ].val + a[ cur ].val ) % mod ) % mod ) % mod ;
}
ans = ( ans + sum ) % mod ;
}
printf( "%lld\n" , ( ans + mod ) % mod ) ;
return ;
}
100pts
考虑使用数学方法优化以上做法
我们现在有什么条件呢,列举一下
1.$x$和$z$的奇偶性相同且颜色相同
2.求和公式为$(x+z)(num_x+num_z)$
从求和公式入手,把括号拆掉,式子变成$x*num_x+x*num_z+z*num_x+z*num_z$
设$x1,x2,x3$的奇偶性相同且颜色相同(那么他们组成的二元组就等同于符合条件的三元组)
考虑$x1$对答案的贡献:
对于$(x1,x2)$,$x1$的贡献为$x1*num_{x1}+x1*num_{x2}$
对于$(x1,x3)$,$x1$的贡献为$x1*num_{x1}+x1*num_{x3}$
那么$x_1$对答案的贡献就是
$x1*(num_{x2}+num_{x3})+2*x1*num_{x1}$
然后多加入一个数$x4$也可以得到类似的贡献(多代几个也就看出来规律了)
推广到$x_n$结论也是一样的
结论:
设$s=\sum num_{xi}$(这里的$xi$均满足条件1)$cnt=cnt_{col}$,$cnt_{col}$为当前颜色奇偶相同的数的个数
则$xi$对答案的贡献为$xi*(s-num_{xi})+(cnt-1)*xi*num_{xi}$($cnt$要$-1$,因为$xi$不能和自己组成二元组)
所以将$s$和$cnt$统计出来就行了,注意要分奇偶统计
复杂度$O(n)$
#include <bits/stdc++.h> #define ll long long
#define inf 0x3f3f3f3f
#define il inline
#define int long long namespace io { #define in(a) a=read()
#define out(a) write(a)
#define outn(a) out(a),putchar('\n') #define I_int int
inline I_int read() {
I_int x = , f = ; char c = getchar() ;
while( c < '' || c > '' ) { if( c == '-' ) f = - ; c = getchar() ; }
while( c >= '' && c <= '' ) { x = x * + c - '' ; c = getchar() ; }
return x * f ;
}
char F[ ] ;
inline void write( I_int x ) {
if( x == ) { putchar( '' ) ; return ; }
I_int tmp = x > ? x : -x ;
if( x < ) putchar( '-' ) ;
int cnt = ;
while( tmp > ) {
F[ cnt ++ ] = tmp % + '' ;
tmp /= ;
}
while( cnt > ) putchar( F[ -- cnt ] ) ;
}
#undef I_int }
using namespace io ; using namespace std ; #define N 100010
const int mod = ; int n , m ;
struct node {
int col , val ;
} a[ N ] ;
int cnt[ N ][ ] , sum[ N ][ ] ; signed main() {
n = read() , m = read() ;
for( int i = ; i <= n ; i ++ ) a[ i ].val = read() ;
for( int i = ; i <= n ; i ++ ) a[ i ].col = read() ;
for( int i = ; i <= n ; i ++ ) {
cnt[ a[ i ].col ][ i& ] ++ ;
sum[ a[ i ].col ][ i& ] = ( sum[ a[ i ].col ][ i& ] + a[ i ].val ) % mod ;
}
ll ans = ;
for( int i = ; i <= n ; i ++ ) {
ll num = cnt[ a[ i ].col ][ i& ] , s = sum[ a[ i ].col ][ i& ] ;
ans = ( ans + ( i*(s-a[i].val) + (num-)*a[i].val*i ) % mod ) % mod ;
}
printf( "%lld\n" , ans % mod ) ;
return ;
}
Luogu 2671 求和 NOIP2015T3的更多相关文章
- luogu 4427 求和
bjoi 2018 求和 唯一一道可能切的题一个数组还没开long long就成0分了 题目大意: 一棵有根树,并且希望多次询问这棵树上一段路径上所有节点深度的k次方和,而且每次的k可能是不同的 此处 ...
- [洛谷2671]求和<前缀和&模拟>
题目链接:https://www.luogu.org/problemnew/show/P2671 这是noip2015普及组的第三题,谁说的普及组的题就一定水的不行,这道题就比较有意思的 这道题的暴力 ...
- Luogu P1625 求和
题意 给定两个整数 \(n,m\),求 \[\sum\limits_{i=1}^{n}\frac{1}{\prod\limits_{j=i}^{i+m-1}j} \] \(\texttt{Data R ...
- luogu P1630 求和(枚举暴力)
题意 题解 可以发现当a=10001时, 和1是等价的. 所以这题就水了. #include<iostream> #include<cstring> #include<c ...
- [Luogu] 余数求和
question: $$\sum_{i=1}^{n} k \bmod i$$$$\sum_{i=1}^{n} k - \lfloor \frac{k}{i} \rfloor i$$$$\sum_{i= ...
- [Luogu 2261] CQOI2007 余数求和
[Luogu 2261] CQOI2007 余数求和 这一定是我迄今为止见过最短小精悍的省选题了,核心代码 \(4\) 行,总代码 \(12\) 行,堪比小凯的疑惑啊. 这题一看暴力很好打,然而 \( ...
- BZOJ 4555 Luogu P4091 [HEOI2016/TJOI2016]求和 (第二类斯特林数)
题目链接 (luogu) https://www.luogu.org/problem/P4091 (bzoj) https://www.lydsy.com/JudgeOnline/problem.ph ...
- [Luogu P2261] [CQOI2007]余数求和 (取模计算)
题面 传送门:https://www.luogu.org/problemnew/show/P2261 Solution 这题显然有一个O(n)的直接计算法,60分到手. 接下来我们就可以拿出草稿纸推一 ...
- 【Bzoj4555】【Luogu P4091】求和(NTT)
题面 Bzoj Luogu 题解 先来颓柿子 $$ \sum_{i=0}^n\sum_{j=0}^iS(i,j)2^jj! \\ =\sum_{j=0}^n2^jj!\sum_{i=0}^nS(i,j ...
随机推荐
- 003-redis-命令-key操作,字符串操作
Redis 键(key) Redis 键命令用于管理 redis 的键. 序号 命令及描述 1 DEL key该命令用于在 key 存在时删除 key. 2 DUMP key 序列化给定 key ,并 ...
- vue中使用better-scroll实现滑动效果
1.安装:npm install better-scroll 2.引入:import BetterScrol from "better-scroll"; 1.滚动效果 better ...
- http协议基础(八)请求首部字段
请求首部字段 定义:请求首部字段是从客户端到服务器发送请求报文中所使用的字段,里面包含了附加信息.客户端信息以及对响应内容相关的优先级等内容 1.Accept 通知服务器用户代理可处理的媒体类型及媒体 ...
- mysql日志详解
日志分类: 一.错误日志. 1.在配置文件中的配置是:log-error="DESKTOP-igoodful.err",查看参数的键值对:show variables like ' ...
- Python: 没有switch-case语句
初学Python语言,竟然很久才发现Python没有switch-case语句 官方的解释说,“用if... elif... elif... else序列很容易来实现 switch / case 语句 ...
- linux服务器---安装samba
安装samba 1.检测samba是否安装,如果没有,那么可以使用yum来安装.至少需要安装3个软件:samba,samba-client.samba-common [root@localhost p ...
- 安全测试工具之Burpsuite
端口即服务,每一个服务对应一个或多个端口.端口扫描即通过一些方法检测到一台主机的一段特定端口是否提供相应的服务.利用这些扫描结果,正常用户可以访问系统所提供的服务,而黑客却可以利用这些服务中的漏洞对系 ...
- python之路----TCP与UDP
TCP import socket #tcp协议 sk = socket.socket() # 买手机 创建一个socket对象 sk.bind(('127.0.0.1',8080)) # 给serv ...
- MIPI接口LCD屏调试心得(转)
源: MIPI接口LCD屏调试心得
- 为自己的网站添加Markdown功能 markedjs
Markdown几个简单的标记可以实现轻量级的代替Word方案 不多说,引入开源库js https://github.com/chjj/marked使用方式简单,如下实例代码: <!DOCTYP ...