ArcGIS 10.0 Python Script

I need a custom python script made for ArcGIS 10.0. Prior experience programming for ArcGIS is necessary. The script must read a point dataset of events that is stored in a file geodatabase. The sample data has four attributes of importance: date/time field, military grid reference system (MGRS) location, lat field, long field (both in DD). Typical data sets will have approximately 100,000 points. The script must calculate "areas of interest." Areas of interest are defined as X number of events within Y kilometers squared within Z days. For example, I would want the script to identify an area of interest as 1 square kilometer that has seen 3 or more events within a 30 day period. The 30 day period is not necessarily the last thirty days, but any thirty day period in the dataset. After identifying all AOIs in the dataset, the script must create a scaled raster image depicting the number of events. I.E. a 1KM AOI that is shaded green for 3 events, yellow for 4 events, red for 5 events, etc. I think the best approach for this would be a density, but I am open to suggestions. Finally, I'd like a toolbox made in ArcGIS that will allow the selection of the layer to process, ability to input X, Y, & Z values, and name of output raster. The systems this script will operate on have ArcGIS 10.0 Desktop installed. The final deliverable must only use Python for processing. No other dependencies are allowed on the system (PostGIS, etc).