关闭

[LibreOJ β Round]ZQC的树列

313人阅读 评论(0) 收藏 举报
分类:

题目描述

见LOJ。

做法

显然对于每一段数字,要么可选要么必选。
假如有t个,若必选,则方案为2^t-1,否则为2^t。
因此只有1、2、3就可以构造所有情况。
尝试把n拆分成2^t-1以及2^t的乘积。
这个用dfs即可。

#include<cstdio>
#include<algorithm>
#include<cmath>
#define fo(i,a,b) for(i=a;i<=b;i++)
using namespace std;
typedef long long ll;
ll a[100000+10],sta[100000+10],b[100000+10];
ll i,j,k,l,t,m,tot,top,now,cnt;
ll n;
bool czy;
ll random(ll x){
    ll t=rand()%10000;
    t=t*10000+rand()%10000;
    t=t*10000+rand()%10000;
    t=t*10000+rand()%10000;
    return t%x;
}
ll mul(ll a,ll b,ll p){
    if (!b) return 0;
    ll t=mul(a,b/2,p);
    t=(t+t)%p;
    if (b%2) t=(t+a)%p;
    return t;
}
ll qsm(ll x,ll y,ll mo){
    if (!y) return 1;
    ll t=qsm(x,y/2,mo);
    t=mul(t,t,mo);
    if (y%2) t=mul(t,x,mo);
    return t;
}
bool Miller_Rabin(ll n){
    if (n==1) return 0;
    int s=5,t=0,i;
    ll a,p,k=n-1;
    while (k%2==0) k/=2,t++;
    while (s--){
        a=random(n-1)+1;
        p=a=qsm(a,k,n);
        fo(i,1,t){
            a=mul(a,a,n);
            if (a==1&&p!=1&&p!=n-1) return 0;
            p=a;
        }
        if (a!=1) return 0;
    }
    return 1;
}
ll gcd(ll a,ll b){
    return b?gcd(b,a%b):a;
}
ll Pollard_Rho(ll n){
    ll k,x,y,c,d,i=1;
    while (1){
        c=random(n-1);
        k=2;y=x=random(n);
        i=1;
        while (1){
            y=(mul(y,y,n)+c)%n;
            d=gcd(abs(x-y),n);
            if (i==k) x=y,k<<=1;
            i++;
            if (d!=1) break;
        }
        if (d!=n) return d;
    }
}
void work(ll n){
    if (Miller_Rabin(n)){
        sta[++top]=n;
        return;
    }
    ll p=Pollard_Rho(n);
    ll q=n/p;
    work(p);work(q);
}
bool pd(ll n){
    while (n%2==0) n/=2;
    if (n==1) return 1;else return 0;
}
void dfs(ll n){
    if (pd(n)){
        cnt=0;
        while (n%2==0) cnt++,n/=2;
        czy=1;
        return;
    }
    ll t=4,l=2;
    while (t-1<=n){
        if (n%(t-1)==0){
            top++;
            b[top]=l;
            dfs(n/(t-1));
            if (czy) return;
            top--;
        }
        t*=2;
        l++;
    }
}
int main(){
    czy=0;
    scanf("%lld",&n);
    if (n==1){
        printf("1\n");
        printf("1\n");
        return 0;
    }
    dfs(n);
    /*t=0;
    a[tot=1]=1;
    now=3;
    fo(i,1,top)
        if (sta[i]==2) t++;
        else{
            k=sta[i]+1;
            l=0;
            while (k>1){
                if (k%2==1) break;
                k/=2;
                l++;
            }
            if (k>1){
                czy=0;
                break;
            }
            while (l){
                a[++tot]=now;
                l--;
            }
            if (now==1) now=3;else now=1;
        }
    if (!czy){
        printf("qnq\n");
        return 0;
    }
    while (t){
        a[++tot]=2;
        t--;
    }
    a[++tot]=now;*/
    if (!czy){
        printf("qnq\n");
        return 0;
    }
    a[tot=1]=1;
    now=3;
    fo(i,1,top){
        k=b[i];
        while (k){
            a[++tot]=now;
            k--;
        }
        if (now==1) now=3;else now=1;
    }
    while (cnt){
        a[++tot]=2;
        cnt--;
    }
    a[++tot]=now;
    printf("%lld\n",tot);
    fo(i,1,tot) printf("%lld ",a[i]);
}
0
0
查看评论
发表评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场

「LibreOJ β Round #2」贪心只能过样例 [bitset]【STL】

题目链接:https://loj.ac/problem/515 ——————————————————————————————————515. 「LibreOJ β Round #2」贪心只能过样例内...
  • qq_33184171
  • qq_33184171
  • 2017-07-03 14:51
  • 386

LibreOJ β Round #2 A. 模拟只会猜题意【前缀和】

A. 模拟只会猜题意 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 转到题库提交记录返回比赛 ...
  • mengxiang000000
  • mengxiang000000
  • 2017-07-03 17:56
  • 453

LibreOJ β Round F. ZQC 的游戏【思维建图+最大流】

F. ZQC 的游戏 内存限制: 256 MiB时间限制: 1000 ms 标准输入输出 转到题库提交记录返回比赛 题目描述 Agar.io 是一款流行的游戏,每个玩家在二维...
  • mengxiang000000
  • mengxiang000000
  • 2017-06-17 11:57
  • 387

「LibreOJ β Round #2」DP 一般看规律

set,map,启发式合并
  • KikiDMW
  • KikiDMW
  • 2017-07-04 17:39
  • 179

[数论] LOJ #510. 「LibreOJ NOI Round #1」北校门外的回忆

这个题跟树状数组没有半毛钱关系 首先这是一个最低位翻倍的过程,如果这一位最终会变成 00 ,那么步数是 O(logn)O(\log n)的? 要是这一位不能变成 00 ,也就是在环上跑了,似乎跑到...
  • u014609452
  • u014609452
  • 2017-07-11 22:29
  • 250

LibreOJ【配对游戏】

有 n 个人排成一排,一开始全部面向前方,然后随机朝左或是朝右转。 然后我们不断审查这个队列,每次选择两个面对面的相邻的人,将他们从队列中取出。 例如(> 表示向右, 表示向左): 队列 >>>...
  • ouqingliang
  • ouqingliang
  • 2017-07-11 15:20
  • 118

LibreOJ β Round #2 B.贪心只能过样例【Dp+bitset优化】

B. 贪心只能过样例 内存限制:256 MiB时间限制:1000 ms标准输入输出 题目类型:传统评测方式:文本比较 转到题库提交记录返回比赛 ...
  • mengxiang000000
  • mengxiang000000
  • 2017-07-03 17:53
  • 345

[数论] LOJ #508. 「LibreOJ NOI Round #1」失控的未来交通工具

这种非简单路 一般转化成任意一条路加上若干环 这里大概是任意一条路加上若干环长的gcd 任意一条路 可以弄出任意一颗生成树? 但是这里实际上只需要一个带权并查集 详见官方题解#include #...
  • u014609452
  • u014609452
  • 2017-07-23 13:26
  • 227

516. 「LibreOJ β Round #2」DP 一般看规律 [set/SPLAY] 【STL/数据结构】

题目链接:https://loj.ac/problem/516 —————————————————————————————————————————— 516. 「LibreOJ β Round #...
  • qq_33184171
  • qq_33184171
  • 2017-07-04 21:01
  • 166

LibreOJ NOIP Round #1 —— 蒟蒻总结&部分题解

LibreOJ NOIP Round #1 Day 1 Day 1 T1. DNA 序列 水题,最直观的解法是hash,发现k 定位:noip d1t1,送分题 #include...
  • youhavepeople
  • youhavepeople
  • 2017-11-08 17:23
  • 112
    个人资料
    • 访问:327186次
    • 积分:11634
    • 等级:
    • 排名:第1512名
    • 原创:814篇
    • 转载:4篇
    • 译文:0篇
    • 评论:201条
    最新评论
    文章分类