杭电多校第七场-J-Sequence
题目描述

Your job is simple, for each task, you should output Fn module 109+7.
输入
Then, for the next T lines, each line consists of 6 integers, A , B, C, D, P, n.
1≤T≤20 0≤A,B,C,D≤109 1≤P,n≤109
输出
样例输入
2
3 3 2 1 3 5
3 2 2 2 1 4
样例输出
36
24
矩阵快速幂
问题在于p/i
上图为p=100的表,可以发现p/i相同的连续的一段的最后一个数的坐标为p/(p/i) ,然后分段做就可以了
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int p=1e9+;
int n,P,T;
ll a,b,c,d;
struct Mat
{
ll v[][]; Mat() {
memset(v, , sizeof(v));
}
void init()
{
for (int i=;i<;i++)
v[i][i]=(ll);
} };
Mat operator *(Mat a,Mat b)
{
Mat c;
for (int i=;i<;i++)
{
for (int j=;j<;j++)
{
c.v[i][j]=;
for (int k=;k<;k++)
c.v[i][j]+=((a.v[i][k]%p*b.v[k][j]%p)%p+p)%p;
}
}
return c;
}
Mat qmod(Mat a,int k)
{
Mat c;
c.init(); while (k)
{
if (k&) c=c*a;
a=a*a;
k>>=;
}
return c;
}
void solve()
{
scanf("%lld%lld%lld%lld%d%d",&a,&b,&c,&d,&P,&n);
Mat mp;
mp.v[][]=d; mp.v[][]=c;
mp.v[][]=; mp.v[][]=; for (int i=;i<=n;)
{ if(P/i==)
{
Mat w=mp;
w=qmod(w,n-i+);
cout<<(w.v[][]*b+w.v[][]*a+w.v[][])%p<<endl;
return;
} int j=min(n,P/(P/i));
Mat w=mp;
w.v[][]=P/i;
w=qmod(w,j-i+);
ll ans_a=(w.v[][]*b+w.v[][]*a+w.v[][])%p;
ll ans_b=(w.v[][]*b+w.v[][]*a+w.v[][])%p;
a=ans_a;b=ans_b;
i=j+; }
printf("%lld\n",b);
}
int main()
{
scanf("%d",&T);
while (T--) solve();
return ;
}
杭电多校第七场-J-Sequence的更多相关文章
- 杭电多校第七场 1010 Sequence(除法分块+矩阵快速幂)
Sequence Problem Description Let us define a sequence as below f1=A f2=B fn=C*fn-2+D*fn-1+[p/n] Your ...
- 2019杭电多校第七场 HDU - 6656 Kejin Player——概率&&期望
题意 总共有 $n$ 层楼,在第 $i$ 层花费 $a_i$ 的代价,有 $pi$ 的概率到 $i+1$ 层,否则到 $x_i$($x_i \leq 1$) 层.接下来有 $q$ 次询问,每次询问 $ ...
- 【杭电多校第七场】A + B = C
原题: Given a,b,c, find an arbitrary set of x,y,z such that a*10^x+b*10^y=c*10^z and 0≤x,y,z≤10^6. 给你三 ...
- [2019杭电多校第七场][hdu6646]A + B = C(hash)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6646 题意为求a*10x+b*10y=c*10z满足公式的任意一组解x,y,z. 因为c有可能会由a+ ...
- hdu61272017杭电多校第七场1008Hard challenge
Hard challenge Time Limit: 4000/2000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) ...
- 2017杭电多校第七场1011Kolakoski
Kolakoski Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) Tota ...
- 2017杭电多校第七场1005Euler theorem
Euler theorem Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 524288/524288 K (Java/Others) ...
- [2019杭电多校第七场][hdu6656]Kejin Player
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6656 题意为从i级花费a元有p的概率升到i+1级,有1-p的概率降到x级(x<i),查询从L级升 ...
- [2019杭电多校第七场][hdu6655]Just Repeat
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6655 题意是说两个人都有一些带有颜色的牌,两人轮流出牌,但是不能出对面出过的颜色的牌,最后谁不能出牌谁 ...
随机推荐
- Centos7 下nginx nginx-1.13.4 安装
环境:CentOS Linux release 7.3.1611 (Core) Linux localhost.localdomain 3.10.0-514.26.2.el7.x86_64 #1 S ...
- Android 开发错误集锦
1. eclipse的Device中不显示手机 在eclipse中连接不上手机,出现adb server didn't ACK fail to start daemon 错误. 出现这种原因是因为a ...
- python学习第二天-基本数据类型常用方法
1.直入主题 python中基本的数据类型有 数字(整形,长整形,浮点型,复数) 字符串 字节串:在介绍字符编码时介绍字节bytes类型 列表 元组 字典 集合 下面我们直接将以下面几个点进行学习 # ...
- 转 Js 跨域CORS报错 Response for preflight has invalid HTTP status code 405
转自:http://www.cnblogs.com/SilenceTom/p/6697484.html 调用接口遇到Response for preflight has invalid HTTP st ...
- nginx的平滑升级,不间断服务
nginx的平滑升级,不间断服务 Nginx更新真的很快,最近nginx的1.0.5稳定版,nginx的0.8.55和nginx的0.7.69旧的稳定版本已经发布.我一项比较喜欢使用新版本的软件, ...
- django使用ajax提交表单数据报403错解决方法
只需要在.ajaxSetup方法中设置csrfmiddlewaretoken即可 $.ajaxSetup({ data: {csrfmiddlewaretoken: '{{ csrf_token }} ...
- 异步执行任务SimpleAsyncTaskExecutor详解
SimpleAsyncTaskExecutor 异步执行用户任务的SimpleAsyncTaskExecutor.每次执行客户提交给它的任务时,它会启动新的线程,并允许开发者控制并发线程的上限(con ...
- 第52天:offset家族、scroll家族和client家族的区别
一.offset家族 1.offsetWidth offsetHeight offsetLeft offsetTop offsetParent共同组成了offset家族,用来获取元素尺寸. offse ...
- ASP.NET MVC下使用文件上传和IIS7下的默认设置限制了上传大小的方法
不多说了,直接用别人的 http://www.cnblogs.com/jiekzou/p/4491505.html
- RT-thread内核之线程调度算法
一个操作系统如果只是具备了高优先级任务能够“立即”获得处理器并得到执行的特点,那么它仍然不算是实时操作系统.因为这个查找最高优先级线程的过程决定了调度时间是否具有确定性,例如一个包含n个就绪任务的系统 ...