BZOJ5131: [CodePlus2017年12月]可做题2

BZOJ没有题面,差评

洛谷的题目链接

题解

其实这题很久之前就写了,也想写个题解但是太懒了,咕到了今天

BZOJ5131: [CodePlus2017年12月]可做题2

在typora写完题解不想copy过来再改格式了,于是直接贴截图qwq

#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 int a1 , l , r , k , p , m ;
struct matrix {
int m[][] ;
matrix() { memset(m,,sizeof(m)); }
matrix operator * (const matrix &x) {
matrix ans;
for(int i = ; i < ; i ++) {
for(int j = ; j < ; j ++) {
for(int k = ; k < ; k ++) {
ans.m[i][j] = (ans.m[i][j] + m[i][k] * x.m[k][j]) % p ;
}
}
}
return ans ;
}
} base , ans ; void power(int b) {
base.m[][] = base.m[][] = base.m[][] = ; base.m[][] = ;
ans.m[][] = ans.m[][] = ; ans.m[][] = ans.m[][] = ;
while(b) {
if(b&) ans = ans * base ;
base = base * base ;
b >>= ;
}
} int x , y ;
int exgcd(int a , int b) {
if(b == ) { x = ; y = ; return a ; }
int Ans = exgcd(b , a % b) , t = x ;
x = y ; y = t - (a / b) * y ;
return Ans ;
} int find(int x , int t) {
int l = , r = t / p + ;
while(l <= r) {
int mid = (l + r) >> ;
if(x + p * mid >= t) r = mid - ;
else l = mid + ;
}
return l ;
} signed main() {
int T = read() ;
while(T--) {
a1 = read() , l = read() , r = read() , k = read() , p = read() , m = read() ;
a1 %= p ; power(k - ) ;
int mod = (m - a1 * ans.m[][] % p + p) % p ;
int gcd = exgcd(ans.m[][] , p) ;
if(mod % gcd != ) { puts("") ; continue ; }
x = x * (mod/gcd) ; p /= gcd ; x = (x % p + p) % p ;
outn( find(x , r+) - find(x , l) ) ;
}
}
上一篇:c++ string c_str() 和data()区别


下一篇:Netty 聊天小程序