題目鏈接：點(diǎn)擊查看

題目大意：給出二維平面上的 $n$ 個(gè)點(diǎn)，任意選擇出三個(gè)點(diǎn)可以構(gòu)成一個(gè)三角形，現(xiàn)在問(wèn)滿足下面條件的三角形的個(gè)數(shù)：

三角形面積為整數(shù)

三角形包含的（不包括邊界）整數(shù)點(diǎn)為奇數(shù)

其中所有坐標(biāo)點(diǎn)的 $x$ 和 $y$ 都是偶數(shù)

題目分析：因?yàn)樗凶鴺?biāo)點(diǎn)都是偶數(shù)，所以先猜一波任意三個(gè)點(diǎn)都是滿足第一個(gè)條件的，也就是圍成的三角形面積一定是偶數(shù)。

所以現(xiàn)在問(wèn)題就轉(zhuǎn)換為了，如何求解滿足第二個(gè)條件的三角形的個(gè)數(shù)。賽中經(jīng)隊(duì)友提醒，得知了“皮克定理”這一關(guān)鍵的結(jié)論，即 $S=a+\frac{b}{2}?1$，其中 $S$ 為三角形的面積，$a$ 為三角形內(nèi)部的點(diǎn)數(shù)，$b$ 為三角形邊界上的點(diǎn)數(shù)。$b$ 非常好求，每條邊上的答案就是 $\gcd (dx,dy)+1$，放在三角形上就是 $\gcd (|x1?x2|,|y1?y2|)+\gcd (|x1?x3|,|y1?y3|)+\gcd (|x3?x2|,|y3?y2|)+3?3$，因?yàn)槿齻€(gè)頂點(diǎn)分別計(jì)算了兩次，所以減去三之后就是這個(gè)答案了

皮克定理移項(xiàng)之后得到 $a=\frac{b}{2}?1?S$，如果 $a$ 想要為奇數(shù)，那么 $a+1=\frac{b}{2}?S$ 需要是偶數(shù)，因?yàn)樯厦嬉呀?jīng)保證了 $S$ 是偶數(shù)，所以現(xiàn)在只需要保證 $\frac{b}{2}$ 是偶數(shù)

所以問(wèn)題轉(zhuǎn)換為了，有 $n$ 個(gè)點(diǎn)，任意兩個(gè)點(diǎn)之間的邊權(quán)為 $gcd(dx,dy)/2$，現(xiàn)在要求滿足條件的三元環(huán)個(gè)數(shù)，需要滿足三條邊權(quán)為：

奇數(shù)、奇數(shù)、偶數(shù)

偶數(shù)、偶數(shù)、偶數(shù)

然后到這里就尬住了，比賽時(shí)以為是什么典中典的三元環(huán)計(jì)數(shù)問(wèn)題，就下班了

第二天看了題解發(fā)現(xiàn)，邊權(quán)和坐標(biāo)有著密切的聯(lián)系，又因?yàn)槊總€(gè)點(diǎn)的坐標(biāo)都是偶數(shù)，所以我們不妨一開(kāi)始就將坐標(biāo)除以二

然后我們就很神奇的發(fā)現(xiàn)，一共就只有四類點(diǎn)了，分別是 $x$ 和 $y$ 分別取 $0$ 或 $1$ 的時(shí)候。這是點(diǎn)權(quán)，邊權(quán)的話當(dāng)且僅當(dāng) $x1==x2$ and $y1==y2$ 時(shí)，原先點(diǎn)的 $\gcd$ 才是偶數(shù)，其他情況都是奇數(shù)

從頭開(kāi)始再推一遍，發(fā)現(xiàn)我們之前的得到的結(jié)論沒(méi)有變！所以 $n^3$ 去枚舉計(jì)數(shù)就好了

你說(shuō) $n=6000$？不不，現(xiàn)在的 $n$ 實(shí)質(zhì)上只有 $4$ 種情況，用組合數(shù)學(xué)統(tǒng)計(jì)一下答案就好了

代碼：

// Problem: F1. Gregor and the Odd Cows (Easy) // Contest: Codeforces - Codeforces Round #736 (Div. 2) // URL: https://codeforces.com/contest/1549/problem/F1 // Memory Limit: 256 MB // Time Limit: 6000 ms // // Powered by CP Editor (https://cpeditor.org)// #pragma GCC optimize(2) // #pragma GCC optimize("Ofast","inline","-ffast-math") // #pragma GCC target("avx,sse2,sse3,sse4,mmx") #include<iostream> #include<cstdio> #include<string> #include<ctime> #include<cmath> #include<cstring> #include<algorithm> #include<stack> #include<climits> #include<queue> #include<map> #include<set> #include<sstream> #include<cassert> #include<bitset> #include<list> #include<unordered_map> #define lowbit(x) (x&-x) using namespace std; typedef long long LL; typedef unsigned long long ull; template<typename T> inline void read(T &x) {T f=1;x=0;char ch=getchar();while(0==isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(0!=isdigit(ch)) x=(x<<1)+(x<<3)+ch-'0',ch=getchar();x*=f; } template<typename T> inline void write(T x) {if(x<0){x=~(x-1);putchar('-');}if(x>9)write(x/10);putchar(x%10+'0'); } const int inf=0x3f3f3f3f; const int N=1e6+100; LL cnt[2][2]; LL C(int n,int m) {if(m==1) {return n;} else if(m==2) {return 1LL*n*(n-1)/2;} else if(m==3) {return 1LL*n*(n-1)*(n-2)/6;} } int cal(int x1,int y1,int x2,int y2) {if(x1==x2&&y1==y2) {return true;} else {return false;} } int main() { #ifndef ONLINE_JUDGE // freopen("data.in.txt","r",stdin); // freopen("data.out.txt","w",stdout); #endif // ios::sync_with_stdio(false);int n;read(n);for(int i=1;i<=n;i++) {int x,y;read(x),read(y);cnt[x%4/2][y%4/2]++;}LL ans=0;for(int s1=0;s1<4;s1++) {for(int s2=s1;s2<4;s2++) {for(int s3=s2;s3<4;s3++) {int x1=s1/2,y1=s1%2;int x2=s2/2,y2=s2%2;int x3=s3/2,y3=s3%2;int even=0;//2odds + 1even or 3 eveneven+=cal(x1,y1,x2,y2);even+=cal(x1,y1,x3,y3);even+=cal(x2,y2,x3,y3);if(even==2||even==0) {continue;}if(s1==s2&&s2==s3) {ans+=C(cnt[x1][y1],3);} else if(s1==s2) {ans+=C(cnt[x1][y1],2)*cnt[x3][y3];} else if(s2==s3) {ans+=C(cnt[x2][y2],2)*cnt[x1][y1];} else if(s1==s3) {ans+=C(cnt[x1][y1],2)*cnt[x2][y2];} else {ans+=cnt[x1][y1]*cnt[x2][y2]*cnt[x3][y3];}}}}cout<<ans<<endl;return 0; }

總結(jié)

以上是生活随笔為你收集整理的CodeForces - 1549F1 Gregor and the Odd Cows (Easy)(几何+数论)的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。