topcoder srm 697 div1 -3

topcoder srm 697 div1 -3

1、给定长度为$n$ 的数组$b$,构造长度为$n$ 的且没有重复元素的数组$a$,令$p_{i}$表示$a$中除$a_{i}$外其他元素的乘积。构造出的$a$满足$a_{i}^{b_{i}}$能够被$p_{i}$整除。这样的数组$a$存在否?

思路:因为$a_{i}^{b_{i}}$中所有素因子的指数大于等于$p_{i}$。所以可以假设所有的$a_{i}$只含有素因子2。令$a_{i}=2^{A_{i}},s=\sum_{i=1}^{n}A_{i}$,那么有$A_{i}b_{i}\geq s-A_{i}$,即$A_{i}\geq \frac{s}{1+b_{i}}$。所有这些加起来约去$s$得到$\sum_{i=1}^{n}\frac{1}{1+b_{i}}\leq 1$。令$K=\sum_{i=1}^{n}\frac{1}{1+b_{i}}$。所以$K>1$无解。$K=1$但是存在相同的$b$时无解。因为$K=1$说明$A_{i}= \frac{s}{1+b_{i}}$,但是相同的$b$ 导致出现了相同的$A$。

#include <stdio.h>
#include <string>
#include <stack>
#include <vector>
#include <string.h>
#include <algorithm>
#include <assert.h>
using namespace std;



class DivisibleSetDiv1
{
public:
	string isPossible(vector <int> b)
	{
	    const int n=(int)b.size();
	     sort(b.begin(),b.end());
	     double ans=0;
	    for(int i=0;i<n;++i) ans+=1.0/(b[i]+1);

	    if(ans>1+1e-10) {
            return "Impossible";
	    }
	    else if(ans>1-1e-10) {
            for(int i=1;i<n;++i) if(b[i]==b[i-1]) {
                return "Impossible";
            }
	    }
        return "Possible";
	}
};

  

2、有$n$匹马,每个马有个各不相同的值$a_{i}$。有$2^{m}$场比赛。第$i$匹马在第$j$场比赛中的用时为$a_{i}$^$j$(抑或)。每场比赛所有的马排名为$[0,n-1]$。排名为$k$的马得到的罚时为$k^{2}$。所有比赛结束后每匹马的所有罚时加起来得到第$i$匹马的总罚时$p_{i}$,计算所有$p_{i}$%1000000007的抑或值。

思路:将每匹马的值$a_{i}$建立一课二叉树,并在节点上记录子树中插入的马的个数。然后依次计算每一匹马的$p_{i}$值。对于第$i$匹马,设比赛$j$的二进制高$t$位等于$a_{i}$的高$t$位的所有比赛中它的罚时为$g=\sum_{r=1}^{2^{m-t}}s_{r}^{2}$。$2^{m-t}$是因为已经进行了这么多长比赛。那么在进行比赛的二进制第$m-t$位不同于$a_{i}$的第$m-t$位但是它们的高$t-1$ 位相同时的所有比赛中,它的罚时为$g^{‘}=\sum_{r=1}^{2^{m-t}}(s_{r}+x)^{2}$。其中$x$为与$a_{i}$ 的第$m-t$位不同的一侧所有马的个数。那么$g^{‘}=g+2^{m-t}x^{2}+2x\sum_{r=1}^{2^{m-t}}s_{r}$。这样记录已经遍历的子树的前缀和$\sum_{r=1}^{2^{m-t}}s_{r}$即可。

#include <stdio.h>
#include <string>
#include <stack>
#include <vector>
#include <string.h>
#include <algorithm>
#include <assert.h>
using namespace std;

const int N=200005;
const int mod=1000000007;

int a[N];

int T[N*32][2];
int S[N*32];
int e;
int m;


void add(int x)
{
    int p=1;
    for(int i=m-1;i>=0;--i)
    {
        int t=(x>>i)&1;
        if(!T[p][t]) T[p][t]=++e;
        p=T[p][t];
        ++S[p];
    }
}

int c[32];

void cal(int x)
{
    int p=1;
    for(int i=m-1;i>=0;--i)
    {
        int t=(x>>i)&1;
        c[i]=S[T[p][t^1]];
        p=T[p][t];
    }
}

class XorOrderDiv1
{
public:
    int get(int M,int n,int a0,int b)
    {
        m=M;
        a[1]=a0;
        for(int i=2;i<=n;++i) a[i]=(a[i-1]+b)&((1<<m)-1);
        e=1;
        for(int i=1;i<=n;++i) add(a[i]);

        int ans=0;
        for(int i=1;i<=n;++i)
        {
            cal(a[i]);
            int tmp=0,sum=0;
            for(int j=0;j<m;++j)
            {
                tmp=(tmp+(long long)(1<<j)*c[j]%mod*c[j]%mod+(long long)2*c[j]%mod*sum%mod+tmp)%mod;
                sum=(sum*2+(long long)c[j]*(1<<j))%mod;
            }
            ans^=tmp;
        }
        return ans;
    }
};

  

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/101196.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • plsqldev 日期格式

    plsqldev 日期格式之前装win7+oracle11R2(64)+ instantclient_11_2(32)+PLSQL(32)费了很大力气,见前一个帖子,后果就是plsql启动时读的环境变量位置是五花八门,可能是注册表中oraclehone下的,也可能是instantclient下的或者是电脑高级属性中环境变量,当然start.bat中的设置优先。  plsql中执行以下语…

    2022年5月10日
    83
  • 黎曼猜想和素数分布的关系_黎曼公式和素数的关系

    黎曼猜想和素数分布的关系_黎曼公式和素数的关系自然数简化到素数:黎曼猜想RiemannHypothesis及其解释(公号回复“黎曼猜想”下载PDF经典收藏版彩标资料)原创:秦陇纪数据简化DataSimp今天数据简化DataSimp导读:科学大院《黎曼猜想RiemannHypothesis简介》来自黄逸文(中国科学院数学与系统科学研究院),介绍了黎曼猜想RiemannHypothesis大概。知乎译文《黎曼猜想Riemann…

    2022年8月11日
    42
  • Kali linux新手入门视频教程|Kali linux安装

    Kali linux新手入门视频教程|Kali linux安装一、 Kalilinux是什么?KaliLinux是基于Debian的Linux发行版,设计用于数字取证操作系统。KaliLinux面向专业的渗透测试和安全审计.因此,KaliLinux已经进行了如下的多处核心的修改。单用户,设计成root权限登录:由于安全审计的本质,KaliLinux被设计成使用”单用户,root权限“方案。二、 Kalilinux新手入门教程目录(视频教程)…

    2022年5月26日
    65
  • 网上书城项目总结(servlet_jsp+javaBean)

    网上书城项目总结(servlet_jsp+javaBean)          网上书城项目总结1项目大纲设计:需求分析系统设计详细设计权限设计2技术选型:Servlet+jsp+javaBeanListener+Filter+jstl+fileupload+c3p0+dbutils+mysql3开发顺序:  从dao层到service层再到web层网上书城需求分析:分别对管理员,普通用户,系统三个用户…

    2022年7月27日
    17
  • 提升进程权限-OpenProcessToken等函数的用法[通俗易懂]

    提升进程权限-OpenProcessToken等函数的用法[通俗易懂]提升进程权限文章一:在枚举/结束系统进程或操作系统服务时,会出现自己权限不足而失败的情况,这时就需要提升自己进程到系统权限,其实提升权限的代码很简单的,看到过的最经典的应该是《WINDOWS核心编程》第四章中操作进程给出的那个函数了,如果我们真的不了解它的操作也不要紧,因为只要在你需要的地方调用下面这个函数就是了,以下是它的代码:BOOLEnablePriv(){HAND

    2022年6月25日
    48
  • 角点检测方法_什么叫五点取样法

    角点检测方法_什么叫五点取样法目录原理讲解【1】为何选取角点作为特征?【2】角点的定义:【3】判断角点的方法:【4】Harris角点检测法示例Opencv自带函数:cornerHarris()函数示例程序1示例程序2原理讲解【1】为何选取角点作为特征?角点是一种局部特征。角落上的可区分性特别强,边缘次之,平滑区域则基本没有区分性。【2】角点的定义:【3】判断角点的方法:这里有个细节:将计算的所有方向上的变化……

    2022年10月21日
    3

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号