【bzoj1027】合金

分析

数形结合+计算几何+Floyd最小环。

http://blog.csdn.net/popoqqq/article/details/40539273

虽然这样占大家的很不好,但是我的确需要这么做。

小结

(1)最小环的写法

(2)凸包的一种简易算法

RT。

(3)关于计算几何的练习

通过做这道题,计算几何的模板也熟悉了一些。

计算几何如何调试?

这个么...好像只能瞪眼法了。

代码

#include <cstdio>
#include <climits>
#include <cmath>
#include <algorithm>
using namespace std;

#define rep(i,a,b) for (int i=(a);i<=(b);i++)

const int M=512;
const int N=512;
const int MAX=INT_MAX>>1;
const double EPS=1e-15;

int cmp(double x);

int m,n;
struct point
{
    double x,y;
    point(double _x=0,double _y=0) {
        x=_x,y=_y;
    }
    friend point operator - (point a,point b) {
        return point(a.x-b.x,a.y-b.y);
    }
    friend point operator * (point a,point b) {
        return point(a.x*b.x,a.y*b.y);
    }
    friend int operator == (point a,point b) {
        return !cmp(a.x-b.x)&&!cmp(a.y-b.y);
    }
    friend int operator != (point a,point b) {
        return !(a==b);
    }
}mer[M],prd[N];

int f[M][M];
int res;

void Init(void) {
    scanf("%d%d",&m,&n);
    rep(i,1,m) {
        double x,y; scanf("%lf%lf%*lf",&x,&y);
        mer[i]=point(x,y);
    }
    rep(i,1,n) {
        double x,y; scanf("%lf%lf%*lf",&x,&y);
        prd[i]=point(x,y);
    }
}

int cmp(double x)
{
    if(fabs(x)<EPS)  return 0;
    return x>0?1:-1;
}

double det(point a,point b) {
    return a.x*b.y-a.y*b.x;
}

int Left(point k,point a,point b) {
    double t=det(b-a,k-a);
    return t>0;
}
int On(point k,point a,point b) {
    if (cmp(det(k-a,b-a))!=0) return 0;
    if (!(cmp(min(a.x,b.x)-k.x)<=0&&cmp(k.x-max(a.x,b.x))<=0)) return 0;
    if (!(cmp(min(a.y,b.y)-k.y)<=0&&cmp(k.y-max(a.y,b.y))<=0)) return 0;
    return 1;
}

/*
bool On(const point &a,const point &b,const point &c) {
    return !cmp(det(a-b,c-b))&&cmp((a.x-b.x)*(a.x-c.x))<=0&&cmp((a.y-b.y)*(a.y-c.y))<=0;
}
*/

int Check(point a,point b) {
    rep(i,1,n) {
        int t1=Left(prd[i],a,b);
        int t2=On(prd[i],a,b);
        if (!t1&&!t2) return 0;
    }
    return 1;
}

void MakeGraph(void) {
    rep(i,1,m) rep(j,1,m) f[i][j]=MAX;
    rep(i,1,m) rep(j,1,m)
        if (Check(mer[i],mer[j]))
            f[i][j]=1;
}

int Floyd(void) {
    rep(i,1,m) {
        int mrk=1;
        rep(j,1,n) if (mer[i]!=prd[j]) {
            mrk=0;
            break;
        }
        if (mrk) return 1;
    }

    rep(i,1,m) rep(j,1,m) {
        int mrk=1;
        rep(k,1,n)
            if (!On(prd[k],mer[i],mer[j])) {
                mrk=0;
                break;
            }
        if (mrk) return 2;
    }

    /*
    rep(i,1,m) {
        rep(j,1,m)
            if (f[i][j]!=MAX)
                printf("%d ",f[i][j]);
            else printf("0 ");
        printf("\n");
    }
    */

    int ans=MAX;
    rep(k,1,m) rep(i,1,m) rep(j,1,m)
        f[i][j]=min(f[i][j],f[i][k]+f[k][j]);
    rep(i,1,m) ans=min(ans,f[i][i]);
    return ans!=MAX?ans:-1;
}

int main(void) {
    #ifndef ONLINE_JUDGE
    freopen("bzoj1027.in","r",stdin);
    freopen("bzoj1027.out","w",stdout);
    #endif

    Init();
    MakeGraph();

    res=Floyd();
    printf("%d\n",res);

    return 0;
}

【bzoj1027】合金的更多相关文章

  1. bzoj1027【JSOI2007】合金

    题目描述 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金 ...

  2. bzoj1027 [JSOI2007]合金

    1027: [JSOI2007]合金 Time Limit: 4 Sec  Memory Limit: 162 MBSubmit: 2671  Solved: 703[Submit][Status][ ...

  3. bzoj千题计划123:bzoj1027: [JSOI2007]合金

    http://www.lydsy.com/JudgeOnline/problem.php?id=1027 因为x+y+z=1,所以z=1-x-y 第三维可以忽略 将x,y 看做 平面上的点 简化问题: ...

  4. BZOJ1027 [JSOI2007]合金 【计算几何 + floyd】

    题目 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的 原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金的 ...

  5. [bzoj 1027][JSOI2007]合金(解析几何+最小环)

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1027 分析: 首先因为一个合金的和为1,所以考虑2个材料合金能否合成一个需求合金的时候 ...

  6. 合金装备V 幻痛 制作技术特辑

    合金装备V:幻痛 制作特辑 资料原文出自日版CGWORLD2015年10月号   在[合金装备4(Metal Gear Solid IV)]7年后,序章作品[合金装备5 :原爆点 (Metal Gea ...

  7. 【BZOJ】【1027】【JSOI2007】合金

    计算几何/凸包/Floyd Orz rausen大爷太强辣 计算几何题目果然不会做>_> 这个题……虽然他给了3个坐标,但实际上是个二维的计算几何题= =因为第三维坐标可以直接用前两维坐标 ...

  8. 1027: [JSOI2007]合金 - BZOJ

    Description 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的 ...

  9. BZOJ_1027_[JSOI2007]_合金_(计算几何+Floyd求最小环)

    描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1027 共三种金属,\(m\)种材料,给出每种材料中三种金属的占比. 给出\(n\)种合金的三种 ...

随机推荐

  1. winform 发邮件

    System.Net.Mail.SmtpClient client = new SmtpClient("smtp.163.com");   client.UseDefaultCre ...

  2. 抓取Js动态生成数据且以滚动页面方式分页的网页

    代码也可以从我的开源项目HtmlExtractor中获取. 当我们在进行数据抓取的时候,如果目标网站是以Js的方式动态生成数据且以滚动页面的方式进行分页,那么我们该如何抓取呢? 如类似今日头条这样的网 ...

  3. SqlSever基础 lower函数 返回字符串的小写形式

    镇场诗:---大梦谁觉,水月中建博客.百千磨难,才知世事无常.---今持佛语,技术无量愿学.愿尽所学,铸一良心博客.------------------------------------------ ...

  4. activiti学习总结

    Activiti界面元素的使用总结 一.图形设计中元素的使用 1.SequenceFlow:连接线,可以连接两个任务,来管理流程实例的流向 -----General -----id:流程的id,用与程 ...

  5. BASE64,MD5,SHA,HMAC加密與解密算法(java)

    package com.ice.webos.util.security; import java.io.UnsupportedEncodingException; import java.math.B ...

  6. 我的android学习经历15

    利用Intent实现有返回结果的页面跳转 主要用的方法: (1)Intent的构造方法:intent(当前界面对象,要跳转的界面.class); (2)接受结果的方法onActivityResult( ...

  7. BZOJ 3514 Codechef MARCH14 GERALD07加强版

    题目链接:http://www.lydsy.com:808/JudgeOnline/problem.php?id=3514 题意:给出一个图m条边.每次询问只加入编号在区间[L,R]之内的边有多少连通 ...

  8. Using Pre-Form Trigger In Oracle Forms

    Pre-Form trigger in Oracle Forms fires during the form start-up, before forms navigates to the first ...

  9. Hibernate 查询语句

    Hibernate查询语言(HQL)是一种面向对象的查询语言,类似于SQL,但不是对表和列操作,HQL适用于持久对象和它们的属性. HQL查询由Hibernate转换成传统的SQL查询,这在圈上的数据 ...

  10. 07.常用的SQL语句

    常用SQL语句,仿照着写(只改动字段和表名,库名) 查看数据库引擎                SHOW ENGINES;(默认InnoDB) 用数据库生成UUID            selec ...