题目描述

喜欢焦叔叔的小妹妹实在是太多了，以至于焦叔叔不得不给她们编号1-N。焦叔叔对于编号为质数的小妹妹有种特别的好感，由于有好感的小妹妹数目太多，焦叔叔自己都数不清楚了，请你来帮忙数一下。

输入格式

多组数据，每行一个正整数N，N不大于1000000000。

数据组数小于15。

输出格式

每行一个数字，有好感的小妹妹数。

输入样例

5

输出样例

3

有技巧的打表，分成根号n块，先找到在哪块里，把之前块的素数和用打表结果弄出来，再爆搜块内有多少小于它的素数。。

（当时一直以为有意想不到的公式，还是太年轻。。）

-----------------------------------------------------------------------------

abcdabcd

5

最长重复子串（后缀树的）+1,10的6次可以，直接上模板

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <algorithm>
#define maxn 1000010using namespace std;string s;
int t;
char inps[maxn];void radix(int *str,int *a,int *b,int n,int m)
{static int count[maxn];memset(count,0,sizeof(count));for(int i=0;i<n;++i) ++count[str[a[i] ] ];for(int i=1;i<=m;++i) count[i]+=count[i-1];for(int i=n-1;i>=0;--i) b[--count[str[a[i] ] ] ]=a[i];
}void suffix_array(int *str,int *sa,int n,int m)
{static int rank[maxn],a[maxn],b[maxn];for(int i=0;i<n;++i) rank[i]=i;radix(str,rank,sa,n,m);rank[sa[0] ]=0;for(int i=1;i<n;++i) rank[sa[i] ]=rank[sa[i-1] ]+(str[sa[i] ]!=str[sa[i-1] ]);for(int i=0;1<<i <n;++i){for(int j=0;j<n;++j){a[j]=rank[j]+1;b[j]=j+(1<<i)>=n?0:rank[j+(1<<i) ]+1;sa[j]=j;}radix(b,sa,rank,n,n);radix(a,rank,sa,n,n);rank[sa[0] ]=0;for(int j=1;j<n;++j){rank[sa[j] ]=rank[sa[j-1] ]+(a[sa[j-1] ]!=a[sa[j] ] || b[sa[j-1] ]!=b[sa[j] ]);}}
}
int substr(string str)
{string rev;static int s[maxn],sa[maxn],rank[maxn],h[maxn];int n=str.length();copy(str.begin(),str.end(),s);suffix_array(s,sa,n,256);for(int i=0;i<n;i++) rank[sa[i]]=i;int k=0;int ans1=0,pos1=0;for(int i=0;i<n;++i){k= k==0?0:k-1;while(rank[i]>0 && s[i+k]==s[sa[rank[i]-1]+k ]) ++k;h[rank[i] ]=k;if(h[rank[i] ]>ans1){ans1=h[rank[i] ];pos1=i;}}return ans1;
}
int main()
{scanf("%s",inps);s.assign(inps);printf("%d\n",substr(s)+1);return 0;
}

---------------------------------------------------------------------------------

1
2

1.000000
3.000000

概率题，最简单的想法，假设需要n种蘑菇，当前已踩到i种，则下次采有(n-i)/n的可能踩到不一样的，取个倒数，也就是采n/(n-i)次能踩到1种不一样的

把i从1到n的值加起来就好了~

（错位相减什么的推出来的，我表示高中那玩意儿就没怎么算对过，实在不敢用）

#include <iostream>
#include <cstdio>using namespace std;int n;
double ans;int main()
{while(~scanf("%d",&n)){ans=0;for(int i=1;i<=n;i++) ans+=((double)n)/((double)i);printf("%.6lf\n",ans);}return 0;
}

-----------------------------------------------------

3 3 30 20 1
0 0
0 50
50 0
50 50
0 1000
1000 0

91.500000

二分匹配问题，总以为求最优的要用km算法是，丁神于是推荐了poj 3343这个神题，明白了要先二分时间，再判断在时间k内的最大匹配，匈牙利一下，看是否等于敌人数

但是，如果炮比敌人少怎么办呢？菊苣教育我们，把炮拆点，打的每炮拆成一个点，共n*m个点，就好了~

-----------------------------------------------------------------

5 0 0 1 1
SESNW

4

模拟就可以，如果顺风，就动，如果逆风就扔锚停住

（开始没看懂题以为那个（x，y+1）只是代表方向向量，wa了后换成了每次移动一格的理解方式，就过了。。）

#include <iostream>
#include <cstdio>
#include <cmath>
#define maxt 100005using namespace std;char input[maxt];
int t,sx,sy,ex,ey,ans;int counting(char dir1,char dir2)
{int flag1=0,flag2=0;int i;for(i=0;i<t;i++){if(input[i]==dir1) flag1++;if(input[i]==dir2) flag2++;if(flag1==abs(sy-ey) && flag2==abs(sx-ex)) break;}if(i==t) return -1;else return i+1;
}int main()
{while(~scanf("%d%d%d%d%d",&t,&sx,&sy,&ex,&ey)){scanf("%s",input);if(ex>sx){if(ey>sy){ans=counting('N','E');}else{if(ey==sy){ans=counting('0','E');}else{ans=counting('S','E');}}}else{if(ex==sx){if(ey>sy){ans=counting('N','0');}else{if(ey==sy){printf("0\n");continue;}else{ans=counting('S','0');}}}else{if(ey>sy){ans=counting('N','W');}else{if(ey==sy){ans=counting('0','W');}else{ans=counting('S','W');}}}}printf("%d\n",ans);}return 0;
}