首页 > HDU 4913 Least common multiple

HDU 4913 Least common multiple

/*
hdu4913 Least common multiple
http://acm.hdu.edu.cn/showproblem.php?pid=4913
离散化 线段树 统计逆序数思想
tips:
1、线段树中一定要到处都取模,否则wa。。。
2、lazy是乘积的形式出现,不是加和*/
#include 
#include 
#include 
using namespace std;
const long long Nmax=100005LL;
const long long mod=1000000007LL;
long long ans;
int n;
vector<long long> v;
struct Node
{long long a;long long b;
};
Node num[Nmax];
long long a[Nmax],b[Nmax],ab[Nmax];
bool cmp(Node a,Node b)
{if(a.a==b.a)return a.b<b.b;return  a.a<b.a;
}
int get_id(long long x)
{return lower_bound(v.begin(),v.end(),x)-v.begin()+1;
}
struct Tree
{int l;int r;long long sum;long long times;long long lazy;
}tree[Nmax*4];
void build(int root,int l,int r)
{tree[root].l=l;tree[root].r=r;tree[root].sum=tree[root].times=0LL;tree[root].lazy=1LL;if(l==r)return;int mid=(l+r)>>1;build(root<<1,l,mid);build(root<<1|1,mid+1,r);
}void push_down(int root)
{if(tree[root].l==tree[root].r){tree[root].lazy=1LL;return;}if(tree[root].lazy!=1LL){tree[root<<1].lazy=(tree[root<<1].lazy*tree[root].lazy)%mod;tree[root<<1|1].lazy=(tree[root<<1|1].lazy*tree[root].lazy)%mod;tree[root<<1].sum=(tree[root<<1].sum*tree[root].lazy)%mod;tree[root<<1|1].sum=(tree[root<<1|1].sum*tree[root].lazy)%mod;tree[root].lazy=1LL;}
}void insert(int root,int l,int r,int data)
{if(tree[root].l==l&&tree[root].r==r){tree[root].times+=data;return;}push_down(root);int mid=(tree[root].l+tree[root].r)>>1;if(mid>=l)insert(root<<1,l,r,data);if(mid<r)insert(root<<1|1,l,r,data);tree[root].times=tree[root<<1].times+tree[root<<1|1].times;
}
void insert_sum(int root,int l,int r,int data)
{if(tree[root].l>=l&&tree[root].r<=r){tree[root].sum=(tree[root].sum+data)%mod;return;}push_down(root);int mid=(tree[root].l+tree[root].r)>>1;if(mid>=l)insert_sum(root<<1,l,r,data);if(mid<r)insert_sum(root<<1|1,l,r,data);tree[root].sum=(tree[root<<1].sum+tree[root<<1|1].sum)%mod;
}long long query_sum(int root,int l,int r)
{if(tree[root].l>=l&&tree[root].r<=r){return tree[root].sum%mod;}push_down(root);int mid=(tree[root].l+tree[root].r)>>1;long long ans=0LL;if(mid>=l)ans=( ans+query_sum(root<<1,l,r) )%mod;if(mid<r)ans=( ans+query_sum(root<<1|1,l,r) )%mod;while(ans<0)ans+=mod;return ans;
}long long query(int root,int l,int r)
{if(tree[root].l>=l&&tree[root].r<=r){return tree[root].times;}int mid=(tree[root].l+tree[root].r)>>1;long long ans=0LL;if(mid>=l)ans+=query(root<<1,l,r);if(mid<r)ans+=query(root<<1|1,l,r);return ans;
}void init()
{v.clear();ans=0LL;build(1,1,n);
}int qpow(long long base,long long n)
{base=base%mod;long long ans=1LL;while(n>0){if(n&1)ans=(ans*base)%mod;base=base*base%mod;n>>=1;}while(ans<0)ans+=mod;return ans;
}void update(int root,int l,int r)
{if(tree[root].l>=l && tree[root].r<=r){tree[root].lazy=(tree[root].lazy*2LL)%mod;tree[root].sum=(tree[root].sum*2LL)%mod;return;}push_down(root);int mid=(tree[root].l+tree[root].r)>>1;if(mid>=l)update(root<<1,l,r);if(mid<r)update(root<<1|1,l,r);tree[root].sum=(tree[root<<1].sum+tree[root<<1|1].sum)%mod;
}void watch(int root,int l,int r)
{printf("tree[%d]: l:%d,r:%d,sum:%lld,times:%lld,lazy:%lld
",root,tree[root].l,tree[root].r,tree[root].sum,tree[root].times,tree[root].lazy);if(l==r)return;int mid=(l+r)>>1;watch(root<<1,1,mid);watch(root<<1|1,mid+1,r);
}
int main()
{//freopen("hdu4913.in","r",stdin);while(scanf("%d",&n)==1){init();for(int i=1;i<=n;i++){scanf("%lld%lld",&num[i].a,&num[i].b);v.push_back(num[i].b);}sort(v.begin(),v.end());v.erase(unique(v.begin(),v.end()),v.end());sort(num+1,num+1+n,cmp);int m=v.size();for(int i=1;i<=n;i++){//printf("times[%d]:
",i);int b=get_id(num[i].b);//printf("b:%d
",b);//printf("before:
");//watch(1,1,n);//printf("after:
");//watch(1,1,n);//printf("end!!!!!!!!!!!!!!!!!
");insert(1,b,b,1);long long x=query(1,1,b);//printf("%lld:%lld
",query(1,1,b),query(1,b,b)); //printf("x:%lld
",x); x=qpow(2,x-1);//printf("x:%lld
",x);
x= (x*qpow(3,num[i].b)) %mod;//printf("x:%lld
",x);ans=( ans+ qpow(2,num[i].a)*x )%mod;//printf("ans:%lld
",ans);//if(kkk==0)//if(b!=n)if(b!=n)ans=(ans+(qpow(2,num[i].a)* query_sum(1,b+1,n)) %mod)%mod;//printf("ans:%lld
",ans);if(b!=n)update(1,b+1,n);insert_sum(1,b,b,x);//watch(1,1,n);
        }while(ans<0)ans+=mod;//printf("$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$


");printf("%lld
",ans);}return 0;
}

 

转载于:https://www.cnblogs.com/BBBob/p/6522847.html

更多相关:

  • 说明如下: (1)所有操作最好使用root操作,以尽可能避免权限问题 (2)crtmpserver和web服务器apache(Ngnix亦可)被装在同一台服务器上,客户端使用其他PC的桌面浏览器和android手机,下面是测试环境示意图: (3)测试环境位于局域网内,由一台路由器接入联通运营商提供的10M包年网络服务 (...

  • 1. free 2. top 3. vmstat 4. slabtop; 5. pmap 6. dmesg 7. /proc/meminfo 8. /proc/sys/vm 目录下的文件 9. sync 10./proc/zoneinfo  11./proc/pagetypeinfo 查看内存工具:1.free free - Dis...

  • 文章目录前言创建二叉树先序遍历中序遍历后序遍历获取叶子节点个数获取树的高度测试代码 前言 现有如下二叉树: 关于二叉树的相关操作,我们能够发现二叉树从根节点到子节点,以及每个中间节点基本都能够拆分为若干个子节点的操作,且每个子节点的操作都和其头节点操作一致。 所以我们针对二叉树的操作都可以使用分治算+回溯/归并算法进行...

  • 今天需要部署一个ceph L 版本12.2.12的环境,无奈最近公司网络无法访问到ceph官网,只能使用之前下载好的ceph-deploy-1.5.39版本,安装上之后一口老血喷了出来,没有mgr的部署选项。 无奈之下只能自己制作一个1.5.38版本的ceph-deploy包,借用从ceph-deploy-1.5.39-0.src....

  • root 权限进入MySQL: mysql –uroot 查看当前MySQL用户: select user,host from mysql.user;     此处以User为root,Host为localhost的账户为例,将密码改为password的命令:   SET PASSWORD FOR 'root'@'localhost...

  • L3-010. 是否完全二叉搜索树 时间限制 400 ms 内存限制 65536 kB 代码长度限制 8000 B 判题程序 Standard 作者 陈越 将一系列给定数字顺序插入一个初始为空的二叉搜索树(定义为左子树键值大,右子树键值小),你需要判断最后的树是否一棵完全二叉树,并且给出其层序遍历的结果。 输入格...

  •  一个用于实现初始化指定个数的完全二叉树,以及两个非递归的深度优先遍历,和广度优先遍历 package fifth;  import java.util.Random;  public class Tool{     public static Random rand= new Random(); }  --------------...

  •         Apache POI是一个开源的利用Java读写Excel,WORD等微软OLE2组件文档的项目。        我的需求是对Excel的数据进行导入或将数据以Excel的形式导出。先上简单的测试代码:package com.xing.studyTest.poi;import java.io.FileInputSt...

  • 要取得[a,b)的随机整数,使用(rand() % (b-a))+ a; 要取得[a,b]的随机整数,使用(rand() % (b-a+1))+ a; 要取得(a,b]的随机整数,使用(rand() % (b-a))+ a + 1; 通用公式:a + rand() % n;其中的a是起始值,n是整数的范围。 要取得a到b之间的...

  • 利用本征图像分解(Intrinsic Image Decomposition)算法,将图像分解为shading(illumination) image 和 reflectance(albedo) image,计算图像的reflectance image。 Reflectance Image 是指在变化的光照条件下能够维持不变的图像部分...

  • 题目:面试题39. 数组中出现次数超过一半的数字 数组中有一个数字出现的次数超过数组长度的一半,请找出这个数字。 你可以假设数组是非空的,并且给定的数组总是存在多数元素。 示例 1: 输入: [1, 2, 3, 2, 2, 2, 5, 4, 2] 输出: 2 限制: 1 <= 数组长度 <= 50000 解题: cl...

  • 题目:二叉搜索树的后序遍历序列 输入一个整数数组,判断该数组是不是某二叉搜索树的后序遍历结果。如果是则返回 true,否则返回 false。假设输入的数组的任意两个数字都互不相同。 参考以下这颗二叉搜索树:      5     /    2   6   /  1   3示例 1: 输入: [1,6,3,2,5] 输出...