码迷,mamicode.com
首页 > 其他好文 > 详细

几何向量gcd+暴力枚举——cf552

时间:2019-08-29 17:49:20      阅读:83      评论:0      收藏:0      [点我收藏+]

标签:几何   style   define   ||   lse   暴力   else   mes   gcd   

两两枚举点(xi,yi)(xj,yj)对于第三个点(xk,yk),只要向量ik和ij方向不同即可

只要(xj-xi,yj-yi) 和 (xk-xi,yk-yi)化简后不相同即可,用每个点开个map存下一个gcd的pair即可

#include<bits/stdc++.h>
using namespace std;
#define N 2005
map<pair<int,int>,int>mp[N];
pair<int,int>p[N];
int n;
int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        int x,y;
        scanf("%d%d",&x,&y);
        p[i]=make_pair(x,y);
    }
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)if(i!=j){
            int dx=p[j].first-p[i].first;
            int dy=p[j].second-p[i].second;
            if(dx==0||dy==0){
                if(dx)dx=1;if(dy)dy=1;
                mp[i][make_pair(dx,dy)]++;
            }
            else {
                int d=__gcd(dx,dy);
                mp[i][make_pair(dx/d,dy/d)]++;
            }
        }
    long long ans=0;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)if(i!=j){
            int dx=p[j].first-p[i].first;
            int dy=p[j].second-p[i].second;
            if(dx==0||dy==0){
                if(dx==0){
                    dy=1;
                    ans+=n-2-(mp[i][make_pair(dx,dy)]-1);
                }
                else{
                    dx=1;
                    ans+=n-2-(mp[i][make_pair(dx,dy)]-1);
                }
            }
            else {
                int d=__gcd(dx,dy);
                ans+=n-2-(mp[i][make_pair(dx/d,dy/d)]-1);
            }
        }
    cout<<ans/6<<\n;
} 

 

几何向量gcd+暴力枚举——cf552

标签:几何   style   define   ||   lse   暴力   else   mes   gcd   

原文地址:https://www.cnblogs.com/zsben991126/p/11430979.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!