题目描述

Blinker 有非常多的仰慕者,他给每个仰慕者一个正整数编号。而且这些编号还隐藏着特殊的意义,即编号的各位数字之积表示这名仰慕者对Blinker的重要度。 现在Blinker想知道编号介于某两个值A,B之间,且重要度为某个定值K的仰慕者编号和。


输入格式

输入的第一行是一个整数N,表示Blinker想知道的信息个数。 
接下来的N行,每行有三个数,A,B,K。表示 Blinker想知道编号介于A和B之间的,
重要度为K的仰慕者的编号和。


输出格式

输出N行,每行输出介于A和B之间,重要度为 K的仰慕者编号和。结果可能很大,
模上20120427。


提示

【数据范围】 
对于20%的数据,保证:  2<=A<=B<=1000000000,1<=N<=30 
对于50%的数据,保证:2<=A<=B<=1000000000000000000,1<=N<=30 
对于100%的数据,保证:  2<=A<=B<=1000000000000000000,1<=N<=5000


  • 题解:

    • 由于因子只有2 3 5 7 , 状态大概有3e4+,不超过4e4;
    • 用HASH_TABLE存下来;
    • 多组询问,先预处理$i$位的答案;
    • 对每组询问$dfs$统计答案,如果没有前导零并且没有抵达上界则直接返回,否则进一步$dfs$
    • 代码不太好些,注释里有关于变量的说明;
    • 注意$K=0$和$K \neq 0$最好分开统计;
    • $n$为位数上界,时间复杂度:$O(T*10n +n*4e4*10)$
 #include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=,M=4e4,sz=,mod=;
ll A,B,K;
int cnt,a[N],pw[N];
int o,hd[sz],nt[sz];ll v[sz];
void upd(int&x,int y){x+=y;if(x>=mod)x-=mod;}
void add(ll x){
int u=x%sz;
for(int i=hd[u];i;i=nt[i])if(v[i]==x)return;
nt[++o]=hd[u],hd[u]=o,v[o]=x;
}
int find(ll x){
int u=x%sz;
for(int i=hd[u];i;i=nt[i])if(v[i]==x)return i;
return ;
}//HASH表
void dfs(int x,int y,ll z){
if(x==)add(z);
else{
if(!y)return;
for(int i=x;i<=;++i){
dfs(i,y-,z);
z*=y;
}
}
}//dfs方案
namespace Solve{
int f[N][M],g[N][M];//f[位数][乘积]
void init(){
f[][]=;
for(int i=;i<;++i)
for(int j=;j<=o;++j)
if(f[i][j]){
for(int l=;l<=;++l){
int y=find(v[j]*l);
upd(f[i+][y],f[i][j]);
upd(g[i+][y],(g[i][j]*+f[i][j]*l)%mod);
}
}
}//预处理i位无限制的答案
int cal1(int i,int j,int k,ll x,ll y){//已经枚举好的位数,前导零,上界,前i位的值,当前K
if(i==cnt)return y==?x:;
if(!k&&!j)return (1ll*x*pw[cnt-i]%mod*f[cnt-i][find(y)]%mod+g[cnt-i][find(y)])%mod;
int re=,l=j?:,r=k?a[i+]:;
for(int i1=l,t;i1<=r;++i1)if(!i1||y%i1==){
ll y1=!i1?y:y/i1;
upd(re, t=cal1(i+,j&&!i1,k&&i1==a[i+],(x*+i1)%mod,y1));
}
return re;
}
int cal0(int i,int j,int k,ll x,ll y){//前四个同cal1,y为前i为是否出现0
if(i==cnt)return !y?x:;
if(!k&&!j){
if(!y)return (1ll*x*pw[cnt-i]%mod*pw[cnt-i]%mod+1ll*(pw[cnt-i]-)*pw[cnt-i]/%mod)%mod;//后面cnt-i任意;
else return 1ll*x*pw[cnt-i]%mod*f[cnt-i][o]%mod+g[cnt-i][o];//后面cnt-i为乘积为0;
}
int re=,l=,r=k?a[i+]:;
for(int i1=l;i1<=r;++i1){
int y1=j&&!i1?:y*i1;
upd(re, cal0(i+,j&&!i1,k&&i1==a[i+],(x*+i1)%mod,y1));
}
return re;
}
//注意理解一下两个方面:
//是如何统计前i位达到上界,i+1位不是上界的答案的;(状态k)
//是如何统计位数小于上界位数cnt的答案的; (状态j)
//上界是如何被统计的(line44)
}
int main(){
#ifndef ONLINE_JUDGE
freopen("T2.in","r",stdin);
freopen("T2.out","w",stdout);
#endif
dfs(,,);add();
Solve::init();
for(int i=pw[]=;i<=;++i)pw[i]=pw[i-]*%mod;
int T;scanf("%d",&T);
while(T--){
int ans=;
scanf("%lld%lld%lld",&A,&B,&K);A--;
cnt=;while(A)a[++cnt]=A%,A/=;
for(int i=;i<=cnt>>;++i)swap(a[i],a[cnt-i+]);
ans-= K?Solve::cal1(,,,,K) : Solve::cal0(,,,,);
cnt=;while(B)a[++cnt]=B%,B/=;
for(int i=;i<=cnt>>;++i)swap(a[i],a[cnt-i+]);
ans+= K?Solve::cal1(,,,,K) : Solve::cal0(,,,,);
printf("%d\n",(ans%mod+mod)%mod);
}
return ;
}

bzoj2757【scoi2012】Blinker的仰慕者的更多相关文章

  1. BZOJ2757 : [SCOI2012]Blinker的仰慕者

    BZOJ AC900题纪念~~ 若K>0,则 设f[i][j]表示i位数字,积为j的数字的个数 g[i][j]表示i位数字,积为j的数字的和 DP+Hash预处理 查询时枚举LCP然后统计贡献 ...

  2. BZOJ2758 : [SCOI2012]Blinker的噩梦

    首先将包含关系建树. 方法是将每个图形拆成上半边和下半边,从左往右扫描线,用Splay从下到上维护扫描线上所有图形. 每次加入一个新的图形$x$的时候,看看它下方第一个图形$y$,如果$y$是上半边, ...

  3. [BZOJ2758] [SCOI2012]Blinker的噩梦 扫描线+set

    题目大意:有n个圆或凸多边形,这些图形不会相交,每当走入或走出一个图形时需要异或上一个代价,有m组操作: 询问操作,每次询问从一个点走到另一个点时,需要的代价(初始代价为0) 修改操作,每次修改一个图 ...

  4. bzoj AC倒序

    Search GO 说明:输入题号直接进入相应题目,如需搜索含数字的题目,请在关键词前加单引号 Problem ID Title Source AC Submit Y 1000 A+B Problem ...

  5. 近期概况&总结

    下午考完英语的学考就要放假啦,是衡中的假期啊QAQ 所以灰常的激动,一点也不想写题(我不会告诉你其实假期只有一个晚上.. 自从CTSC&APIO回来之后就一直在机房颓颓颓,跟着zcg学了很多新 ...

  6. bzoj2758【scoi2012】Blinker的的噩梦

    题目描述 一天Blinker醒来,发现自己成为了一个二维世界的点,而且被标记上了一个奇怪的值. 这个世界是由N个边界互不相交(且不相切)的图形组成,这里图形仅包括圆和凸多边形.每个图形还有一个权值.每 ...

  7. BZOJ 2756: [SCOI2012]奇怪的游戏 [最大流 二分]

    2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 3352  Solved: 919[Submit][Stat ...

  8. Bzoj2756 [SCOI2012]奇怪的游戏

    2756: [SCOI2012]奇怪的游戏 Time Limit: 40 Sec  Memory Limit: 128 MBSubmit: 3220  Solved: 886 Description ...

  9. 【BZOJ】【2756】【SCOI2012】奇怪的游戏

    网络流-最大流 这题……建模部分先略过 这道题是会卡时限的T_T俺的Dinic被卡了,在此放几篇很棒的讲网络流算法的文章,至于大家耳熟能详的论文就不放了…… http://www.cppblog.co ...

随机推荐

  1. django-simple_tag、filter

    simple_tag与filter的用法 1.支持自定义函数处理方法 2.支持模板调用 创建步骤: a.在app目录下创建templatetags文件夹 b.在templatetags中创建任意名称. ...

  2. 【Docker】第三篇 Docker容器管理

    一.Docker容器概述: 简单理解容器是镜像的一个实例. 镜像是静态的只读文件,而容器的运行需要可写文件层. 二.创建容器 [root@web130 ~]# docker create -it ub ...

  3. 从零开始的Python学习Episode 22——多线程

    多线程 线程 线程是操作系统能够进行运算调度的最小单位.它被包含在进程之中,是进程中的实际运作单位.一条线程指的是进程中一个单一顺序的控制流,一个进程中可以并发多个线程,每条线程并行执行不同的任务. ...

  4. [SimHash] the Hash-based Similarity Detection Algorithm

    The current information explosion has resulted in an increasing number of applications that need to ...

  5. Python基础_可迭代的/迭代器/生成器

    介绍 可迭代的:内部实现了__iter__方法 迭代器:内部实现了__iter__,__next__方法 生成器:yield,yield from 使用 __iter__() __next__() _ ...

  6. 部署Maven项目到tomcat报错:java.lang.ClassNotFoundException: org.springframework.web.context.ContextLoaderListener【转】

    部署Maven项目到tomcat报错:java.lang.ClassNotFoundException: org.springframework.web.context.ContextLoaderLi ...

  7. Bing词典vs有道词典比对测试报告——功能篇之核心功能

    必应词典vs有道词典 核心功能对比 从应用的UI布局来看,这两款软件的功能如下: 相同 不同 必应词典 词典.例句.翻译 百科 有道词典 词典.例句.翻译 应用 就词典类软件来说,词典是最核心的功能. ...

  8. 08慕课网《进击Node.js基础(一)》事件events

    引用events模块中的EventEmitter 事件的监听和发射 相同的事件发射数量有限,可以通过setMaxListeners设置峰值 var EventEmitter = require('ev ...

  9. vue 里面的watch 选项详解

    早就想好好写写这个watch了,一直顾不上,因为想深刻的记录一下,其实这些东西只要是好好看看官网的说明,都在里面呢. 虽然计算属性在大多数情况下更合适,但有时也需要一个自定义的侦听器.这就是为什么 V ...

  10. vue2.0 keep-alive 最佳实战(转载)

    1.基本用法 vue2.0提供了一个keep-alive组件用来缓存组件,避免多次加载相应的组件,减少性能消耗 <keep-alive> <component> <!-- ...