Premier Bob的算法模板库（II）

最新推荐文章于 2023-09-14 11:14:47 发布
最新推荐文章于 2023-09-14 11:14:47 发布
阅读量616 收藏
点赞数
分类专栏： 数据结构 算法导论 文章标签： 算法 namespace 头文件 C++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/ggn_2015/article/details/70051520
版权

返回上一章 Premier Bob的算法模板库（I）

Geometry_Calculate_Basic .hpp （计算几何基础运算）

#include<cmath>
#include<iostream>
#include<cstdlib>
#include<cstdio>
using namespace std;

namespace Geometry_Calculate_Basic
{
    const double PI=acos(-1.0);

    struct Point{double x,y;Point(double X=0,double Y=0){x=X;y=Y;}};
    typedef Point Vector;

    Vector operator+(Vector A,Vector B){return Vector(A.x+B.x,A.y+B.y);}
    Vector operator-(Point A,Point B){return Point(A.x-B.x,A.y-B.y);}
    Vector operator*(Vector A,double p){return Vector(A.x*p,A.y*p);}
    Vector operator/(Vector A,double p){return Vector(A.x/p,A.y/p);}

    bool operator<(const Point& a, const Point& b){return a.x<b.x || (a.x==b.x && a.y<b.y);}

    const double eps=1e-10;
    int dcmp(double x){if(fabs(x)<eps)return 0;else return x<0 ? -1 : 1;}
    bool operator==(const Point& a,const Point& b){return dcmp(a.x-b.x)==0 && dcmp(a.y-b.y)==0;}

    double Dot(Vector A,Vector B){return A.x*B.x+A.y*B.y;}
    double Length(Vector A){return sqrt(Dot(A,A));}
    double Angle(Vector A,Vector B){return acos(Dot(A,B)/Length(A)/Length(B));}
    double Cross(Vector A,Vector B){return A.x*B.y-A.y*B.x;}

    double Area2(Point A,Point B,Point C){return Cross(B-A,C-A);}

    typedef double ArcAngle;
    Vector Rotate(Vector A,ArcAngle rad){return Vector(A.x*cos(rad)-A.y*sin(rad),A.x*sin(rad)+A.y*cos(rad));}
    Vector Normal(Vector A){double L=Length(A);return Vector(-A.y/L,A.x/L);}

    Point GetLineIntersection(Point P,Vector v,Point Q,Vector w)
    {
        Vector u=P-Q;
        double t=Cross(w,u)/Cross(v,w);
        return P+v*t;
    }

    double DistanceToLine(Point P,Point A,Point B)
    {
        Vector v1=B-A,v2=P-A;
        return fabs(Cross(v1,v2))/Length(v1);
    }

    double DistanceToSegment(Point P,Point A,Point B)
    {
        if(A==B)
            return Length(P-A);
        Vector v1=B-A,v2=P-A,v3=P-B;
        if(dcmp(Dot(v1,v2))<0)return Length(v2);
        else if(dcmp(Dot(v1,v3))>0)
            return Length(v3);
        else return fabs(Cross(v1,v2))/Length(v1);
    }

    Point GetLineProjection(Point P,Point A,Point B)
    {
        Vector v=B-A;
        return A+v*(Dot(v,P-A)/Dot(v,v));
    }

    bool SegmentProperIntersection(Point a1,Point a2,Point b1,Point b2)
    {
        double c1=Cross(a2-a1,b1-a1),c2=Cross(a2-a1,b2-a1),
            c3=Cross(b2-b1,a1-b1),c4=Cross(b2-b1,a2-b1);
        return dcmp(c1)*dcmp(c2)<0 && dcmp(c3)*dcmp(c4)<0;
    }

    double PolygonArea(Point* p,int n)
    {
        double area=0;
        for(int i=1;i<n-1;i++)
            area+=Cross(p[i]-p[0],p[i+1]-p[0]);
        return area/2;
    }
}

但愿这个不是使用样例：

#include"Geometry_Calculate_Basic.hpp"
int main()
{
    int T;cin>>T;
    for(int i=1;i<=T;i++)
    {
        Geometry_Calculate_Basic::Point A,B,C;
        cin>>A.x>>A.y>>B.x>>B.y>>C.x>>C.y;
        Geometry_Calculate_Basic::ArcAngle c_a=Geometry_Calculate_Basic::Angle(A-B,C-B);
        Geometry_Calculate_Basic::ArcAngle a_b=Geometry_Calculate_Basic::Angle(B-C,A-C);
        Geometry_Calculate_Basic::ArcAngle b_c=Geometry_Calculate_Basic::Angle(C-A,A-B);
        Geometry_Calculate_Basic::Vector a1,a2,b1,b2,c1,c2;
        a1=Geometry_Calculate_Basic::Rotate(C-A,b_c/3*2);
        a2=Geometry_Calculate_Basic::Rotate(C-A,b_c/3);

        b1=Geometry_Calculate_Basic::Rotate(A-B,c_a/3*2);
        b2=Geometry_Calculate_Basic::Rotate(A-B,c_a/3);

        c1=Geometry_Calculate_Basic::Rotate(B-C,a_b/3*2);
        c2=Geometry_Calculate_Basic::Rotate(B-C,a_b/3);

        Geometry_Calculate_Basic::Point D,E,F;

        D=Geometry_Calculate_Basic::GetLineIntersection(B,b1,C,c2);
        E=Geometry_Calculate_Basic::GetLineIntersection(A,a1,B,b2);
        F=Geometry_Calculate_Basic::GetLineIntersection(C,c1,A,a2);

        cout<<D.x<<" "<<D.y<<" "<<E.x<<" "<<E.y<<" "<<F.x<<" "<<F.y<<endl;
    }
    return 0;
}

所以说应该写成这样：

#include"Geometry_Calculate_Basic.hpp"
using namespace Geometry_Calculate_Basic;
//using namespace!
int main()
{
    int T;cin>>T;
    for(int i=1;i<=T;i++)
    {
        Point A,B,C;
        cin>>A.x>>A.y>>B.x>>B.y>>C.x>>C.y;
        ArcAngle c_a=Angle(A-B,C-B);
        ArcAngle a_b=Angle(B-C,A-C);
        ArcAngle b_c=Angle(C-A,A-B);
        Vector a1,a2,b1,b2,c1,c2;
        a1=Rotate(C-A,b_c/3*2);
        a2=Rotate(C-A,b_c/3);

        b1=Rotate(A-B,c_a/3*2);
        b2=Rotate(A-B,c_a/3);

        c1=Rotate(B-C,a_b/3*2);
        c2=Rotate(B-C,a_b/3);

        Point D,E,F;

        D=GetLineIntersection(B,b1,C,c2);
        E=GetLineIntersection(A,a1,B,b2);
        F=GetLineIntersection(C,c1,A,a2);

        cout<<D.x<<" "<<D.y<<" "<<E.x<<" "<<E.y<<" "<<F.x<<" "<<F.y<<endl;
    }
    return 0;
}

NlogN_LIS .hpp （nlogn时间复杂度求解LIS长度）

#include<algorithm>
using namespace std;

namespace NlogN_LIS
{
    const int INF=2147483647;
    const int maxn=1048576;
    int g[maxn],d[maxn];
    int LIS(int* Array,int n)
    {
        for(int i=1;i<=n;i++)
            g[i]=INF;
        int ans=0;
        for(int i=1;i<=n;i++)
        {
            int k=lower_bound(g+1,g+n+1,Array[i])-g;
            d[i]=k;
            g[k]=Array[i];
            ans=max(ans,d[i]);
        }
        return ans;
    }
}

UnionFindSet .hpp （并查集模板）

//made by GGN -from NECY
#pragma once

#include<cstdio>
#include<cstdlib>

namespace UnionFindSet
{
    const int UFS_Size=1048576;
    //you can use this to change the size of UFS
    class UFS
    {
        int Father[UFS_Size];
        bool NeedInit;
        int UseSize;

        public:
            UFS(){NeedInit=1;}

            void Init(int NodeCount)
            {
                NeedInit=0;
                for(int i=1;i<=NodeCount;i++)
                    Father[i]=i;
                UseSize=NodeCount;
            }

            int Find(int NodeNum)
            {
                if(NeedInit)
                {
                    printf("Error \"UnionFindSet\":Find  you have to init before use!\n");
                    system("pause>nul");
                    return -1;
                }
                if(NodeNum<=0 || NodeNum>UseSize)
                {
                    printf("Error \"UnionFindSet\":Find NodeNum=%d (<=0 or >UseSize=%d) is not allowed!\n",NodeNum,UseSize);
                    system("pause>nul");
                    return -1;
                }
                if(NodeNum!=Father[NodeNum])
                    Father[NodeNum]=Find(Father[NodeNum]);
                return Father[NodeNum];
            }

            int Union(int NodeX,int NodeY)
            {
                if(NeedInit)
                {
                    printf("Error \"UnionFindSet\":Union  you have to init before use!\n");
                    system("pause>nul");
                    return -1;
                }
                if(NodeX<=0 || NodeX>UseSize)
                {
                    printf("Error \"UnionFindSet\":Union NodeX=%d (<=0 or >UseSize=%d) is not allowed!\n",NodeX,UseSize);
                    system("pause>nul");
                    return -1;
                }
                if(NodeY<=0 || NodeY>UseSize)
                {
                    printf("Error \"UnionFindSet\":Union NodeY=%d (<=0 or >UseSize=%d) is not allowed!\n",NodeY,UseSize);
                    system("pause>nul");
                    return -1;
                }
                Father[Find(NodeY)]=Find(NodeX);
            }

            int Judge(int NodeX,int NodeY)
            {
                if(NeedInit)
                {
                    printf("Error \"UnionFindSet\":Judge  you have to init before use!\n");
                    system("pause>nul");
                    return -1;
                }
                if(NodeX<=0 || NodeX>UseSize)
                {
                    printf("Error \"UnionFindSet\":Judge NodeX=%d (<=0 or >UseSize=%d) is not allowed!\n",NodeX,UseSize);
                    system("pause>nul");
                    return -1;
                }
                if(NodeY<=0 || NodeY>UseSize)
                {
                    printf("Error \"UnionFindSet\":Judge NodeY=%d (<=0 or >UseSize=%d) is not allowed!\n",NodeY,UseSize);
                    system("pause>nul");
                    return -1;
                }
                return Find(NodeX)==Find(NodeY);
            }

            int Ancestor(int Node)
            {
                if(NeedInit)
                {
                    printf("Error \"UnionFindSet\":Ancestor  you have to init before use!\n");
                    system("pause>nul");
                    return -1;
                }
                if(Node<=0 || Node>UseSize)
                {
                    printf("Error \"UnionFindSet\":Ancestor Node=%d (<=0 or >UseSize=%d) is not allowed!\n",Node,UseSize);
                    system("pause>nul");
                    return -1;
                }
                return Father[Node];
            }
    };
}

TernarySearch .hpp （三分法求单峰函数最小值）

//made by GGN -from NEYC
#pragma once

#include<cstdio>
#include<cstdlib>
#include<cmath>
namespace TernarySearch
{
    double TSmin(double L,double R,double (*function)(double x),double eps=1e-7)
    {
        //make a ternary search for min value and return x
        if(R==L)
            return L;
        if(R<L)
        {
            printf("Error \"Ternary Search\":TSmin R=%d (<L=%d) is not allowed!\n",R,L);
            system("pause>nul");
            return L;
        }
        if(function==NULL)
        {
            printf("Error \"Ternary Search\":TSmin f=[NULL] is not allowed!\n");
            system("pause>nul");
            return L;
        }
        while(fabs(L-R)>eps)
        {
            double m1=L+(R-L)/3;
            double m2=R-(R-L)/3;
            if(function(m1)<function(m2))
                R=m2;
            else
                L=m1;
        }
        return L;
    }

    double TSmax(double L,double R,double (*function)(double x),double eps=1e-7)
    {
        //make a ternary search for max value and return x
        if(R==L)
            return L;
        if(R<L)
        {
            printf("Error \"Ternary Search\":TSmin R=%d (<L=%d) is not allowed!\n",R,L);
            system("pause>nul");
            return L;
        }
        if(function==NULL)
        {
            printf("Error \"Ternary Search\":TSmin f=[NULL] is not allowed!\n");
            system("pause>nul");
            return L;
        }
        while(fabs(L-R)>eps)
        {
            double m1=L+(R-L)/3;
            double m2=R-(R-L)/3;
            if(function(m1)>function(m2))
                R=m2;
            else
                L=m1;
        }
        return L;
    }
}

Treap的基本板子

详见 Treap树堆的基本模板（无讲解）

后记

算法模板库仍在扩充中…

GGN_2015
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
大模型理论基础--大模型的数据
Runnymmede的博客
06-27 20
到目前为止，我们已经讨论了大型语言模型的行为（能力和损害）。现在，我们要剥开洋葱的第一层，开始讨论这些模型是如何构建的。任何机器学习方法的起点都是训练数据，因此这就是我们开始的地方。附录：通常在机器学习中，训练数据和测试（评估）数据是相似的，或者至少是同一类型的。但对于大型语言模型来说，训练数据就是“原始文本”。
XMind总结C++标准IO库
06-02
C++标准IO库，通过C++标准输入输出流对控制窗口，磁盘文件字符串流进行控制。内容也包含条件状态，文件模式等内容。
Premier Bob的算法模板库
总理同学的编程尝试
04-07 1130
我是Premier Bob，这是我的在线算法模板库。热烈欢迎NEYC和NEYCBS的同学们来此参观。
Premier Bob的博客文章目录
总理同学的编程尝试
03-28 988
前言 同学们好，我是总理，一名普通的YCOIer，欢迎NEYC及NEYCBS的同学们访问我的博客。 另外，感谢带我入OI之门的Ingu老师，感谢邱老师，感谢培训期间机房的学长们，感谢所有支持我参加OI的老师和同学。 特别鸣谢：我的爷爷、奶奶，以及双语初中班主任王老师。 此页面仍在继续完善，敬请期待。 当前网址：http://blog.csdn.net/ggn_2015/ar
Spark Graphx图计算之二跳邻算法实战！
邓佑权的博客
11-17 3452
Spark Graphx图计算之二跳邻算法实战！ def sendMsgFunc(edge:EdgeTriplet[Int, Int]) = { if(edge.srcAttr <= 0){ if(edge.dstAttr <= 0){ // 如果双方都小于0，则不发送信息 Iterator.empty
看书
bb8cbb8c的专栏
12-02 1万+
<br />驱动器 D 中的卷是 新加卷<br /> 卷的序列号是 FC1A-1944<br /> D:/编程/整理编程文档/80 O'Reilly Books 的目录<br />2008/04/20  21:21       245,106,438 80.本.rar<br />               1 个文件    245,106,438 字节<br /> D:/编程/整理编程文档/ajax 的目录<br />2007/05/14  23:49         2,597,610 Ajax In
Datawhale9月组队学习Day4
最新发布
royi_0825的博客
09-14 130
昨天晚上突然发现自己之前参加datawhale拿了优秀学习者的比赛，自己得了二等奖……证书忘记申请了。我当时真没想到自己能得奖，毕竟CV好卷。昨晚过去瞅了一眼发现自己竟然得奖了。说回这边的大模型学习，昨天做的笔记链接，我最后都没有放入打卡里面，因为实在是很丢人，把原本的笔记给照抄照搬了一份过来。今天早上看新一章节的时候就在想要怎么写笔记。我想应该还是“尝试用自己的话来概括”，第二是仔细思考能从里面学习到，或者有什么疑问吧。
与大家分享一些计算机方面的电子书籍
angang6606的博客
08-31 1004
---update at 2005-12-7 地址已经改变，请重新来mail索取 ---update at 2006-3-13 个人的服务器有点问题，书籍暂时停止下载。 ---update at 2006-3-23 服务器恢复正常，欢迎大家下载，可以用留言或者mail的方式索取下载地址。 --update at 2006-7-10 近来发现有人开着多个线程从服务器下载，...
installshield 2015 Premier
12-18
installshield 2015 Premier + Crack
avast premier 高级版 授权文件
04-25
avast premier 高级版 授权文件 avast premier 高级版 授权文件 有效期到2021
InstallShield 2015 Premier 破解版
05-29
installshield 2015 Premier包含普通版本所没有的功能，如：可以设置安装文件夹权限... 压缩包包含使用说明
Primer Premier5.0
11-22
Primer Premier 5.0是一款广泛使用的分子生物学软件，主要用于设计PCR（聚合酶链式反应）引物和探针。这款软件集成了多种生物信息学工具，帮助科研人员在基因克隆、基因表达分析、SNP（单核苷酸多态性）检测等领域...
FFT详解
热门推荐
总理同学的编程尝试
03-31 9万+
一直想学FFT，当时由于数学基础太差，导致啥都学不懂。请教了机房里的几位学长大神，结果还是没太明白。因此下定决心写一篇关于“FFT”的文章，一篇起码我能看得懂的“FFT”。
（初学者）素数计算法
总理同学的编程尝试
03-29 8302
前言记得我刚开始学习OI的时候接触到的最早的算法莫过于素数算法。当然了，素数这个知识小学的时候我们就已经有所学习了。当时我一直有一个疯狂的设想，就是编写一个程序把我能算出来的素数都算出来。
平面最近点对的分治做法及其证明
总理同学的编程尝试
06-23 7381
2018.6.23 好久没写博客了，做了一道有趣的分治题，写个博客。 题目传送门： P1429 平面最近点对（加强版） 题意 题目描述 给定平面上n个点，找出其中的一对点的距离，使得在这n个点的所有点对中，该距离为所有点对中最小的 输入格式 第一行：n；2≤n≤200000 接下来n行：每行两个实数：x y (0≤x,y≤10^9)，表示一个点的行坐标和列坐标，中间用一个空格...
补充——FFT中的二进制翻转问题
总理同学的编程尝试
04-07 6033
注在我之前的FFT博客中好像并没有具体地介绍FFT中的“二进制翻转是如何实现的”，在这里我们简短的讨论一下。
线性时间复杂度排序
总理同学的编程尝试
04-05 5174
前言 我的博客之所以叫做“总理同学的暴力编程尝试…”说来话长。在我无数次的暴力编程尝试中，其中有一次就是给500万个整数进行排序，快排还有归并排序显然是卡爆了。快排一直在运转没有头，而归并排序直接就RE程序崩溃了。后来我听说了一种叫做“计数排序”的排序方式，得到了令我十分惊讶的成果，我用了4秒钟生成了500万个随机数，而给这些随机数排序却只用了3秒钟。因此怀着对“线性时间发杂度排序”的敬仰之情，
C++函数模板详解：提升算法设计效率
本章节探讨的是C++中的函数模板，这是C++编程中的一个重要概念，特别是在算法设计中。函数模板是一种通用编程工具，允许程序员为不同数据类型定义一个函数，而无需为每种类型分别编写函数。原始标题提到的min()函数...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值