题目:
题目描述
老虎和蒜头是好朋友。
秋天到了,森林中学又迎来了一年一度的招生季。
为了与科技时代接轨,森林中学今年即将使用全新的答题卡。
只要考试时学生对答题卡进行规范性的填涂,考试后不用多久就能得到大家考试的成绩了。
然而,被聘请为顾问的老虎却犯了难。原来,在数学考试中,往往会出现如 2/3,√23这样的答案,
而这样的答案在答题卡上没有填涂的方案——答题卡上只能够填涂在 0 至 999 之间的整数。
于是老虎灵机一动,想出了一个绝妙的主意:
如果答案是一个分数 x/y 的话,如果该分数已经不可约分,那么学生需要在答题卡上填涂 x+y ;
否则,学生需要约分后再填涂。
例如,答案为 4/6 时应该填涂 5 ,而答案是 233/1 时应该填涂 234 。
不过,这样似乎并不能解决所有的问题,虽然答案是 1996/2 时可以填涂 999 ,
但答案是 1998/2 时应当填涂 1000 ,但答题卡已经不允许填涂这么大的数字了。
因此,在出题时一定要谨慎,避免这种情况的发生。
一天,老虎和蒜头聊起这事。于是蒜头产生了一个奇怪的问题:给定 A≤B,C≤D ,那么对于所有 A≤x≤B,C≤y≤DA≤x≤B,C≤y≤D ,有多少分数 x/y 是可以被填涂进答题卡的?你现在的任务就是解决这样的问题。
输入格式
共一行四个正整数 A,B,C,DA,B,C,D 。
输出格式
共一行一个整数,表示可以被填涂进答题卡的数字个数。由于答案可能很大,你只要输出其对 109+7 取模的答案即可。
样例
样例一
input
5 8 3 6
output
16
样例二
input
2018 2019 2018 2019
output
2
数据范围及限制
对于 100% 的数据,1≤A,B,C,D≤1012,A≤B,C≤D1≤A,B,C,D≤1012,A≤B,C≤D 。
对于 20% 的数据,1≤A,B,C,D≤101≤A,B,C,D≤10 。
对于 40% 的数据,1≤A,B,C,D≤10001≤A,B,C,D≤1000 。
对于 60% 的数据,1≤(B−A+1)×(D−C+1)≤1061≤(B−A+1)×(D−C+1)≤106 。
对于另外 20% 的数据,1≤A,B,C,D≤10^6 1≤A,B,C,D≤10^6 。
时间限制: 1s
空间限制: 512MB
思路:
枚举最简分数i/j,求出一个k的范围,使得 A<=ik<=B,C<=jk<=D。
代码:
//by juruoZAR
#pragma GCC optimize(2)
#include<bits/stdc++.h>
using namespace std;
#define maxn 20000
#define read(x) scanf("%d",&x);
int n;
char a[maxn+5];
char b[maxn+5];
int ans[maxn+5],cnt=0;
int main() {
scanf("%s",a+1);
n=strlen(a+1);
for(int i=1;i<=n;i++) b[i]=a[i];
for(int i=1;i<n;i++) {
int s=0;
for(int j=1;j<=n;j++) a[j]=b[j];
for(int j=1;j<=n-i;j++) {
if(a[j]!=a[j+i]) {
s++;
if(s>=2) goto END;
a[j+i]=a[j];
}
}
ans[++cnt]=i;
continue;
END:;
s=0;
for(int j=1;j<=n;j++) a[j]=b[j];
for(int j=1+i;j<=n-i;j++) {
if(a[j]!=a[j+i]) {
goto FLAG;
}
}
for(int j=1;j<=i;j++) {
if(j+i>n) break;
if(a[j]!=a[j+i]) s++;
}
if(s<2) ans[++cnt]=i;
FLAG:;
}
int c=(cnt-1)/(1e5)+1;
printf("%d\n",cnt);
for(int i=1;i*c<=cnt;i++) printf("%d ",ans[i]);
return 0;
}