cogs 141. [USACO Jan08] 奶牛的选举

★   输入文件:elect.in   输出文件:elect.out   简单对比
时间限制:1 s   内存限制:16 MB

在推翻了Farmer John这个残暴的统治者后,奶牛们举行了她们的第一次总统大选,贝茜也是N(1 <= N <= 50,000)头候选奶牛之一。不过,作为一头有远见的奶牛,贝茜想在选举开始前就计算出,哪头奶牛最有可能在竞争中胜出。

选举分两轮进行。第一轮中,得票最多的K(1 <= K <= N)头奶牛晋级到下一轮,在第二轮选举中得票最多的奶牛成为最终的总统。

现在,贝茜告诉了你奶牛i在第一轮投票中的期望得票数A_i(1 <= A_i <= 1,000,000,000)以及她在第二轮投票中的期望得票数B_i(1 <= B_i <= 1,000,000,000)(如果奶牛i能成功晋级的话),她希望你帮她计算一下,如果这些数据无误,那么哪头奶牛将成为总统。任何数值都不会在A_i 列表中出现两次,在B_i列表中也是如此。

程序名: elect

输入格式:

  • 第1行: 2个用空格隔开的整数:N 和 K
  • 第2..N+1行: 第i+1为2个用空格隔开的整数:A_i 和 B_i

输入样例 (elect.in):

5 3
3 10
9 2
5 6
8 4
6 5

输入说明:

一共有5头奶牛参加选举,在第一轮中得票最多的3头奶牛可以晋级至第二轮。奶牛们在第一轮中的得票期望分别为3,9,5,8,6,第二轮中,分别为10,2,6,4,5。

输出格式:

  • 第1行: 输出1个整数,为将被选为总统的奶牛的编号

输出样例 (elect.out):

5

输出说明:

奶牛2,4,5晋级到第二轮。奶牛5在第二轮投票中得到了最多的5票,赢得了选举的最终胜利。

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<algorithm>
 4 #define ll long long
 5 
 6 using namespace std;
 7 const ll N=50001;
 8 
 9 ll a[N];
10 
11 struct node{
12     ll me;
13     ll cmaq;
14     ll cmah;
15 }E[N];
16 
17 inline void read(ll & x) 
18 {
19     char c=getchar();
20     x=0;
21     while(c<'0'||c>'9')c=getchar();
22     while(c>='0'&&c<='9')x=x*10+c-'0',c=getchar();
23 }
24 
25 inline bool cmp1(node a,node b)
26 {
27     return a.cmaq>b.cmaq;
28 }
29 
30 inline bool cmp2(node a,node b)
31 {
32     return a.cmah>b.cmah;
33 }
34 
35 int main()
36 {
37     freopen("elect.in","r",stdin);
38     freopen("elect.out","w",stdout);
39     ll n,k;
40     read(n);
41     read(k);
42     for(int i=1;i<=n;i++)
43     {
44         read(E[i].cmaq);
45         read(E[i].cmah);
46         E[i].me=i;
47     }
48     sort(E+1,E+n+1,cmp1);
49     sort(E+1,E+k+1,cmp2);
50     
51     printf("%d",E[1].me);
52     return 0;
53 }

 

转载于:https://www.cnblogs.com/lyqlyq/p/6863342.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是将代码修改为cot平滑的方法: 1. 首先,需要使用边界角的cot权重计算每个顶点的权重。 2. 然后,使用cot权重对每个顶点的邻域点进行加权计算,得到平滑后的坐标。 3. 最后,根据平滑后的坐标更新每个顶点的位置。 修改后的代码如下: float smooth() { float err = -1; cogs.clear(); v_end = mesh.vertices_end(); //cot平滑 for (v_it = mesh.vertices_begin(); v_it != v_end; ++v_it) { cog[0] = cog[1] = cog[2] = weight_sum = 0.0; for (vv_it = mesh.vv_iter(*v_it); vv_it.is_valid(); ++vv_it) { double cot_weight = 0.0; MyMesh::HalfedgeHandle heh = mesh.find_halfedge(*v_it, *vv_it); if (!mesh.is_boundary(heh)) { MyMesh::HalfedgeHandle prev_heh = mesh.prev_halfedge_handle(heh); MyMesh::HalfedgeHandle next_heh = mesh.next_halfedge_handle(heh); MyMesh::VertexHandle prev_vh = mesh.to_vertex_handle(prev_heh); MyMesh::VertexHandle next_vh = mesh.to_vertex_handle(next_heh); MyMesh::Point prev_p = mesh.point(prev_vh); MyMesh::Point curr_p = mesh.point(*v_it); MyMesh::Point next_p = mesh.point(next_vh); double cot_alpha = cot(prev_p - curr_p, next_p - curr_p); double cot_beta = cot(curr_p - prev_p, next_p - prev_p); cot_weight = cot_alpha + cot_beta; } cog += cot_weight * mesh.point(*vv_it); weight_sum += cot_weight; } cogs.push_back(cog / weight_sum); } for (v_it = mesh.vertices_begin(), cog_it = cogs.begin(); v_it != v_end; ++v_it, ++cog_it) { if (!mesh.is_boundary(*v_it)) { MyMesh::Point p = mesh.point(*v_it); err = max(err, (p - *cog_it).norm()); mesh.set_point(*v_it, *cog_it); } } return err; } 其中cot函数的定义如下: double cot(MyMesh::Point a, MyMesh::Point b) { return dot(a, b) / cross(a, b).norm(); } 注意,这里使用的是边界角的cot权重,因此在计算cot权重时需要判断当前边是否为边界。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值