F - Expected Square Beauty

思路:
https://codeforces.com/blog/entry/68111

代码:

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#pragma GCC optimize(4)
#include<bits/stdc++.h>
using namespace std;
#define y1 y11
#define fi first
#define se second
#define pi acos(-1.0)
#define LL long long
//#define mp make_pair
#define pb push_back
#define ls rt<<1, l, m
#define rs rt<<1|1, m+1, r
#define ULL unsigned LL
#define pll pair<LL, LL>
#define pli pair<LL, int>
#define pii pair<int, int>
#define piii pair<pii, int>
#define pdd pair<double, double>
#define mem(a, b) memset(a, b, sizeof(a))
#define debug(x) cerr << #x << " = " << x << "\n";
#define fio ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
//head const int MOD = 1e9 + ;
const int N = 2e5 + ;
int l[N], r[N], n;
LL p[N], q[N], sum[N];
LL q_pow(LL n, LL k) {
LL res = ;
while(k) {
if(k&) res = (res * n) % MOD;
n = (n * n) % MOD;
k >>= ;
}
return res;
}
int main() {
scanf("%d", &n);
for (int i = ; i <= n; ++i) scanf("%d", &l[i]);
for (int i = ; i <= n; ++i) scanf("%d", &r[i]);
for (int i = ; i <= n; ++i) {
int L = max(l[i-], l[i]), R = min(r[i-], r[i]);
if(L > R) p[i] = , q[i] = ;
else q[i] = (R-L+)*q_pow((r[i]-l[i]+)*1LL*(r[i-]-l[i-]+)%MOD, MOD-)%MOD, p[i] = (-q[i]+MOD)%MOD;
sum[i] = (sum[i-] + p[i]) % MOD;
}
LL ans = sum[n];
for (int i = ; i <= n ;++i) {
LL tot = sum[n];
tot = (tot - p[i])%MOD;
if(i- >= ) tot = (tot - p[i-]) % MOD;
if(i+ <= n) tot = (tot - p[i+]) % MOD;
tot = (tot + MOD) % MOD;
ans = (ans + p[i]*tot%MOD) % MOD;
}
for (int i = ; i < n; ++i) {
LL tot = ((-q[i]-q[i+])%MOD+MOD)%MOD;
if(i- >= ) {
int L = max(l[i-], max(l[i], l[i+])), R= min(r[i-], min(r[i], r[i+]));
if(L <= R)tot = (tot + (R-L+)*q_pow((r[i]-l[i]+)*1LL*(r[i-]-l[i-]+)%MOD*(r[i+]-l[i+]+)%MOD, MOD-)%MOD)%MOD;
}
ans = (ans + *tot) % MOD;
}
printf("%lld\n", ans%MOD);
return ;
}

最新文章

  1. python 2.7 学习笔记--文件的基本操作
  2. Java内存区域与内存溢出异常
  3. R语言多项式回归
  4. SVN提交.a文件的方法
  5. C# 线程(二):关于线程的相关概念
  6. JS中的集中页面跳转的方法
  7. sql 批量操作(存在的更新,不存在的插入)
  8. Jquery中去除左右空格
  9. Linux以KB显示内存大小
  10. mysql判断某个字符串是否包含另一个
  11. tp可用的超强第三方图表类库-JpGraph
  12. trangle
  13. [转]Zend Studio 文件头和方法注释设置
  14. Xamarin android SwipeRefreshLayout入门实例
  15. bzoj 5016: [Snoi2017]一个简单的询问
  16. Anaconda安装及使用
  17. 【原】Java学习笔记020 - 面向对象
  18. ORACLE11g 重装系统后根据dbf恢复数据库
  19. js 分页问题
  20. 2019-03-18-day013-装饰器与内置函数

热门文章

  1. 金仓Kingbase数据库网页数据维护分析工具
  2. 最新 蚂蚁金服java校招面经 (含整理过的面试题大全)
  3. vue中的$listeners属性作用
  4. 【C++札记】动态分配内存(malloc,free)
  5. 标准Trie、压缩Trie、后缀Trie
  6. Mobile Phone Network CodeForces - 1023F (最小生成树)
  7. (一)JNDI基础
  8. (七)Redis之持久化之RDB方式
  9. (二)发布第一个WebService服务与DSWL文档解析
  10. IDEA如何本机调试springboot应用打的jar包