博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
BZOJ3240 NOI2013矩阵游戏(数论)
阅读量:4952 次
发布时间:2019-06-12

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

  必修五题。

// luogu-judger-enable-o2#include
#include
#include
#include
#include
#include
using namespace std;int read(){ int x=0,f=1;char c=getchar(); while (c<'0'||c>'9') { if (c=='-') f=-1;c=getchar();} while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar(); return x*f;}#define N 1000010#define P 1000000007int n,m,a[N],b[N],c,d,u,v,p,q,ans,x,y,k,tmp1,tmp2;int ksm(int a,int k){ if (k==0) return 1; int tmp=ksm(a,k>>1); if (k&1) return 1ll*tmp*tmp%P*a%P; else return 1ll*tmp*tmp%P;}int inv(int x){ return ksm((x+P)%P,P-2);}int main(){ char ch=getchar(); while (ch>='0'&&ch<='9') k=(10ll*k+(ch^48))%(P-1),tmp1=(10ll*tmp1+(ch^48))%P,ch=getchar(); while (ch<'0'||ch>'9') ch=getchar(); while (ch>='0'&&ch<='9') p=(10ll*p+(ch^48))%(P-1),tmp2=(10ll*tmp2+(ch^48))%P,ch=getchar(); u=read(),v=read(),c=read(),d=read(); if (u>1) p=(p+P-2)%(P-1),p=ksm(u,p),q=1ll*v*inv(u-1)%P*(p+P-1)%P; else tmp2=(tmp2+P-1)%P,q=1ll*tmp2*v%P,p=1; x=p,y=q; p=1ll*p*c%P;q=(d+1ll*c*q%P)%P; if (p>1) k=(k+P-2)%(P-1),q=1ll*q*inv(p-1)%P,p=ksm(p,k),q=1ll*q*(p+P-1)%P,ans=(p+q)%P; else tmp1=(tmp1+P-1)%P,ans=(1+1ll*tmp1*q)%P; ans=(1ll*ans*x%P+y+P)%P; cout<

 

转载于:https://www.cnblogs.com/Gloid/p/9662750.html

你可能感兴趣的文章
前端性能优化集【持续更新】
查看>>
第二章:webdriver 控制浏览器窗口大小
查看>>
四则运算2初步构思
查看>>
Break the Chocolate(规律)
查看>>
C#jbox小节
查看>>
结构体指针释放的问题
查看>>
C#枚举Enum[轉]
查看>>
第三百五十七天 how can I 坚持
查看>>
【动态规划】流水作业调度问题与Johnson法则
查看>>
startActivityForResult不起作用
查看>>
Python&Selenium&Unittest&BeautifuReport 自动化测试并生成HTML自动化测试报告
查看>>
活现被翻转生命
查看>>
POJ 1228
查看>>
SwaggerUI+SpringMVC——构建RestFul API的可视化界面
查看>>
springmvc怎么在启动时自己执行一个线程
查看>>
流操作的规律
查看>>
Python基础学习15--异常的分类与处理
查看>>
javascript运算符的优先级
查看>>
React + Redux 入门(一):抛开 React 学 Redux
查看>>
13位时间戳和时间格式化转换,工具类
查看>>