B. Ela’s Fitness and the Luxury Number
思路:
对于一个因子x的都有三个数能整除它,分别为x*x x*x+x x*x+2*x
。然后去讨论它的边界问题。
一个注意事项:不能用sqrt直接开方,要用二分逼近取下边界。
#include <bits/stdc++.h>
#define int long long
#define ios cin.tie(0),cout.tie(0),ios::sync_with_stdio(0);
#define endl '\n'
#define ULL unsigned long long
#define down 0.996
using namespace std;
const double eps=1e-8;
const double inf=1e18;
const int mod=998244353;
const int P=131;
const int N=1e7+5;
int l,r;
int f(int x){
int l =0,r=inf;
while(l<r){
int mid = l+r+1>>1;
if((__int128)mid*mid<=x)l=mid;
else r=mid-1;
}
return r;
}
void solve()
{
cin>>l>>r;
int a=f(l);
int b=f(r);
int ans=(b-a-1)*3;
int a1=a*a,b1=b*b;
int a2=a1+a,a3=a1+2*a;
int b2=b1+b,b3=b1+2*b;
if(l>a2&&l<=a3) ans++;
else if(l>a1&&l<=a2) ans+=2;
else if(l==a1) ans+=3;
if(r==b3) ans+=3;
else if(r<b3&&r>=b2) ans+=2;
else if(r<b2&&r>=b1) ans+=1;
cout<<ans<<endl;
}
signed main()
{
//ios;
int T;cin>>T;
while(T--)
solve();
return 0;
}
G. Shinyruo and KFC
好废物啊。
思路:
1.O(n*m)的思路肯定过不去。明知道不去想想怎么优化,净想了些乱七八糟的东西。
2.若题目中给出数组a[N]的和不超过一定的量级,则说明数组中不同的数不超过O(logN)
3.相同的数归为一组,然后用快速幂优化就好。笨啊
#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define For(i,a,b) for(i=(a);i<=(b);++i)
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
using namespace std;
const int N=2e5+5;
const int inf=1e18;
const int mod=998244353;
int n,m,a[N],num[N],cnt;
map<int,int>mp;
int fac[N],inv[N];
int qpow(int a,int b){
int res=1;
while(b)
{
if(b&1) res=res*a%mod;
a=a*a%mod;
b>>=1;
}
return res;
}
int getinv(int x){return qpow(x,mod-2);}
int C(int a,int b)
{
return (fac[a]*inv[a-b]%mod)*inv[b]%mod;
}
void solve()
{
fac[0]=1;
for(int i=1;i<=50000;i++) fac[i]=fac[i-1]*i%mod;
inv[50000]=getinv(fac[50000]);
for(int i=50000;i>=1;i--) inv[i-1]=inv[i]*i%mod;
cin>>n>>m;
int mx=0;
for(int i=1;i<=n;i++) cin>>a[i],mp[a[i]]++,mx=max(mx,a[i]);
for(auto x:mp)
num[++cnt]=x.first;
for(int i=1;i<=m;i++)
{
if(i<mx)
{
cout<<0<<endl;continue;
}
int ans=1;
for(int j=1;j<=cnt;j++)
ans=ans*qpow(C(i,num[j]),mp[num[j]])%mod;
cout<<ans<<endl;
}
}
signed main()
{
//ios;
//int T;cin>>T;
//while(T--)
solve();
return 0;
}
J. Circular Billiard Table
gcd最大公约数这个点没有想到。
得出公式:n*2*(a/b)=360*k
,一次碰撞构成的圆心角为2*(a/b)
,回到原点可能经过k圈。
进一步化简:n=k*(b*180/a)
,k保证了该式子会是一个整数,那么如何使得n的值最小,
因此可以将式子化简为:180*b/gcd(180*b,a)
#include<bits/stdc++.h>
#define int long long
#define endl '\n'
#define For(i,a,b) for(i=(a);i<=(b);++i)
#define ios (ios::sync_with_stdio(false),cin.tie(0),cout.tie(0))
using namespace std;
const int N=2e6+5;
const int inf=1e18;
const int mod=998244353;
//int fac[N],inv[N];
//int qpow(int a,int b){
// int res=1;
// while(b)
// {
// if(b&1) res=res*a%mod;
// a=a*a%mod;
// b>>=1;
// }
// return res;
//}
//int getinv(int x){return qpow(x,mod-2);}
//int C(int a,int b)
//{
// return (fac[a]*inv[a-b]%mod)*inv[b]%mod;
//}
int a,b;
void solve()
{
cin>>a>>b;
cout<<b*180/__gcd(180*b,a)-1<<endl;
}
signed main()
{
ios;
int T;cin>>T;
while(T--)
solve();
return 0;
}