题意
求满足如下不等式的非负整数 (x,y) 的对数
[ax+byleq c
]
Sol
a,b,c 都是非负的,那么先随便变个形:
[yleqfrac{c-ax}{b}
]
显然答案就是:
[sum_{x=0}^{lfloorfrac{c}{a}
floor}igg(lfloor frac{c-ax}{b}
floor+1igg)
]
类欧板子题了,但是这个斜率居然是一个负数 ,不能直接套用以前的做法。
一开始想着先把这条直线往 (x) 轴对称然后向上平移,但是这样太麻烦了。
令 (n=lfloor frac{c}{a}
floor) ,直接把直线沿着 (n/2) 对称就好了,这样在下方的点不受影响,左上方的点也被对应到了右上方,并且整点依然是整点。
code:
#include<bits/stdc++.h>
using namespace std;
template<class T>inline void init(T&x){
x=0;char ch=getchar();bool t=0;
for(;ch>'9'||ch<'0';ch=getchar()) if(ch=='-') t=1;
for(;ch>='0'&&ch<='9';ch=getchar()) x=(x<<1)+(x<<3)+(ch-48);
if(t) x=-x;return;
}typedef long long ll;
ll a,b,c;
inline ll gcd(ll a,ll b){return b? gcd(b,a%b):a;}
inline ll likegcd(ll n,ll a,ll b,ll c){
if(!n) return b/c;if(!a) return b/c*(n+1);
ll d=abs(gcd(gcd(a,b),c));
if(d>1) a/=d,b/=d,c/=d;
if(a>=c||b>=c) return b/c*(n+1)+a/c*(n+1)*n/2+likegcd(n,a%c,b%c,c);
ll m=(a*n+b)/c;
return n*m-likegcd(m-1,c,c-b-1,a);
}
int main()
{
init(a),init(b),init(c);
ll n=c/a;ll ans=n+1+likegcd(n,a,c-a*n,b);
cout<<ans<<endl;
return 0;
}