#line 1 "library/Template/template.hpp" #include using namespace std; #define rep(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define ALL(v) (v).begin(),(v).end() #define UNIQUE(v) sort(ALL(v)),(v).erase(unique(ALL(v)),(v).end()) #define MIN(v) *min_element(ALL(v)) #define MAX(v) *max_element(ALL(v)) #define LB(v,x) lower_bound(ALL(v),(x))-(v).begin() #define UB(v,x) upper_bound(ALL(v),(x))-(v).begin() using ll=long long int; using ull=unsigned long long; const int inf = 0x3fffffff; const ll INF = 0x1fffffffffffffff; templateinline bool chmax(T& a,T b){if(ainline bool chmin(T& a,T b){if(a>b){a=b;return 1;}return 0;} templateT ceil(T x,U y){assert(y!=0); if(y<0)x=-x,y=-y; return (x>0?(x+y-1)/y:x/y);} templateT floor(T x,U y){assert(y!=0); if(y<0)x=-x,y=-y; return (x>0?x/y:(x-y+1)/y);} templateint popcnt(T x){return __builtin_popcountll(x);} templateint topbit(T x){return (x==0?-1:63-__builtin_clzll(x));} templateint lowbit(T x){return (x==0?-1:__builtin_ctzll(x));} #line 2 "library/Utility/fastio.hpp" #include class FastIO{ static constexpr int L=1<<16; char rdbuf[L]; int rdLeft=0,rdRight=0; inline void reload(){ int len=rdRight-rdLeft; memmove(rdbuf,rdbuf+rdLeft,len); rdLeft=0,rdRight=len; rdRight+=fread(rdbuf+len,1,L-len,stdin); } inline bool skip(){ for(;;){ while(rdLeft!=rdRight and rdbuf[rdLeft]<=' ')rdLeft++; if(rdLeft==rdRight){ reload(); if(rdLeft==rdRight)return false; } else break; } return true; } template::value,int> =0>inline bool _read(T& x){ if(!skip())return false; if(rdLeft+20>=rdRight)reload(); bool neg=false; if(rdbuf[rdLeft]=='-'){ neg=true; rdLeft++; } x=0; while(rdbuf[rdLeft]>='0' and rdLeft=rdRight)reload(); bool neg=false; if(rdbuf[rdLeft]=='-'){ neg=true; rdLeft++; } x=0; while(rdbuf[rdLeft]>='0' and rdLeft=rdRight)reload(); x=0; while(rdbuf[rdLeft]>='0' and rdLeft::value,int> =0>inline bool _read(T& x){ if(!skip())return false; if(rdLeft+20>=rdRight)reload(); bool neg=false; if(rdbuf[rdLeft]=='-'){ neg=true; rdLeft++; } x=0; while(rdbuf[rdLeft]>='0' and rdbuf[rdLeft]<='9' and rdLeft='0' and rdbuf[rdLeft]<='9' and rdLeft=rdRight)reload(); x=rdbuf[rdLeft++]; return true; } inline bool _read(string& x){ if(!skip())return false; for(;;){ int pos=rdLeft; while(pos' ')pos++; x.append(rdbuf+rdLeft,pos-rdLeft); if(rdLeft==pos)break; rdLeft=pos; if(rdLeft==rdRight)reload(); else break; } return true; } templateinline bool _read(vector& v){ for(auto& x:v){ if(!_read(x))return false; } return true; } char wtbuf[L],tmp[50]; int wtRight=0; inline void flush(){ fwrite(wtbuf,1,wtRight,stdout); wtRight=0; } inline void _write(const char& x){ if(wtRight>L-32)flush(); wtbuf[wtRight++]=x; } inline void _write(const string& x){ for(auto& c:x)_write(c); } template::value,int> =0>inline void _write(T x){ if(wtRight>L-32)flush(); if(x==0){ _write('0'); return; } else if(x<0){ _write('-'); if (__builtin_expect(x == std::numeric_limits::min(), 0)) { switch (sizeof(x)) { case 2: _write("32768"); return; case 4: _write("2147483648"); return; case 8: _write("9223372036854775808"); return; } } x=-x; } int pos=0; while(x!=0){ tmp[pos++]=char((x%10)|48); x/=10; } rep(i,0,pos)wtbuf[wtRight+i]=tmp[pos-1-i]; wtRight+=pos; } inline void _write(__int128_t x){ if(wtRight>L-40)flush(); if(x==0){ _write('0'); return; } else if(x<0){ _write('-'); x=-x; } int pos=0; while(x!=0){ tmp[pos++]=char((x%10)|48); x/=10; } rep(i,0,pos)wtbuf[wtRight+i]=tmp[pos-1-i]; wtRight+=pos; } inline void _write(__uint128_t x){ if(wtRight>L-40)flush(); if(x==0){ _write('0'); return; } int pos=0; while(x!=0){ tmp[pos++]=char((x%10)|48); x/=10; } rep(i,0,pos)wtbuf[wtRight+i]=tmp[pos-1-i]; wtRight+=pos; } templateinline void _write(const vector& v){ rep(i,0,v.size()){ if(i)_write(' '); _write(v[i]); } } public: FastIO(){} ~FastIO(){flush();} inline void read(){} template inline void read(Head& head,Tail&... tail){ assert(_read(head)); read(tail...); } templateinline void write(){if(ln)_write('\n');} template inline void write(const Head& head,const Tail&... tail){ if(space)_write(' '); _write(head); write(tail...); } }; /** * @brief Fast IO */ #line 3 "sol.cpp" #line 2 "library/Geometry/geometry.hpp" using T=double; const T eps=1e-12; using Point=complex; using Poly=vector; #define X real() #define Y imag() inline bool eq(const T& a,const T& b){ return fabs(a-b)=0?0:1));}; if(sub(a)!=sub(b))return sub(a)>(istream& is,Point& p){ T x,y; is>>x>>y; p=Point(x,y); return is; } ostream& operator<<(ostream& os,Point& p){ os<eps)return 1; //ccw if(cross(b,c)<-eps)return -1; //cw if(dot(b,c)<0)return 2; //c,a,b if(norm(b)a.r+b.r+eps)return 4; if(eq(d,a.r+b.r))return 3; if(eq(d,abs(a.r-b.r)))return 1; if(da.r+eps)return res; Point h=Projection(b,a.p); if(eq(d,a.r)){ res.push_back(h); return res; } Point e=unit(b.b-b.a); T ph=sqrt(a.r*a.r-d*d); res.push_back(h-e*ph); res.push_back(h+e*ph); return res; } Poly Intersection(const Circle& a,const Segment& b){ Line c(b.a,b.b); Poly sub=Intersection(a,c); double xmi=min(b.a.X,b.b.X),xma=max(b.a.X,b.b.X); double ymi=min(b.a.Y,b.b.Y),yma=max(b.a.Y,b.b.Y); Poly res; rep(i,0,sub.size()){ if(xmi<=sub[i].X+eps and sub[i].X-eps<=xma and ymi<=sub[i].Y+eps and sub[i].Y-eps<=yma){ res.push_back(sub[i]); } } return res; } Poly Intersection(const Circle& a,const Circle& b){ Poly res; int mode=isIntersect(a,b); T d=abs(a.p-b.p); if(mode==4 or mode==0)return res; if(mode==3){ T t=a.r/(a.r+b.r); res.push_back(a.p+(b.p-a.p)*t); return res; } if(mode==1){ if(b.rc.r-eps)return c.r*c.r*arg(vb*conj(va)); auto u=Intersection(c,Segment(p1,p2)); Poly sub; sub.push_back(p1); for(auto& x:u)sub.push_back(x); sub.push_back(p2); rep(i,0,sub.size()-1)res+=self(self,c,sub[i],sub[i+1]); return res; }; T res=.0; rep(i,0,n)res+=rec(rec,b,a[i],a[(i+1)%n]); return fabs(res/2.); } T Area(const Circle& a,const Circle& b){ T d=abs(a.p-b.p); if(d>=a.r+b.r-eps)return .0; if(d<=abs(a.r-b.r)+eps){ T r=min(a.r,b.r); return M_PI*r*r; } T ath=acos((a.r*a.r+d*d-b.r*b.r)/d/a.r/2.); T res=a.r*a.r*(ath-sin(ath*2)/2.); T bth=acos((b.r*b.r+d*d-a.r*a.r)/d/b.r/2.); res+=b.r*b.r*(bth-sin(bth*2)/2.); return fabs(res); } bool isConvex(const Poly& a){ int n=a.size(); int cur,pre,nxt; rep(i,0,n){ pre=(i-1+n)%n; nxt=(i+1)%n; cur=i; if(ccw(a[pre],a[cur],a[nxt])==-1)return 0; } return 1; } int isContained(const Poly& a,const Point& b){ // 0:not contain,1:on edge,2:contain bool res=0; int n=a.size(); rep(i,0,n){ Point p=a[i]-b,q=a[(i+1)%n]-b; if(p.Y>q.Y)swap(p,q); if(p.Yeps)res^=1; if(eq(cross(p,q),.0) and dot(p,q)=2 and cross(res[k-1]-res[k-2],a[i]-res[k-1])<-eps)k--; } for(int i=n-2,t=k+1;i>=0;res[k++]=a[i--]){ while(k>=t and cross(res[k-1]-res[k-2],a[i]-res[k-1])<-eps)k--; } res.resize(k-1); return res; } T Diam(const Poly& a){ int n=a.size(); int x=0,y=0; rep(i,1,n){ if(a[i].Y>a[x].Y)x=i; if(a[i].YT{ if(rb-lb<=1)return (T)INF; int mid=(lb+rb)>>1; auto x=a[mid].X; T res=min(self(self,lb,mid),self(self,mid,rb)); inplace_merge(a.begin()+lb,a.begin()+mid,a.begin()+rb, [&](auto p,auto q){return p.Y=res)continue; rep(j,0,ptr){ auto sub=a[i]-buf[ptr-1-j]; if(sub.Y>=res)break; chmin(res,abs(sub)); } buf[ptr++]=a[i]; } return res; }; return rec(rec,0,n); } Circle Incircle(const Point& a,const Point& b,const Point& c){ T A=abs(b-c),B=abs(c-a),C=abs(a-b); Point p(A*a.X+B*b.X+C*c.X,A*a.Y+B*b.Y+C*c.Y); p/=(A+B+C); T r=Dist(Line(a,b),p); return Circle(p,r); } Circle Circumcircle(const Point& a,const Point& b,const Point& c){ Line l1((a+b)/2.,(a+b)/2.+(b-a)*Point(0,1)); Line l2((b+c)/2.,(b+c)/2.+(c-b)*Point(0,1)); Point p=Intersection(l1,l2); return Circle(p,abs(p-a)); } Poly tangent(const Point& a,const Circle& b){ return Intersection(b,Circle(a,sqrt(norm(b.p-a)-b.r*b.r))); } vector tangent(const Circle& a,const Circle& b){ vector res; T d=abs(a.p-b.p); if(eq(d,0))return res; Point u=unit(b.p-a.p); Point v=u*Point(0,1); for(int t:{-1,1}){ T h=(a.r+b.r*t)/d; if(eq(h*h,1)){ res.push_back(Line(a.p+(h>0?u:-u)*a.r, a.p+(h>0?u:-u)*a.r+v)); } else if(1>h*h){ Point U=u*h,V=v*sqrt(1-h*h); res.push_back(Line(a.p+(U+V)*a.r,b.p-(U+V)*(b.r*t))); res.push_back(Line(a.p+(U-V)*a.r,b.p-(U-V)*(b.r*t))); } } return res; } /** * @brief Geometry */ #line 5 "sol.cpp" FastIO io; int main(){ int a,b,c; io.read(a,b,c); Point A(0,0),B(b,0); rep(x,1,b+c)rep(y,1,b+c){ double cosa=double(b*b+c*c-(x+y+a)*(x+y+a))*1./(2*b*c); if(cosa>1.-eps or fabs(cosa)